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Introduction 


The  principal  objective  of  this  proposal  is  to  obtain  physico-chemical  information  on  estrogen 
derivatives  and  correlate  this  knowledge  to  their  known  biological  functionality.  Slight  chemical 
variation  in  these  molecules  can  change  the  carcinostatic  potentials  from  agonistic  to  inhibitory. 
Relating  these  biological  reactions  to  physical  properties  such  as  point  charges  of  atoms  and  the 
electrostatic  potential  is  a  logical  first  step  in  the  intelligent  design  of  therapeutic  drugs  [1-4].  We 
are  obtaining  information  about  the  electronic  properties  of  estrogen  derivatives  from 
experimental  determination  of  their  electron  density  using  high  quality  single  crystal  X-ray 
crystallography  at  low  temperatures. 


Body 

Task  1.  Preliminary  studies  on  a  series  of  crystals  of  estrogen  derivatives 

•  Feasibility  study  of  already  available  crystals  of  estrogen  derivatives 

Feasibility  studies  on  crystals  of  6  commercial  estrogen  derivatives  (17a-estradiol,  17^- 
estradiol,  4,3,17-estriol,  2,3,17-estriol,  3,16a,  17P-estriol  and  estrone)  demonstrated  their 
unsuitability  for  charge  density  studies.  These  studies  require  high  quality  crystals  of  at 
least  0.015  mm^.  All  the  crystals  under  examination  proved  to  be  either  too  small  or  not 
of  a  high  enough  quality.  The  2,3,17-estriol  crystals  also  appeared  to  be  air  and  light 
sensitive. 


•  Development  of  crystallization  methods  for  the  derivatives  not  vet  available  as  high 
quality  single  crystals 

We  have  developed  crystallization  methods  for  obtaining  high  quality  crystals  of  several 
different  estrogen  derivatives.  We  succeeded  in  the  preparation  of  high  quality  crystals  of 


17a-estradiol»l/2  HjO 

17a-estradiol 

4,3,17-estriol 

3,16a,  np-estriol 

17p-estradiol«l/2H20 

17P-estradiol»urea 

17P-estradiol»  %  MeOH*  ‘/3H2O 

17p-estradiol*i/2MeOH 

estrone 


(crystallized  from  methanol) 
(crystallized  from  methanol) 
(crystallized  from  acetonitrile) 
(crystallized  from  acetonitrile) 
(crystallized  from  isopropanol) 
(crystallized  from  benzene+urea) 
(crystallized  from  methanol+water) 
(crystallized  from  methanol) 
(crystallized  from  ethanol+ethyl  acetate) 


For  the  17/?-estradiol  or  17flr-estradiol,  the  best  conditions  found  involve  dissolving  the 
powder  in  “wet”  methanol.  Allowing  the  solvent  to  evaporate  as  slowly  as  possible  (~2-4 
weeks)  will  yield  charge  density  quality  crystals.  This  often  results  in  the  formation  of 
different  crystal  systems  (such  as  17y^estradiol  •  V2  MeOH  or  17y5-estradiol  •  %  MeOH  • 
H2O).  This  is  probably  due  to  the  extremely  small  difference  in  packing  energy  of  the 
different  systems.  On  the  contrary,  the  best  estriol  crystals  were  obtained  from 
acetonitrile,  and  the  best  estrone  crystals  -  from  ethanol:  ethyl  acetate  system. 
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•  Temperature  studies  on  each  sample  to  define  the  appropriate  temperature  for  the 
measurement 


Variable  temperature  studies  were  performed  on  the  diffractometer  on  all  of  the  above 
mentioned  crystals.  As  our  criterion  for  stability  we  chose  the  change  in  mosaicity  of  the 
respective  crystal  as  identified  by  the  broadening  of  the  profiles  of  the  diffraction 
intensities  in  comparison  to  a  room  temperature  measurement.  In  our  in-house 
experiments  using  Ag-  (X;=0.5609  A)  and  Mo-  (X.=0.71073  A)  radiation,  all  of  the 
samples  proved  to  be  stable  under  nitrogen  cooling  to  90  K. 

Further  studies  with  liquid  helium  cooling  (ca.  17  K)  were  performed  with  synchrotron 
radiation  at  the  Argonnne  (A.=0.40663  A)  and  Brookhaven  (^0.643  A)  National 
Laboratories.  Surprisingly,  these  studies  showed  deviating  results:  while  investigations 
on  estrone  at  Argonne  National  Laboratory  proved  its  suitability  for  He-temperature 
charge  density  studies,  crystals  of  17a-estradiol»  V2U2O,  17|3-estradiol*urea,  17P- 
estradioNViMeOH,  17P-estradiol*y3MeOH»14H20  and  estrone,  studied  at  Brookhaven 
National  Lab  showed  a  vast  increase  in  their  mosaicity  after  cooling.  Modifying  the 
crystal  mounting  technique  and  extending  the  time  period  during  which  a  crystal  was 
cooled  to  the  final  temperature  of  ca.  17  K  proved  to  be  successful  for  a  crystal  of 
3,16a,17P-estriol. 

Several  estrogen  crystals  were  tested  in-house  with  Ag-  radiation  using  helium  cooling. 
Contrary  to  our  previous  synchrotron  results,  we  found  no  significant  increase  in  the 
crystals'  mosaicities  after  slow  cooling  using  the  modified  mounts. 

These  experiments  have  shown  a  vast  increase  of  scattering  power  of  the  respective 
crystals  at  He-temperature  compared  to  the  room  temperature  diffraction  experiment.  A 
typical  example  is  shown  in  the  Appendix  A;  the  two  pictures  compare  the  scattering 
ability  of  the  same  crystal  of  17  a-estradiol«l/2  H2O  at  room  temperature  and  at  17  K. 

As  a  conclusion  of  these  studies,  it  is  preferable  to  use  He  temperatures  for  the  charge 
density  experiments,  assuring  that  the  crystal  quality  remains  good  after  cooling.  Use  of 
synchrotron  radiation  is  also  preferable  since  it  gives  a  big  increase  in  the  scattering 
power.  Finally,  a  high  quality  charge  density  data  set  has  been  successfully  measured  at 
the  Advanced  Photon  Source  at  Argonne  National  Lab  at  about  18  K  on  the  estrone 
crystal. 


•  Routine  X-rav  crystal  structure  determination  on  uncharacterized  derivatives 

In  the  course  of  our  crystallization  efforts  we  crystallized  and  determined  three  new 
crystal  modifications  of  estradiol  derivatives: 

173-estradiol«2/3  MeOH«l/3  H2O  (Appendix  B),  reported  in  Acta  Cryst.  (2003) 

-  17P-estradiol«l/2MeOH  (Appendix  H),  reported  in  Acta  Cryst.  (1999) 

-  1 7a-estradiol»  1/2  H2O  (Appendix  D). 
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Task!.  Electron  density  studies  on  the  estrogen  analogs 


•  Electron  density  studies  at  liquid  nitrogen  temperatures 

We  have  collected  complete  data  sets  for  charge  density  analyses  on  crystals  of 

-  17a-estradiol»‘/2H20 

-  17a-estradiol 

-  3,16a,17p-estriol 

-  17P-estradiol»‘/2H20 

-  17P-estradiol*urea 

-  17P-estradiol*Vi  MeOH 

-  estrone 

The  first  measured  data  set  for  the  17P-estradiol«l/2H20  crystal  appeared  to  be  of 
insufficient  quality  due  to  several  instrumentation  problems.  This  data  set  had  to  be  re¬ 
measured.  Two  data  sets  for  the  17p-estradiol»urea  crystal  were  collected  due  to 
reflection  overlapping  problems  in  the  first  data  set.  The  other  five  data  sets  were  of  the 
high  quality  necessary  for  electron  density  analysis. 

To  develop  a  procedure,  how  to  get  fast  (in  several  hours)  high  quality  charge  density 
data  for  estrogens,  the  estrone  and  17P-estradiol»urea  crystals  have  also  been  measured  at 
the  Advanced  Photon  Source  at  Argonne  National  Laboratory.  The  experiments  were 
unsuccessful  due  to  an  instrumentation  problem  that  was  discovered  after  the 
experiments. 


•  Electron  density  studies  at  helium  temperatures 

A  complete  and  high  quality  charge  density  measurement  has  been  performed  on  the 
estrone  crystal  at  the  Argonne  National  Lab  at  about  18  K.  The  data  analysis  is  currently 
proceeding  but  incomplete. 

•  Analysis  of  the  experimental  data,  comparison  of  the  results,  preparation  of  manuscripts 

After  processing  several  data  sets  of  estrogen  derivatives,  we  faced  an  unexpected 
problem:  the  outcome  partially  depends  on  the  methodology  of  the  data  treatment. 
Therefore,  we  had  to  develop  a  uniform  approach  to  treat  the  experimental  data, 
including  the  CCD  data  integration  and  the  following  model  refinement. 

The  integration  of  the  raw  data  involves  integrating  the  intensity  of  the  reflections  as 
measured  by  a  two-dimensional  CCD  detector.  Several  parameters  must  be  defined  to 
determine  exactly  how  the  software  integrates  the  reflections.  It  was  found  that  different 
detector  settings,  even  for  the  same  data  set,  require  different  box  size  parameters  and 
profile  fitting  limits.  The  simple  sum  perimeter  limit  of  0.02  was  found  to  be  the  best 
value  for  all  data  sets. 
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An  optimized  local  atomic  coordinate  system  has  been  developed  (Appendix  C)  in  order 
to  better  compare  the  results  from  studies  on  different  estrogen  analogs  (the  paper  is  in 
press  in  the  Journal  of  Applied  Crystallography). 

It  was  also  found  that  the  starting  values  for  the  multipole  model  of  the  electron  density 
greatly  influenced  the  path  the  refinement  would  take.  It  was  determined  that  a  specific 
set  of  starting  values  should  be  applied  to  each  structure  to  ensure  consistency.  These 
optimum  values  were  determined  after  preliminary  refinements  of  several  estrogen 
derivatives;  they  are  shown  in  the  Appendix  C.  The  model  parameters  in  the  least- 
squares  process  have  to  be  refined  in  small  groups  in  order  to  avoid  severe  correlation 
problems  associated  with  the  relatively  large  size  of  the  estrogen  molecules. 

All  this  work  was  necessary  in  order  to  associate  the  differences  in  the  electron  density 
and  electrostatic  potentieil  in  estrogens  with  their  chemical  behavior,  and  not  with 
differences  in  the  data  treatment. 

-  17a-estradiol«‘/2H20:  As  it  was  discovered  during  the  electron  density  model 
refinement,  the  hydrogen  atom  in  the  water  molecule  was  disordered.  This 
complicated  the  refinement  significantly,  however  it  was  successfully  completed 
as  well  as  the  full  topological  analysis.  The  electrostatic  potential  has  also  been 
calculated.  The  complete  results  are  in  Appendix  D.  We  are  in  the  process  of 
final  data  analysis  and  preparing  the  manuscript.  The  results  were  reported  at  the 
DOD  Era  of  Hope  Meeting  (2002). 

17a-estradiol:  The  charge  density  study  has  been  completed  (Appendix  E).  We 
are  in  the  process  of  final  data  analysis.  The  results  were  reported  at  the  DOD  Era 
of  Hope  Meeting  (2002). 

-  3,16a,17P-estriol:  A  significant  effort  was  put  in  the  data  re-refinement  during  the 
methodology  development.  The  new  methodology  resulted  in  the  significant 
improvement  of  the  model.  The  final  results  are  reported  in  Appendix  F.  We  are 
in  the  process  of  final  data  analysis  and  preparing  the  manuscript.  The  results 
were  reported  at  the  Gordon  research  conference  (2001),  Pittsburgh  Diffraction 
conference  (2001),  Broker  Users  Group  Meeting  (2001)  and  the  DOD  Era  of 
Hope  Meeting  (2002). 

17P-estradiol»*/2H20:  The  experimental  data  have  been  successfully  integrated, 
and  we  are  currently  refining  the  electron  density  model. 

17P-estradiol»urea:  The  second  measured  data  set  appeared  to  be  of  sufficient 
quality.  This  allowed  the  multipole  refinement  to  be  completed  as  well  as  a  full 
topological  analysis  of  the  electron  density  and  the  electrostatic  potential  analysis. 
The  manuscript  (Appendix  G)  is  currently  under  revision.  The  results  were 
reported  at  the  DOD  Era  of  Hope  Meeting  (2002). 

17P-estradiol»‘/2MeOH:  This  system  was  particular  challenging,  because  the 
crystal  has  the  lowest  possible  symmetry,  resulting  in  low  redundancy  of  the  data. 
However,  the  multipole  refinement  was  successfully  completed  as  well  as  a  full 
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topological  analysis  of  the  electron  density  and  the  electrostatic  potential  analysis. 
The  complete  results  are  in  Appendix  H.  We  are  in  the  process  of  final  data 
analysis  and  preparing  the  manuscript.  The  results  were  reported  at  the  DOD  Era 
of  Hope  Meeting  (2002). 

Estrone:  Based  on  our  new  developed  methodology,  a  significant  amount  of  work 
has  been  done  in  order  to  extract  more  information  firom  the  data.  The  experiment 
has  been  re-integrated  several  times  in  order  to  get  the  most  precise  structure 
factor  list  with  minimum  noise.  Then,  the  model  has  been  re-refined  according  to 
the  new  approach  (see  above).  We  are  in  the  process  of  finishing  the  manuscript 
(Appendix  I).  The  results  were  reported  at  the  ACA  meeting  (2000),  lUCR 
meeting  (2000),  Pittsburgh  Diffraction  conference  (2001),  Bruker  Users  Group 
Meeting  (2001)  and  the  DOD  Era  of  Hope  Meeting  (2002). 


The  comparison  of  the  results  from  the  series  of  estrogen  analogs  showed  that  the  core 
structure  of  the  different  estrogen  molecules  remains  relatively  unchanged  from  system  to 
system.  As  the  same  time,  significant  differences  are  found  at  the  activity-sensitive 
molecular  parts. 

Accurate  experimental  data  and  the  new  methodology  developed  (see  above)  allowed  us 
to  reveal  such  fine  details  as  oxygen  lone  pair  charge  density  concentrations.  It  appeared 
that  each  hydroxyl  oxygen  atom  had  two  lone  pairs  in  approximately  sp^  type  geometry. 
This  demonstrates  that  these  lone  pairs  are  surprisingly  robust,  and  they  do  not 
significantly  change  in  different  hydrogen  bonding  schemes  or/and  when  the  oxygen  is 
bound  to  the  aromatic  system. 

The  critical  point  analysis  of  covalent  bonds  of  the  estrogen  molecules  demonstrated  the 
characteristic  (3,-1)  critical  points,  with  the  expected  electron  density,  Laplacian  and 
ellipticity  values  given  the  atom  hybridization  and  connectivity.  The  average  electron 
density  values  at  the  critical  points  are  1.91  eA'^  for  the  C  -  O  bond,  2.17  eA'^  for  Car  - 
Car,  1.66  eA'^  for  C  -  C,  1.91  eA'^  for  C  -  H  and  2.07  eA'^  for' the  O  -  H  bond  based  on 
the  17yff-estradiol«urea  study.  In  all  compounds,  the  critical  points  are  well  centered  in  the 
homonuclear  C  -  C  bonds,  and  proportionally  displaced  away  from  the  more 
electronegative  atom  in  the  heteronuclear  bonds. 

Although  the  influence  of  the  lattice  is  inherent  in  the  multipole  model  of  the  estrogen 
molecule,  the  electrostatic  potential  is  believed  to  provide  a  reasonable  representation  of 
a  single  molecule  isolated  from  the  crystal  lattice.  In  these  studies,  the  electrostatic 
potentials  of  the  hydroxyl  groups  on  the  aromatic  rings  were  qualitatively  consistent  in 
both  shape  and  magnitude.  Generally  covering  nearly  the  entire  lone  pair  region,  the 
maximum  magnitudes  are  ~  -20  eA  *.  The  C17  hydroxyl  groups  are  also  relatively 
consistent  in  their  lone  pair  regions  and  they  exhibit  an  interesting  trend.  The  potential 
changes  significantly  in  both  shape  and  magnitude  based  on  the  relative  position  of  the 
hydroxyl  hydrogen  with  respect  to  H17.  17>9-estradiol»urea  and  17£ir-estradiol«  l^HaO 
displayed  a  shift  in  the  potential  toward  the  lone  pair  region  opposite  H17,  as  well  as 
increased  magnitudes  of  — 35  and  — 45  eA  *.  H20  and  H20’  were  virtually  opposite 
H17  in  17yff-estradiol»  l^MeOH.  This  resulted  in  a  fairly  uniform  distribution  of  the 
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electrostatic  potential  with  a  maximum  value  of  — 20  eA'*.  Plots  of  electrostatic 
potential  isosurfaces  are  shown  in  Appendices  J-K. 


Tasks  3-6.  Extension  of  the  above  mentioned  studies  to  the  series  of  E?  C-ring  analogues  and 
nonsteroidal  molecules. 


•  Feasibility  study  of  already  available  crystals 

Feasibility  studies  on  commercial  diethylstilbestrol  crystals  demonstrated  their 
unsuitability  for  charge  density  studies.  These  studies  require  high  quality  crystals  of  at 
least  0.015  mm^.  All  the  crystals  under  examination  appeared  to  be  too  small. 

•  Development  of  crystallization  methods  for  the  derivatives  not  vet  available  as  high 
quality  single  crystals 

A  significant  effort  was  made  to  obtain  crystals  of  diethylstilbestrol.  This  is  a  molecule 
known  to  interact  strongly  with  the  estrogen  receptor.  It  crystallizes  in  a  centrosymmetric 
space  group  with  the  molecule  sitting  on  a  special  position,  which  makes  it  a  good 
candidate  for  charge  density  studies.  Unfortunately,  slow  evaporation  studies  and  solvent 
diffusion  studies  using  many  solvents  such  as,  isooctane,  p-chlorophenol,  acetonitrile, 
acetone  and  several  different  types  of  alcohols  were  not  successful.  Only  a  powdery 
residue  was  observed,  or  in  the  best  case,  microcrystals. 


Task  7.  Final  analysis 


Presumably  there  is  a  preferred  orientation  for  the  hydrogen  of  the  C 17  hydroxyl  group  in 
the  active  site  of  the  estrogen  receptor.  Unfortunately,  the  protein  crystal  structures 
determined  to  date  do  not  have  enough  resolution  to  determine  this.  It  seems  possible  that 
this  is  one  mechanism  by  which  the  estrogen  ligands  could  stimulate  different  responses 
from  the  estrogen  receptor.  Binding  of  different  ligands  would  most  likely  result  in 
different  orientations  of  the  C17  hydroxyl  hydrogen,  thereby  presenting  different 
electrostatic  potentials  to  the  receptor.  This  could  help  to  account  for  the  changes  in 
binding  affinity  and  activity  of  the  different  estrogen  derivatives.  Continued  effort  must 
be  made  to  analyze  more  systems  because  a  sample  of  six  structures  is  not  enough  to 
confirm  the  trends  that  are  appearing  in  the  electron  density  and  the  electrostatic 
potential. 


Key  Research  Accomplishments 

♦  Crystallization  of  high  quality  crystals;  the  best  conditions  to  obtain  high  quality 
single  crystals 

♦  Methodology  development  of  the  X-ray  CCD  data  treatment  and  model  least- 
squares  refinement: 
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a)  optimized  temperature  conditions 

b)  best  radiation  type  determination 

c)  best  data  integration  options 

d)  best  local  atomic  coordinate  system 

e)  best  first  approximation  of  the  model 

f)  best  and  uniform  refinement  procedure 

♦  New  crystal  structures  solved  for  17P-estradiol*2/3  MeOH«l/3  H2O;  17^- 
estradiol«l/2MeOH;  17a-estradiol*l/2  H2O 

♦  Charge  density  and  electrostatic  potential  analysis  of  the  17a-estradiol*‘/2H20; 
17a-estradiol;  3,16a,17P-estriol;  17P-estradiol»urea;  17P-estradiol«‘/2  MeOH 
and  estrone 

♦  The  core  estrogen  structure  is  very  consistent  among  different  derivatives 

♦  Parameters  describing  the  activity-sensitive  molecular  parts  are  different 

♦  It  is  possible  to  locate  lone  pair  densities  of  oxygen  in  such  large  systems 

♦  The  deformation  electron  density  distributions  of  all  hydroxyl  oxygen  atoms  are 
close  to  sp^  in  shape,  even  when  oxygens  are  bound  to  aromatic  neighbors 

♦  Hydroxyl  oxygen  lone  pairs  appear  to  be  unaffected  by  completely  different 
hydrogen  bonding  environments 

♦  The  electrostatic  potentials  around  the  oxygen  atoms  are  consistent  in  the 
different  hydrogen  bonding  environments 

♦  The  electrostatic  potential  around  the  C 17  hydroxyl  groups  changes  significantly 
in  both  shape  and  magnitude  based  on  the  relative  position  of  the  hydroxyl 
hydrogen  with  respect  to  H17 


Reportable  Outcomes 
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2.  Pre-doctoral  Traineeship  Award  granted  by  the  Department  of  the  Army  for  Damon 
Parrish.  $66,000  “Measurement  of  the  Electron  Density  Distribution  of  Estrogens  -  a  First 
Step  to  Advanced  Drug  Design” 

3.  Parrish,  D.,  Wu,  N.  &  Pinkerton,  A.  A.,  Charge  Density  Distribution  of  Estrone,  American 
Crystallographic  Association  Meeting,  St.  Paul,  MN,  July  22-27, 2000. 

4.  Parrish,  D.,  Wu,  N.,  Zhurova,  E.A.  &  Pinkerton,  A.  A.,  Preliminary  Results  for  a  Charge 
Density  Study  of  Estrone,  19th  European"  Crystallographic  Meeting,  Nancy,  France, 

August  25-31, 2000. 
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10.  Parrish,  D.  &  Pinkerton,  A.A.,  A  New  Estra-l,3,5(10)-triene-3,17fi-diol  Solvate: 
Estradiol-  Methanol-  Water  (3/2/1 ),  Acta  Cryst.,  2003,  C59,  o80-o82. 

11.  Parrish,  D.,  Comparative  Charge  Density  Studies  of  17a-Estradiol  and  17 fi-Estradiol,  A 
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12.  Chen,  Yu-Sh.,  Kirschbaum,  K.,  Kumaradhas,  P.,  Parrish,  D.,  Pinkerton,  A.  A.  &  Zhurova, 
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Conclusion 


The  ability  to  determine  the  experimental  electron  density  distribution  and  electrostatic 
potential  of  larger  molecular  systems  such  as  estrogens  has  been  demonstrated.  We  have 
developed  the  methodology  of  crystallization,  the  X-ray  CCD  data  treatment  and  least- 
squares  model  refinement  procedure  for  estrogen  crystals  in  order  to  extract  maximum 
reliable  and  comparable  information  from  the  data.  We  completed  charge  density  studies 
of  six  estrogen  derivatives,  performed  the  preliminary  analysis  and  compared  the  results. 
We  found  that  the  deformation  electron  density  distributions  of  all  hydroxyl  oxygen 
atoms  are  near  sp^  in  shape;  their  lone  pair  densities  have  been  reliably  located.  These 
configurations  as  well  as  the  electrostatic  potentials  around  the  oxygen  atoms  are  very 
consistent  in  different  hydrogen  bonding  environments.  The  core  estrogen  structure  is 
also  very  consistent  between  the  derivatives.  As  the  same  time,  the  significant  differences 
are  found  at  the  activity-sensitive  molecular  parts. 

To  be  able  to  reliably  associate  the  molecular  recognition  of  a  drug  molecule  to  receptor 
with  a  electrostatic  potential  distribution  in  the  molecule,  a  significant  data  base  of  charge 
density  for  estrogen  molecules  has  to  be  built.  So  far,  we  have  completed  six  structures, 
and  this  small  sample  is  not  enough  to  make  any  statistical  conclusions.  More  charge 
density  studies  of  different  estrogen  analogues  are  needed  to  complete  this  study. 
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Appendix  A. 


Scattering  of  a  crystal  of  17a-estradiol*l/2  H2O  at  room  temperature  (a) 
and  He  (ca.l7  K)  temperature  (b). 
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The  title  solvate  of  the  steroid  17jS-estradiol  (E2)  with 
methanol  and  water,  Ci8H2402*0.67CH40*0.33H20,  is  the 
first  E2  derivative  to  contain  three  crystallographically 
independent  molecules  in  the  asymmetric  unit.  The  three 
steroid  molecules,  along  with  two  methanol  molecules  and  a 
water  molecule,  create  a  three-dimensional  hydrogen-bonded 
system.  Three-sided  columns  are  formed,  with  the  estradiol 
molecules  aligned  lengthwise  parallel  to  (101),  and  joined  by 
solvent  molecules  at  both  hydrophilic  ends.  The  three 
estradiol  molecules  differ  slightly  in  their  ring-bowing 
angles,  i.e.  the  angle  between  the  mean  plane  of  the  A  ring 
and  that  of  the  BCD  ring;  this  angle  ranges  from  7.1  to  12.2°. 

Comment 

17jS-Estradiol,  E2,  is  a  member  of  the  estrogen  family  of 
hormones.  In  recent  years,  interest  in  these  molecules  has 
focused  primarily  on  understanding  their  biological  role  in 
initiating  breast  cancer.  It  is  well  known  that  their  ability  to 
form  hydrogen  bonds  in  the  active  site  of  the  estrogen 
receptor  (ER)  influences  biological  activity.  This  ability  to 
form  hydrogen  bonds  has  been  clearly  demonstrated  in  an 
array  of  crystal  structures,  especially  in  that  of  E2,  containing 
different  solvent  molecules  or  other  hydrogen-bond  acceptors. 

When  comparing  these  structures,  the  flexibility  of  the 
hydrophobic  region  of  the  molecule,  more  specifically  the  B 
ring,  becomes  very  apparent.  Several  papers  have  already 
discussed  the  ring  bowing  of  this  molecule  (Cooper  et  al,  1969; 
Cody  et  al,  1971;  Busetta  et  al,  1976;  Duax  et  al,  1979).  Weise 
&  Brooks  (1994)  took  it  a  step  further,  performing  molecular- 
modeling  calculations  on  observed  and  novel  confirmations. 
Ivanov  and  co-workers  carried  out  a  related  study  which 
included  a  larger  body  of  compounds  and  more  detailed 
analysis  (Ivanov  et  al,  1998).  They  concluded  that  there  are 
two  possible  conformations  very  close  in  energy,  which  differ 
in  the  B-ring  arrangement,  although  the  chance  of  the 
strained-geometry  binding  is  low,  as  one  quarter  of  the  binding 
energy  (“11.9  kcal  mol“^;  IkcalmoP^  =  4.184  kJ  mol“^)  is 


predicted  to  be  lost  (Anstead  et  al,  1997).  The  point  initially 
postulated  by  Weise  &  Brooks  remains,  i.e.  that  the  flexibility 
to  allow  bending  or  other  conformational  changes  of  the 
ligand  in  the  receptor  appears  energetically  achievable,  and 
this  could  be  an  important  property  in  determining  their 
activity. 

Molecular-dynamics  studies  on  the  ligand-binding  domain 
(LBD)  of  the  ER  demonstrate  that  the  motion  of  the  LBD 
requires  the  A  ring  to  remain  fairly  steady  while  the  CD  ring 
retains  a  higher  degree  of  freedom.  The  range  of  motion  found 
in  that  study  (Maalouf  et  al,  1998)  agrees  very  well  with  the 
reported  crystal  structures.  It  is  well  known  that  the  hydrogen 
bonding  of  the  ligand  in  the  LBD  causes  conformational 
changes  in  the  receptor.  This  seemingly  accommodating 
motion  of  the  LBD  and  the  ligand  supports  the  idea  postu¬ 
lated  by  Weise  &  Brooks  that  the  flexibility  of  E2,  or  of  any 
other  ligand  which  enters  the  LBD,  could  also  effect  the 
activity  of  the  ligand.  A  more  rigid  or  flexible  ligand  could 
change  the  natural  motion  of  the  ER  complex,  resulting  in  a 
change  in  activation  factor  (AF-2)  activity,  and  therefore 
possibly  effecting  co-activator  recruitment.  This  process  is  not 
well  understood;  the  flexibility  of  the  ligand  is  certainly  only 
one  of  many  physical  factors  associated  in  the  activity  of  E2. 
The  structure  of  the  title  solvate,  (I),  continues  the  trend  in 
observing  the  flexibility  of  this  molecule  from  a  structural  and 
hydrogen-bonding  point  of  view,  while  the  molecule  lies  in  the 
lowest  energy  conformation. 


•0.67CH40-0.33H20 


(0 

The  crystal  of  (I)  contains  three  crystallographically  unique 
E2  molecules  (Fig.  1).  The  B  rings  of  the  E2  molecules  adopt 
the  typical  conformation  of  a  distorted  7Q:,8j6-half-chair,  and 
this  is  responsible  for  most  of  the  structural  flexibility. 
Calculation  of  the  ring-bowing  angle,  as  defined  by  Duax  & 
Norton  (1975),  reveals  a  range  of  5.1°  for  the  three  molecules. 
Table  1  compares  the  ring-bowing  angles  of  the  currently 


Figure  1 

A  view  of  molecule  A  of  (I),  with  the  atom-numbering  scheme;  the  other 
two  molecules  are  similarly  labeled,  with  the  corresponding  suffix  B  or  C 
Displacement  ellipsoids  are  plotted  at  the  50%  probability  level.  The 
solvent  molecules  have  been  omitted  for  clarity. 
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known  E2  crystal  structures.  The  spread  of  over  12°,  with 
angles  found  distributed  over  the  entire  range,  indicate  the 
shallowness  of  the  potential  associated  with  the  bowing 
deformation.  In  a  system  as  dynamic  as  the  human  body,  the 
molecule  would  certainly  have  a  wide  range  of  conformations 
easily  available. 

The  packing  of  (I)  results  in  three-sided  cylinders  of 
estradiol  molecules  arranged  parallel  to  (101),  with  the  Cl 8 
methyl  group  pointing  towards  the  center  (Fig.  2).  This,  of 
course,  aligns  the  large  hydrophobic  regions  of  the  molecules, 
as  well  as  the  hydrophilic  hydroxy  groups.  These  cylinders 
then  stack  on  top  of  each  other,  creating  three-sided  columns 
(Fig.  3)  held  together  by  hydrogen  bonds  involving  the 
hydroxy  groups  and  the  three  solvent  molecules.  A  more 
detailed  picture  of  the  hydrogen  bonding  between  the  solvent 
and  E2  molecules  can  be  seen  in  Fig.  4  and  from  the  data  in 
Table  3.  There  is  also  an  intercolumnar  hydrogen  bond,  which 
occurs  between  the  water  molecule  and  a  3-hydroxy  group. 
This,  along  with  the  intercolumnar  hydrogen-bonded  solvent, 


Figure  2 

The  molecular  arrangement  of  (I),  looking  down  the  three-sided  column. 
Displacement  ellipsoids  are  plotted  at  the  50%  probability  level 
[symmetry  codes:  (i)  x,  y,  z;  (ii)  — 2  —  ;c,  y  —  J,  —  1  —  z;  (iii)  ;c  —  1,  y,  z]. 


Figure  4 

The  hydrogen-bonded  region  of  the  column  shown  in  Fig.  3,  including  the 
solvent  molecules  and  only  the  A  or  D  rings  of  E2.  Displacement 
ellipsoids  are  plotted  at  the  50%  probability  level  [symmetry  codes:  (i)  x, 
y,  z;  (ii)  -x,  y  - -z;  (iii)  -1  -  y  ~  i  ^z\  (iv)  1 4-  jc,  y ,  1  +  z;  (v)  -2  -  x, 
y  -  i  -1  -  z;  (vi)  x,y,l  +  z;  (vii)  -1  -  x,  y,  z]. 
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Figure  5 

A  schematic  diagram  of  the  packed  columns  in  (I)  and  their  relation  to 
the  hydrophilic  layer  containing  the  hydroxy  groups  and  solvent 
molecules. 


creates  a  hydrophilic  layer  which  is  approximately  30°  from 
the  perpendicular  to  the  columns  (Fig.  5).  Structurally,  all 
bond  lengths  and  angles  are  in  expected  ranges  (Table  2). 


Figure  3 

The  packing  of  E2  groups  in  (I),  with  the  hydrogen-bonded  solvent 
molecules,  viewed  perpendicular  to  the  columns.  Dotted  lines  indicate  the 
region  detailed  in  Fig.  4.  Displacement  ellipsoids  are  plotted  at  the  50% 
probability  level  [symmetry  codes:  (i)  jc,  y ,  z;  (ii)  1  -h  x,  y,  1  -1-  z;  (iii)  —2  —  x, 
y-i-l-  z;(iv)  -1  -  JC.y  -zKv)^:,^,  1  +  z;(vi)  -x,y-\,  -z]. 


Experimental 

Crystals  were  grown  by  slow  evaporation  from  a  methanol  solution 
open  to  the  air.  It  is  presumed  that  the  methanol  was  either  wet  prior 
to  use  in  this  experiment  or  absorbed  atmospheric  moisture.  Crystals 
of  (I)  were  grown  over  the  course  of  3  d  and  harvested  when  reaching 
the  appropriate  size. 

Crystal  data 

Ci8H2402-0.67CH400.33H20 
Mr  =  299.74 
Monoclinic,  Pl^ 
a  =  11.7152  (4)  A 
b  =  19.6270  (6)  A 
c  =  12.1310  (4)  A 
=  117.978  (1)° 

V  =  2463.34  (14) 

Z  =  6 


1.212  Mg  m-^ 

Ag  Ka  radiation 
Cell  parameters  from  7309 
reflections 
e  =  2.3-27.9° 

-  0.05  mm“^ 
r=298  (1)K 
Cuboid,  colorless 
0.4  X  0.3  X  0.3  mm 
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Data  collection 

Siemens  SMART  Platform  CCD 
area-detector  diffractometer 
(o  scans 

Absorption  correction:  empirical 
via  multipole  expansion 
(Blessing,  1995)  using  SADABS 
(Sheldrick,  1996) 

Tnoin  =  0.840,  =  0.985 

Refinement 

Refinement  on 
R[F^  >  2g{F^)]  =  0.058 
wR{F^)  =  0.145 
S  =  0.99 

21 122  reflections 
698  parameters 


44  195  measured  reflections 
21 122  independent  reflections 
15  597  reflections  with  I  >  2a(J) 
=  0.048 

0max  =  27.9'’ 

/^  =  -18  19 
A:  =  -32  29 
/  =  -20  20 


H  atoms  treated  by  a  mixture  of 
independent  and  constrained 
refinement 

w  =  l/[<r^(F/)  +  (0.0812P)2] 
where  P  =  {Fj^  +  2F^)I3 
(A/or)„^,  =  0.001 
A/Omax  =  0.41eA/ 

^Pmin  —0.23  e  A 


Table  1 

Comparison  of  ring-bowing  angles  (°)  in  reported  E2  crystal  structres. 


Compound 

Molecule  1  Molecule  2 

Molecule  3 

angle 

angle 

angle 

17^-estradiol-0.5H2Ot 

15.6 

17^-estradiol-propanolt 

12.9 

17^-estradiol-ureat 

5.6 

17^-estradiol-0.5MeOH$ 

10.4 

3.4 

(i)§ 

12.2 

11.9 

7.1 

t  Wiese  &  Brooks  (1994). 

t  Parrish  &  Pinkerton  (1999),  §  This  work. 

Table  2 

Selected  bond  lengths  (A). 

01A-C3A 

1.3716  (18) 

024-C174 

1.4354  (19) 

01B-C3S 

1.3792  (17) 

02B-C17B 

1.4305  (19) 

oic-ac 

1.3691  (18) 

02C-C17C 

1.4400  (19) 

Table  3 

Hydrogen-bonding  geometry  (A,  ®). 

D-H--A 

D-H 

H*  •  *4  Z>-  •  -4 

D-n-A 

OIA-HIOA  -024^ 

0.74  (3) 

1.89  (3)  2.6226  (17) 

171  (3) 

01J5-H10B--02B‘ 

0.85  (3) 

1.83  (3)  2.6654  (18) 

171  (3) 

OIC-HIOC-  -03 

0.76  (3) 

1.82(3)  2.5767(18) 

170  (3) 

02A-H20A-  •  045 

0.76  (3) 

1.97(3)  2.7017(18) 

162  (2) 

02B-H20B-  •04B” 

0.73  (2) 

2.06  (2)  2.7823  (19) 

168  (2) 

02C-H20C-  •  OlC" 

0.73  (2) 

1.99  (3)  2.7169  (17) 

175  (3) 

03-H30A-G1A*" 

0.72  (4) 

2.02  (4)  2.743  (2) 

178  (4) 

03-H30B  •  OIB'' 

0.76  (3) 

2.04  (3)  2.7926  (19) 

173  (3) 

04A-H40A  •.02C 

0.83  (3) 

1.85  (3)  2.6809  (18) 

174  (3) 

04B-H40S'  •  *044''* 

0.78  (2) 

1.88  (2)  2.6462  (17) 

166  (2) 

Symmetry  codes:  (i)  1  +  jr,  y,  1  +  z;  (ii)  —3—x,\-^y,—2  —  z;  (iii)  jr  -  1,  y,  ^  —  1;  (iv) 
-JC,5+y,  -z;  (v)  1  -\-x,y,  z;  (vi)  j:-l,y,  z. 


The  absolute  configuration  of  the  E2  molecules  in  (I)  was  known 
from  the  natural  product  starting  material.  The  intensity  data  were 
corrected  for  decay  and  absorption  using  SADABS  (Sheldrick,  1996). 
All  H  atoms  of  the  solvent  molecules  and  the  hydroxy  groups  of  the 
E2  molecules  were  located  in  a  difference  map  and  refined  with 
isotropic  displacement  parameters.  The  remaining  H  atoms  were 
included  with  idealized  geometries  (C— H  =  0.96-0.98  A),  and  their 
isotropic  displacement  parameters  were  refined. 

Data  collection:  SMART  (Bruker,  1998);  cell  refinement:  SAINT 
(Bruker,  1999);  data  reduction:  SAINT;  program(s)  used  to  solve 
structure:  SHELXS91  (Sheldrick,  1990);  program(s)  used  to  refine 
structure:  SHELXL91  (Sheldrick,  1997);  molecular  graphics: 
SHELXTL  (Siemens,  1994);  software  used  to  prepare  material  for 
publication:  SHELXTL. 

We  thank  the  College  of  Arts  and  Sciences  of  the  University 
of  Toledo  for  generous  financial  support  of  the  X-ray 
diffraction  facility. 


Supplementary  data  for  this  paper  are  available  from  the  lUCr  electronic 
archives  (Reference:  SQ1003).  Services  for  accessing  these  data  are 
described  at  the  back  of  the  journal. 
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A  Standard  Local  Coordinate  System  for  Multipole  Refinements  of  the 

Estrogen  Core  Structure 
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Synopsis  The  initiation  of  a  comparative  charge  density  study  on  a  series  of  estrogen  derivatives  revealed  the 
need  for  a  standardized  local  coordinate  system  for  each  atom  in  order  to  facilitate  comparison  of  the  multipole 
refinements.  Herein  we  propose  such  a  system  for  the  core  structure  of  a  generic  estrogen  molecule  and  suggest 
starting  parameters  for  the  population  of  the  most  prominent  multipoles. 

Abstract  A  comparative  charge  density  study  on  a  series  of  estrogen  derivatives  has  been  initiated.  The  study 
utilizes  the  Hansen-Coppens  atom  centered  multipole  model  to  describe  the  valence  electron  density  distribution. 
Direct  comparison  of  the  population  parameters  for  each  estrogen  after  the  respective  multipole  refinements 
requires  standardization  of  the  atom  centered  local  coordinate  systems.  Such  a  standard  coordinate  system  for  the 
common  estrogen  core  is  reported,  taking  advantage  of  the  shape  of  those  multipoles  which  have  the  spatial 
characteristics  of  sp^  and  sp^  hybrid  orbitals.  Additionally,  populating  these  principal  multipoles  at  the  beginning 
stage  of  the  refinements  improves  the  stability  of  these  large,  highly  correlated  calculations. 

Keywords:  Estrogen;  charge  density;  coordinate  system 
1.  Introduction 

Recent  advances  in  computers,  software,  and  X-ray  data  collection  technology  have  made  charge  density 
studies  on  a  series  of  larger  "small  molecules"  in  a  reasonable  amount  of  time  feasible.  A  comparative  charge 
density  analysis  has  been  initiated  to  study  a  series  of  estrogen  derivatives. 

Electron  density  refinements  from  experimental  X-ray  diffraction  data  have  been  performed  using  the  XD 
suite  of  programs  (Koritsanszky  et  al.,  1995),  which  utilizes  the  Hansen-Coppens  model  to  describe  the  valence 
electron  density  (Coppens  1997,  Hansen  &  Coppens  1978).  This  model  is  based  on  multipole  expansions  centred 
on  the  atoms  and  concomitantly  requires  a  user  defined  local  Cartesian  coordinate  system  for  each  atom  in  the 
structure . 
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In  a  comparative  study,  consistency  in  this  model  setup  is  essential  to  enable  direct  comparison  of  the 
parameters  of  the  multipole  refinements  for  different  molecules  in  a  related  series.  This  article  proposes  a 
standard  local  coordinate  system  for  the  atoms  of  the  core  structure  of  a  generic  estrogen  molecule. 

As  well  as  being  responsible  for  the  development  of  secondary  sexual  characteristics,  it  has  been 
determined  that  some  estrogens  are  responsible  for  the  initiation  and  progression  of  certain  types  of  breast  cancer. 
This  study  attempts  characterization  of  the  physical  and  electronic  features  of  these  molecules  and  the  relationship 
of  such  submolecular  properties  to  biological  activity. 

Most  biologically  active  estrogens  are  A-ring  or  D-ring  isomers  of  the  naturally  occurring  estradiol  or 
derivatives  thereof,  containing  the  same  hydrocarbon  ring  core  structure.  With  careful  choice  of  directions,  the 
local  coordinate  system  can  take  advantage  of  the  spatial  characteristics  of  certain  multipoles  that,  if  populated, 
would  mimic  the  shape  of  the  bonding  electron  density  derived  from  sp^  and  sp^  hybrid  orbitals.’  (Pichon-Pesme 
et  al.,  1995;  Coppens,  1997)  This  allows  a  majority  of  the  valence  electron  density  to  be  modelled  by  the  linear 
combination  of  only  a  few  multipoles.  Given  the  number  of  parameters  for  such  large,  non-centrosymmetric 
systems,  typically  greater  than  40  atoms,  this  approach  can  considerably  stabilize  the  initial  multipole 
refinements.  By  using  a  standard  coordinate  system,  upon  completion  of  the  multipole  refinement,  a  rational 
comparison  can  then  be  made  among  various  related  charge  density  studies. 

The  17/?-estradiol  molecule,  along  with  the  labelling  scheme  of  the  core  structure,  is  shown  in  scheme  1. 
For  all  methylene  groups,  hydrogen  atoms  on  file  same  side  of  the  eore  as  the  methyl  group  at  position 
18  are  labelled  b,  those  on  the  other  side  are  labelled  a.  The  C  and  D  rings  projected  down  the  C18-C13 
bond  are  provided  as  an  insert  to  define  the  orientation  of  the  methyl  group  hydrogen  atoms. 


Scheme  1 


’  It  is  inportant  to  distinguish  between  the  multipoles,  which  are  density  functions,  and  orbitals,  which  are  derived  from 
wavefunctions. 
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2.  Method 


The  Cartesian  axes  x,y,z  of  the  above  mentioned  atom  centered  coordinate  systems  are  calculated  in  the  program 
XD  (Koritsanszky  et  al.,  1995)  based  on  two  user  defined  vectors.  In  general,  the  first  vector  (vi)  for  a  given 
atom  (1)  originates  at  its  nucleus,  is  directed  towards  another  atom  (2)  in  the  system  -  or,  if  necessary,  a  "dummy" 
atom  -  and  defines  axis(l).  A  second  vector  (V2)  with  the  same  origin  points  towards  another  atom  (3),  or  dummy 
atom,  describing  the  axis(l)  -  axis(2)  plane.  A  third  vector  (V3)  is  taken  perpendicular  to  this  plane, 

Vi  =  (r2  -  ri)  V2  =  (rs  -  rO  V3  =  Vi  x  V2 

where  ri,  r2,  and  rs  are  the  position  vectors  of  atom(l),  atom(2),  and  atom(3),  respectively. 

An  orthonormal  vector  triplet  (Caxi,  ©3x2,  Caxs)  is  then  formed  which  we  have  chosen  to  be  right  handed. 

Caxl  =  Vi/|V,  I  eax2=(V3  X  Vi)/|(V3  X  Vi)|  63x3  =  V3/IV3I 

Although  axes  1, 2  and  3  can  be  assigned  labels  x,y,z  arbitrarily,  for  the  current  molecules  we  define  axis(l)  as  x, 
axis(2)  as  y,  and  axis(3)  as  z  for  non-hydrogen  atoms,  and  axis(l)  as  z,  axis(2)  as  y  and  axis(3)  as  x  for  hydrogen 
atoms. 

Given  these  assignments,  the  valence  electron  density  of  the  aromatic  sp^  type  carbons  atoms  in  the  A- 
ring  (scheme  1),  can  be  initially  approximated  by  populating  two  multipoles,  the  33+  and  20.  Pictorial 
representations  of  the  angular  function  for  these  multipoles  can  be  seen  in  Figure  1.  Defining  the  xy  plane  of  the 
coordinate  system  to  lie  in  the  plane  of  the  bound  atoms,  i.e.  the  aromatic  plane,  will  allow  the  positive  lobe  of  the 
33+  octupole  to  align  with  the  three  covalent  bond  regions.  The  20  quadrupole  will  tiien  be  positioned 
perpendicular  to  this  plane,  a  negative  populations  adding  to  the  in  plane  a  bonding  and  depopulating  the  p  orbital 
as  required  for  tt  bond  formation. 


Figure  1  Coordinate  system  and  angular  functions  of  the  33+  and  20  multipoles  with  hlue  representing  positive  regions 
and  red,  negative.  Functions  visualized  using  MATLAB  (2000). 

The  covalent  bonding  regions  of  the  sp^  type  carbon  atoms  in  the  B,  C,  and  D-rings  of  the  generic 
estrogen  can  effectively  be  modelled  by  a  single  multipole,  the  32-  octupole.  The  axes  associated  with  the  32- 
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octupole  do  not  correspond  to  the  populated  regions,  hence  two  "dummy"  atoms  must  be  used  to  define  the  local 
coordinate  system  in  this  case.  Their  positions  are  calculated  on  the  bisectors  of  the  bond  angles  (the  two-folds 
axes  of  the  tetrahedron)  at  the  atom  of  interest.  An  example  of  an  sp^  coordinate  system  and  the  32-  octupole  is 
shown  in  Figure  2. 

Approximation  of  the  bond  density  of  the  hydrogen  atoms  can  be  achieved  by  initially  populating  a  single 
dipole  along  z.  A  complete  list  of  vector  definitions  for  carbon  and  hydrogen  atoms  for  the  estrogen  core  are  given 
in  Tables  1  and  2.  Using  this  standardized  set  of  coordinates  systems  for  future  theoretical  and  experimental 
electron  density  studies  on  estrogen  derivatives  (and  possibly  expanding  it  for  other  steroids  like  progesterone)  it 
will  be  possible  to  compare  individual  multipole  populations  for  different  molecules  in  this  series. 


Figure  2  Model  coordinate  system  for  a  sp^  type  atom  and  32-  multipole,  with  blue  representing  positive  density  and  red, 
negative.  (MATLAB,  2000). 

Populating  these  principal  multipoles  at  the  beginning  stage  of  the  refinements  may  improve  the  stability 
of  these  large,  highly  correlated  calculations.  We  have  thus  listed  suitable  values  for  the  core  structure  based  on 
our  own  refinements  (Chen  et  al.,  2003)  of  six  different  estrogen  molecules  (Table  3).  We  note  that  the  sum  of 
the  monopole  populations  do  not  sum  to  zero  for  the  estrogen  core.  To  maintain  electroneutrality,  it  is  essential 
that  appropriate  starting  values  be  assigned  to  other  substituents,  e.g.  for  estradiol,  values  of -0.47  and  0.29  would 
be  appropriate  monopole  populations  for  the  O  and  H  atoms  of  the  hydroxy  groups. 

3.  Conclusion 

Examples  of  typical  input  for  the  XD  program  and  pictorial  representation  of  the  coordinate  schemes  for 
the  hydrocarbon  core  of  estrogens  as  recommended  in  this  paper  are  available  as  supplementary  material. 
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Table  1  Atoms  used  to  define  vectors  Vi  and  V2  for  carbon  atoms® 


Atom(l) 

Atom  (2) 

Atom(3) 

Cl 

C2 

CIO 

C2 

C3 

Cl 

C3 

C4 

C2 

C4 

C5 

C3 

C5 

CIO 

C4 

C6 

Duml  (C5  -  Cl) 

Dum2  (C5  -  H6a) 

C7 

Dum3  (C6  -  C8) 

Dum4  (C6  -  H7a) 

C8 

Dum5  (C7  -  C9) 

Dum6  (C7  -  C14) 

C9 

Dum7(C8-Cil) 

Dum8  (C8  -  CIO) 

CIO 

Cl 

C5 

Cll 

Dum9  (C9  -  C12) 

DumlO(C12-Hlla) 

C12 

Dumll  (Cll  -  Cl 3) 

Duml2(C13-H12a) 

C13 

Duml3(C12-C14) 

Duml4(C14-C17) 

C14 

Duml5(C8-C15) 

Duml6(C8-C13) 

C15 

Duml7(C14-C16) 

Duml8(C14-H15a) 

C16 

Duml9(C15-C17) 

I  Dum20(C15-H16a) 

C17 

Dum21  (C13-C16) 

Dum22  (C16  -  02) 

C18 

Dum23  (H18a-H18b) 

Duni24(H18b-H18c) 

®Bond  lengths  of  atoms  used  to  calculate  dummy  atom  positions  (in  parentheses)  are  adjusted  to  be  equidistant  from  the  atom 
of  interest  to  ensure  a  tme  bisection. 


Table  2.  Atoms  used  to  define  vectors  Vi  and  for  hydrogen  atoms 


Atom(l) 

Atom(2) 

Atom(3) 

HI 

Cl 

C9 

H2 

C2 

HI 

H4 

C4 

01 

H6a 

C6 

H6b 

H6b 

C6 

H6a 

H7a 

C7 

H7b 

H7b 

Cl 

H7a 

H8 

C8 

C13 

H9 

C9 

H14 

Hlla 

Cll 

Hllb 

Hllb 

Cll 

Hlla 

H12a 

C12 

HI  2b 
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Table  3:  Initial  set  of  population  parameters  for  the  most  in^ortant  multipoles. 


sp^ 

sp^ 

Atom 

00 

20 

33+ 

32- 

Atom 

00 

10 

Cl 

-0.22(3) 

-0.22(5) 

0.32(5) 

HI 

0.24(4) 

0.14(4) 

-0.28(7) 

-0.21(5) 

0.34(5) 

H2 

0.24(4) 

0.15(4) 

C3 

0.15(4) 

-0.18(5) 

0.34(7) 

H4 

0.16(2) 

C4 

-0.26(5) 

-0.20(4) 

0.33(4) 

H6x 

0.15(2) 

C5 

-0.13(3) 

0.14(2) 

C6 

-0.27(4) 

0.31(7) 

H8 

0.13(5) 

Cl 

-0.28(6) 

0.31(6) 

H9 

C8 

-0.18(6) 

0.38(5) 

Hllx 

0.18(4) 

C9 

0.18(4) 

H15x 

0.12(4) 

-0.28(5) 

0.32(6) 

H16x 

0.17(5) 

0.15(3) 

-0.21(6) 

0.38(7) 

H17 

0.10(2) 

0.18(2) 

C14 

-0.15(4) 

0.34(3) 

H18x 

0.14(3) 

0.13(4) 

C15 

C16 

-0.35(5) 

C17 

0.16(4) 

0.30(5) 

C18 

-0.39(6) 

0.22(6) 

Hydrogen  labels  such  as  H6x,  represent  both  H6A  and  H6B.  Populations  are  average  values  from  6  studies  on  -  17i3- 
estradiol.!4MeOH,  17i8-estradiol.urea,  17of-estradiol,  17o<-estradiol.!^20,  16Q;17i3  -estriol,  estrone,  /c’s  for  H  atoms  were 
fixed  at  1.4.  Other  values  for  k  would  be  expected  to  introduce  small  changes  to  the  refined  monopole  populations,  but  to 
have  minimal  effect  on  the  higher  poles. 
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Appendix  D. 


Crystal  structure  and  charge  density  analysis  of  17a-estradiol*l/2  H2O 


1 


C{16) 


C(6)  C(7) 

Thermal  ellipsoid  plot  of  ITa-estradiohViHaO  where  ellipsoids  represent  50% 
probability  electron  density  of  the  atom.  Hydrogen  atoms  are  omitted  for  clarity. 


Data  collection  parameters  for  17a-estradiol*V2H20. 


Crystal  Data 

Chemical  Formula 

C18H25O2.5 

Temperature 

100.0(1) K 

Crystal  Dimensions 

0.24  X  0.33  X  0.33  mm 

Space  Group 

C2 

A 

19.0235(5)  A 

B 

7.0653(2)  A 

C 

13.3496(3)  A 

P 

124.0544(10) 

Volume 

1486.56(10)  A^ 

Z  (Crystallographic) 

4 

Integration  Parameters 

Box  Size  (°) 

Profile 
Fitting  (1/<t) 

Simple  Sum 
Perimeter  Limit 

Low  Angle 

1.2  X  1.2  X  0.8 

4010 

0.02 

High  Angle 

1.0  X  1.0  X  0.6 

3010 

0.02 

Reflection  Statistics  (from  SORTAV) 

Total  Reflections 

85540 

Rejected  Outliers 

69 

Unique  Reflections 

14593 

Average  Redundancy 

5.9 

Resolution 

1.319  A’* 

Completeness 

98.2  % 

Ri 

3.76  % 

R2 

4.13  % 

Rw 

13.52  % 

Z  (Refinement) 

1.949 

Table  2.  Selected  crystal,  integration,  and  reflection  data  for  17a-estradiol»*/4H20. 
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2-fold  axis 


Figure  2. 


Coordinate  system  for  the  water  molecule. 
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Rx 
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i?w 

z 

V 

Q<-4 

0 

0 

0.0 

0.0000 

0.0000 

0.0000 

0.000 

0.000 

-4<0<-3 

0 

0 

0.0 

0.0000 

0.0000 

0.0000 

0.000 

0.000 

-3<0<-2 

2 

1 

2.0 

0.3784 

0.3657 

0.3695 

1.086 

0.381 

-2<Q<-1 

43 

14 

3.1 

0.886 

1 

A 

o 

A 

o 

3.6 

■ESZI 

mm 

0<Q<  1 

6088 

3.9 

0.9922 

TOEiEB 

2.031 

mWMM 

1<Q<2 

mmm 

1657 

4.0 

0.5715 

■ingyiaL-y 

2<Q<3 

1299 

4.0 

0.3530 

wmjgm 

3<Q<4 

4586 

4.4 

0.2565 

0.3076 

2.257 

4<Q<6 

7017 

0.1805 

BtfcfclEB 

0.1905 

2.353 

0.195 

6<Q<8 

5775 

1128 

5.1 

0.1290 

0.1406 

2.153 

0.141 

8<Q<10 

5187 

5.8 

BilltfeVB 

0.1095 

1.999 

10<Q<20 

17218 

6.8 

0.0725 

20<Q<30 

16066 

1596 

10.1 

0.0340 

0.0467 

0.0399 

30<Q<50 

10194 

745 

13.7 

0.0244 

0.0340 

0.0283 

BKIttB 

mxnwim 

50  <  Q  <  100 

144 

17 

0.0170 

BiMEB 

100  <Q 

0 

0 

■wmnif 

0.0000 

0.000 

Table  3.  Intensity-Significance  Intervals  where  n  is  the  number  of  reflections,  m  is  the 

number  of  unique  reflections,  <n>  is  the  average  measurement  multiplicity,  and 
Q=I/Max  (oint/cext)  respectively  for  17a-estradiol»‘/2H20. 


n 

m 

<n> 

Rx 

Ri 

Rw 

Z 

V 

D>  1.029 

■on 

0.0273 

BtHltfcB 

0.030 

1.029  >D  >0.817 

BdrEB 

mm 

0.0329 

BtItkfifeB 

1.817 

0.039 

0.817>D>0.713 

6005 

743 

8.1 

0.0362 

0.0381 

0.1209 

2.082 

0.041 

0.713  >D> 0.648 

4332 

750 

5.8 

0.0384 

0.0368 

0.1320 

2.145 

0.044 

0.648  >  D  >  0.602 

4119 

748 

5.5 

0.0480 

0.0463 

0.1379 

2.182 

0.053 

0.602  >  D  >  0.566 

3914 

749 

5.2 

0.0565 

0.0537 

0.1396 

2.137 

0.064 

0.566  >  D  >  0.538 

3668 

735 

5.0 

BtlEKrB 

B»I<Vi?ll 

0.538  >D>  0.514 

4.8 

BiItf:W!iB 

BilfcbltB 

El^ 

0.100 

0.514  >D>  0.495 

3483 

752 

mm 

■iHiM 

0.1058 

0.1609 

2.089 

0.116 

0.495  >  D  >  0.478 

3201 

721 

mm 

BiinTdiB 

0.1015 

0.1629 

2.055 

0.117 

0.478  >D>  0.463 

3092 

720 

4.3 

0.1054 

0.0935 

0.1688 

2.063 

0.114 

0.463  >  D  >  0.449 

3023 

■>lEblB 

0.1327 

0.1945 

2.020 

0.163 

0.449  >D>  0.438 

2820 

msm 

0.2191 

1.970 

0.212 

0.438  >D>  0.427 

BiMEM 

0.2331 

1.947 

0.236 

686 

3.8 

BiliHlrJI 

0.2218 

BiMiitfiB 

2.007 

■owai^ 

BtWrgbB 

0.2656 

1.927 

MikVlM 

0.408  >  D  >  0.400 

2537 

0.3111 

0.2877 

BiWiVM 

1.932 

0.351 

0.400  >D>  0.393 

2331 

3.4 

0.3549 

0.3364 

BtltEMB 

mm 

BigfgB 

0.393  >D> 0.386 

2277 

3.3 

0.3852 

0.3783 

BilWda 

fern 

■»gfcW 

0.386  >  D  >  0.379 

1629 

513 

3.2 

0.4653 

0.4655 

0.4370 

1.769 

0.540 

Table  4.  Equal-Volume  Resolution  Shells  where  n  is  the  number  of  reflections,  m  is  the 
number  of  imique  reflections,  <n>  is  the  average  measurement  multiplicity,  and 
D=sin6>//1  (A"')  respectively  for  17a-estradiol»‘/2H20. 
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Monooole 

_ sp _ 

sp^ 

20 

33+ 

32- 

01 

-0.50 

02 

-0.49 

Cl 

-0.30 

-0.22 

0.34 

C2 

-0.38 

-0.19 

0.37 

C3 

0.27 

-0.21 

0.38 

C4 

-0.33 

-0.17 

0.36 

C5 

-0.18 

-0.22 

0.33 

C6 

-0.26 

0.31 

Cl 

-0.31 

0.34 

C8 

-0.21 

0.39 

C9 

-0.17 

0.31 

CIO 

-0.25 

-0.18 

0.37 

Cll 

-0.31 

0.35 

C12 

-0.28 

0.31 

C13 

-0.16 

0.38 

C14 

-0.20 

1 

0.38 

C15 

-0.26 

0.33 

C16 

-0.35 

0.42 

C17 

0.20 

0.38 

C18 

-0.32 

0.27 

Monopole 

HIO 

0.40 

H20 

0.38 

HI 

0.23 

H2 

0.22 

H4 

0.26 

H6x 

0.20 

H7x 

0.17 

H8 

0.20 

H9 

0.16 

Hllx 

0.17 

H12x 

0.16 

H14 

0.19 

H15x 

0.16 

H16x 

0.18 

H17 

0.13 

H18x 

0.18 

Atoms 

Kappa 

K 

k’ 

01,  02 

1 

0.97 

1.16 

C3 

2 

1.01 

0.92 

C17 

3 

1.02 

0.95 

Cl,  C2,  C4 

4 

0.97 

0.92 

C5,  CIO 

5 

0.98 

0.87 

C6,  Cl,  C8,  C9, 
C11,C12,C13,C14, 
C15,  C16,C17,C18 

6 

0.98 

0.95 

all  C-H  hydrogen  atoms 

7 

1.20 

1.29 

HIO,  H20 

8 

1.20 

1.29 

03 

9 

H30 

10 

Table  5.  Starting  values  entered  into  the  model  for  the  multipole  refinement  for  17a- 
estradiol»V2H20.  Units  for  multipole  populations  are  e'. 
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N 

-0.4373(7) 

0.5221(7) 

-0.0571(6) 

-0.2723(5) 

-0.2743(5) 

-0.0691(6) 

-0.0879(5) 

0.1345(5) 

0.1198(5)1 

0.0777(5) 

0.1432(5) 

0.1144(4) 

0.0792(5) 

0.3184(5) 

0.2897(5) 

0.3080(5) 

0.3488(5) 

0.2710(6) 

0.5271(6) 

0.4497(6) 

0.4791(4) 

0.2813(6) 

0.2676(6) 

0.1507(5) 

-0.4989(8) 

0.4125(12) 

0.2772(12) 

0.1964(10) 

0.2251(10) 

0.7893(9) 

0.9357(10) 

0.7980(8) 

0.8300(8) 

0.8759(9) 

0.5308(8) 

0.5067(8) 

0.1651(8) 

0.1816(9) 

0.2717(8) 

0.0821(8) 

0.5946(8) 

0.8213(9) 

0.6833(10) 
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0.5301(10) 
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Table  6.  Fractional  atomic  coordinates  for  17a-estradiol»54H20. 


Atom 

- - 

U22 

- iP - 

Uli 

iP 

U23 

01 

0.01787(9) 

0.01159(8) 

0.02309(10) 

0.02411(11) 

0.01353(9) 

-0.00461(9) 

0.00994(8) 

■tfiMMWAl 

■Bl 

0.01204(9) 

0.01148(9) 

0.01149(9) 

-0.00201(7) 

0.00414(8) 

KggnggBii 

0.01414(10) 

0.01165(9) 

-0.00233(8) 

0.00150(8) 

0.01535(10) 

0.01028(9) 

-0.00064(7) 

mmasBM 

0.00195(8) 

0.01196(9) 

0.01291(9) 

0.01167(9) 

0.00007(7) 

0.00282(7) 

■Eg 

0.01117(8) 

0.00945(8) 

0.01079(8) 

0.01588(10) 

0.00952(8) 

0.01253(9) 

0.00521(8) 

C7 

0.00788(8) 

0.01220(9) 

-0.00003(7) 

0.00565(8) 

0.00038(7) 

0.00806(7) 

0.01068(8) 

0.00082(6) 

0.00851(8) 

0.01017(8) 

0.00084(6) 

0.01026(8) 

0.00922(8) 

0.01016(8) 

mmsmiGm 

0.00499(7) 

Cll 

0.01317(9) 

0.00837(8) 

0.01196(9) 

0.00511(8) 

mmniMm 

C12 

0.01257(9) 

0.01020(8) 

0.01198(9) 

-0.00066(7) 

0.00537(8) 

0.00224(7) 

C13 

0.01036(8) 

0.01021(8) 

0.01170(9) 

0.00540(7) 

wmmM 

C14 

0.01012(8) 

0.00939(8) 

0.01077(8) 

0.00488(7) 

wmmem 

mm 

0.01331(10) 

0.01594(10) 

-0.00301(8) 

0.00472(8) 

0.00240(8) 

mm 

mmEM 

0.01801(11) 

0.01616(11) 

-0.00385(9) 

0.00231(9) 

0.00290(9) 

C17 

0.01204(9) 

0.01325(9) 

0.01289(9) 

0.00060(7) 

0.00418(8) 

0.00286(8) 

C18 

0.01651(11) 

0.01638(11) 

0.02061(12) 

msmm 

0.01234(10) 

03 

0.01581(12) 

0.01773(13) 

0.02136(14) 

MilililililiM 

0.00928(11) 

Table  7.  Anisotropic  thermal  parameters  of  non-H  atoms  for  17a-estradiol«y2H20. 


Atom 

Uiso 

H14 

0.0407(12) 

H15A 

0.0470(13) 

H15B 

0.0550(15) 

H16A 

0.0586(16) 

H16B 

0.0598(16) 

H17 

0.0426(12) 

H18A 

0.0577(15) 

H18B 

0.0644(17) 

H18C 

0.0538(14) 

H30 

0.1065(51) 

Atom 

Uiso 

HIO 

0.0386(18) 

H20 

0.0471(20) 

HI 

0.0519(16) 

H2 

0.0460(14) 

H4 

0.0447(14) 

H6A 

0.0564(15) 

H6B 

0.0449(13) 

H7A 

0.0413(13) 

H7B 

0.0434(12) 

H8 

0.0380(12) 

H9 

0.0392(12) 

HllA 

0.0436(13) 

HUB 

0.0464(13) 

H12A 

0.0410(12) 

H12B 

0.0445(13) 

Table  8,  Isotropic  thermal  parameters  of  H  atoms  for  17a-estradiol*}/2H20. 
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Atoms 

Bond  Length  (A) 

01 -C3 

1.3724(3) 

02-C17 

1.4365(4) 

C1-C2 

1.3927(4) 

Cl -CIO 

1.4024(4) 

C2-C3 

1.3945(4) 

C3-C4 

1.3951(4) 

C4-C5 

1.4030(3) 

C5-C6 

1.5162(4) 

C5-C10 

1.4087(3) 

C6-C7 

1.5289(4) 

C7-C8 

1.5274(3) 

C8-C9 

1.5457(3) 

C8-C14 

1.5261(3) 

Atoms 

Bond  Length  (A) 

C9-C10 

1.5232(3) 

C9-C11 

1.5390(3) 

C11-C12 

1.5390(4) 

C12-C13 

1.5305(4) 

C13-C14 

1.5432(3) 

C13-C17 

1.5444(4) 

C13-C18 

1.5431(4) 

C14-C15 

1.5382(4) 

C15-C16 

1.5560(4) 

C16-C17 

1.5465(4) 

Table  9.  Bond  distances  of  non-H  atoms  of  1 7a-estradiol»‘/2H20. 
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Table  10.  Bond  angles  of  17a-estradiol«’/2H20. 


Atom 

Monopole 
Population  {Pq.o) 

01 

6.519(12) 

02 

6.526(12) 

Cl 

4.222(23) 

C2 

4.254(22) 

C3 

3.855(20) 

C4 

4.247(22) 

C5 

4.127(21) 

C6 

4.217(22) 

C7 

4.217(21) 

C8 

4.127(21) 

C9 

4.122(21) 

CIO 

4.101(21) 

Cll 

4.226(20) 

C12 

4.236(20) 

C13 

4.189(22) 

C14 

4.121(22) 

CIS 

4.308(21) 

C16 

4.302(21) 

C17 

3.849(19) 

C18 

4.379(22) 

03 

3.274(8) 

Atom 

Monopole 
Population  {Po.o) 

HIO 

0.621(11) 

H20 

0.611(11) 

HI 

0.780(11) 

H2 

0.788(10) 

H4 

0.783(11) 

H6A 

0.853(9) 

H6B 

0.853(9) 

H7A 

0.854(7) 

H7B 

0.854(7) 

H8 

0.818(10) 

H9 

0.821(10) 

HllA 

0.858(8) 

HUB 

0.858(8) 

H12A 

0.852(8) 

H12B 

0.852(8) 

H14 

0.844(11) 

H15A 

0.853(8) 

H15B 

0.853(8) 

H16A 

0.851(9) 

H16B 

0.851(9) 

H17 

0.908(10) 

H18A 

0.879(7) 

H18B 

0.879(7) 

H18C 

0.879(7) 

H30 

0.726(8) 

Table  1 1 .  Monopole  populations  (e')  of  1 7a-estradiol»'/2H20. 
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Multipoles 

01 

02 

03 

Pi.^i 

-0.011(7) 

-0.036(7) 

0.0 

ll-i 

0.020(11) 

0.030(11) 

0.048(10) 

Pi.o 

0.029(8) 

0.014(6) 

0.0 

P2.0 

0.096(7) 

0.083(7) 

0.0 

P2.+1 

-0.037(6) 

-0.016(6) 

0.0 

P2A 

-0.035(7) 

-0.046(7) 

0.0 

P2.+2 

-0.064(7) 

-0.021(7) 

0.0 

P2.-2 

0.022(7) 

0.065(7) 

0.0 

. . P3.0 

0.014(12) 

-0.040(9) 

0.0 

Ps.+l 

-0.015(9) 

-0.081(8) 

0.0 

P3..1 

0.029(12) 

-0.017(13) 

0.0 

P3,+2 

0.026(10) 

-0.016(9) 

0.0 

Ps,-2 

0.022(11) 

0.089(13) 

0.0 

P3.+3 

0.124(8) 

0.077(9) 

-0.100(7) 

P3.-3 

0.020(11) 

-0.042(14) 

0.0 

P4.0 

0.0 

-0.080(11) 

0.0 

P4.-¥I 

0.0 

0.034(10) 

0.0 

P4rl 

0.0 

0.0 

0.0 

P4.42 

0.0 

0.011(10) 

0.0 

P4..2 

-0.018(10) 

0.021(11) 

0.0 

P4.+3 

-0.062(10) 

0.032(10) 

0.0 

P4..3 

0.0 

0.047(11) 

0.0 

P4.+4 

0.018(9) 

-0.023(10) 

-0.059(6) 

P4.-4 

-0.090(10) 

-0.071(11) 

0.0 

Table  12.  Multipole  populations  (e‘)  of  Oxygen  and  Nitrogen  atoms  of  17a-estradiol«V2H20. 
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Table  13.  Multipole  populations  (e  )  of  Carbon  atoms  of  17a-estradiol»‘/2H20. 


Atoms 

Pi.o 

P2.0 

HIO 

0.155(14) 

0.021(18) 

H20 

0.290(15) 

0.033(19) 

HI 

0.128(15) 

0.023(19) 

H2 

0.186(14) 

0.038(18) 

H4 

0.161(14) 

0.041(17) 

H6A 

0.192(9) 

0.026(11) 

H6B 

0.192(9) 

0.026(11) 

H7A 

0.140(8) 

0.041(10) 

H7B 

0.140(8) 

0.041(10) 

H8 

0.139(11) 

0.0 

H9 

0.107(11) 

0.015(15) 

HllA 

0.135(8) 

0.057(12) 

HUB 

0.135(8) 

0.057(12) 

H12A 

0.159(8) 

0.043(10) 

H12B 

0.159(8) 

0.043(10) 

H14 

0.107(13) 

0.038(17) 

H15A 

0.165(9) 

0.0 

H15B 

0.165(9) 

0.0 

H16A 

0.183(9) 

0.031(11) 

H16B 

0.183(9) 

0.031(11) 

H17 

0.193(13) 

0.0 

H18A 

0.134(7) 

-0.018(9) 

H18B 

0.134(7) 

-0.018(9) 

H18C 

0.134(7) 

-0.018(9) 

H30 

0.141(16) 

0.035(23) 

Table  1 4.  Multipole  populations  (e*)  of  Hydrogen  atoms  of  1 7a-estradiol»*/4H20. 
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Bond _ p(rc) _ V  p(rc) _ _ ^ _ 

C11-C12  1.550  -8.036  1.5406  0.7751  0.7655 


Table  17.  Topological  properties  of  bond  critical  points  in  the  hydrogen  bonds  of  17a-estradiol»14H20. 
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Figure  7.  The  Laplacian  of  the  total  electron  density  of  atoms  at  rest  in  the 

H12A-C12-H12B  plane  of  ITa-estradioh'AHaO.  Contour  intervals  are  5  eA'^ 
starting  at  5  eA'^  (solid  blue  lines),  -2  eA'^  starting  at  -2  eA"^  (dotted  red  lines), 
and  the  dashed  line  plots  0  eA'^. 
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Figure  1 1 .  Electrostatic  potential  isosurfaces  of  the  01  hydroxy  group  of  aE2-l/2water,  The 
blue  surface  in  each  picture  represents  +1.0  eA’*.  The  red  surfaces  represent  the 
following  potentials:  A  -0.05  eA'\  B  -0.10  eA‘\  C  -0.15  eA'*,  D  -0.20  eA"'. 


Electrostatic  potential  isosurfaces  of  the  02  hydroxy  group  of  aE2*V2water.  The  blue  surface  in 
each  picture  represents  +1.0  eA'\  The  red  surfaces  represent  the  following  potentials:  A  -0.05 
eA-\  B  -0.10  eA"^  C  -0.15  eA  \  D  -0.20  eA  \  E  -0.25  eA  ^  F  -0.30  G  -0.35  ek\  H  -0.40 
eA-\l-0.45  eA-\ 


Figure  12. 


Appendix  E. 


Crystal  structure  and  charge  density  analysis  of  17a-estradiol 
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Monopole  populations  from  the  multipole  refinement 
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Bond  critical  points  in  a-  estradiol 


Bond 

P 

e 

_ 

1^ 

£ 

C3-01 

2.10 

-21.9 

0.12 

— 

mm 

-23.3 

0.10 

C17-02 

1.8S 

-12.6 

0.11 

1.97 

-11.7 

0.06 

A-ring 

Ci-C2 

2.23 

-22.4 

0.32 

2.30 

-23.4 

0.19 

C2-C3 

2.32 

-25.6 

0.29 

128 

-23.4 

0.31 

C3-C4 

2.24 

-24.2 

0.35 

134 

-26.2 

0.22 

C5-C10 

2.20 

-22.6 

0.34 

123 

-23.7 

0.32 

Cl-ClO 

2.21 

-23.3 

0.33 

124 

-23.5 

0.20 

B-ring 

C5-C6 

1.82 

-14.2 

0.07 

1.75 

-118 

0.22 

C6-C7 

1.69 

-11.2 

0.12 

1.73 

-13.2 

0.03 

C7-C8 

1.65 

-11.6 

0.10 

1.73 

-13.6 

0.04 

C8-C9 

1.62 

-12.4 

0.06 

1.63 

-11.1 

0.04 

C9-C10 

1.76 

-14.9 

0.14 

1.69 

-115 

0.16 

C&D-ring 

C9-C11 

1.64 

-11.0 

0.13 

1.67 

-11.8 

0.12 

C11-C12 

1.67 

-11.5 

0.03 

1.67 

-10.8 

0.08 

C12-C13 

1.70 

-12.1 

0.09 

1.80 

-112 

0.03 

C13-C14 

1.63 

-11.2 

0.03 

1.70 

-11.4 

0.06 

C14-C15 

1.66 

-10.7 

EliEI 

1.65 

-10.5 

0.07 

C15-C16 

1.52 

-9.6 

Em 

1.64 

-10.4 

0.09 

C16-C17 

1.57 

-9.1 

0.06 

1.66 

-11.1 

0.01 

C17-C13 

1.69 

-11.9 

0.15 

1.60 

-10.4 

0.02 

C13-C18 

1.75 

-10.4 

0.09 

1.81 

13.1 

0.06 

p 

8 

P 

1^ 

£ 

01-HlA 

2.45 

-32.8 

0.01 

. 

141 

-34.3 

0.03 

02-H2A 

2.65 

-37.5 

0.04 

176 

-47.6 

0.04 

Cl-Hi 

1.86 

0.10 

1.84 

-18.8 

0.08 

C2-H2 

1.87 

-18.0 

0.09 

1.88 

-18.6 

0.09 

C4-H4 

1.88 

-19.6 

0.12 

1.90 

-19.6 

0.11 

C6-H6A 

1.80 

-16.2 

0.06 

1.84 

-17.6 

0.08 

C6-H6B 

1.74 

Em 

0.07 

1.77 

-16.2 

0.06 

C7-H7A 

1.82 

-16.9 

0.02 

1.89 

-18.0 

0.06 

C7-H7B 

1.74 

-14.9 

0.11 

1.76 

-16.3 

0.07 

C8-H8 

1.78 

-13.2 

0.04 

1.77 

-17.0 

0.01 

C9-H9 

1.81 

-15.6 

0.01 

1.83 

-16.0 

0.02 

Cll-HllA 

1.73 

-14.3 

0.05 

1.80 

-15.0 

0.02 

Cll-HllB 

1.78 

-14.7 

0.08 

1.89 

-17.7 

0.06 

C12-H12A 

1.74 

-14.9 

0.04 

1.87 

-18.0 

0.02 

C12-H12B 

1.81 

-12.2 

0.06 

1.82 

-17.3 

0.03 

C14-H14 

1,87 

-18.2 

0.06 

1.83 

-16.4 

0.03 

C15-H15A 

1.70 

-14.3 

0.05 

1.82 

-15.9 

0.10 

C15-H15B 

1.76 

-15.2 

0.08 

1.74  : 

-14.7 

0.08 

C16-H16A 

1.88 

-18.2 

0.06 

1.79  ' 

-14.8 

0.19 

J 

1.85 

-17.6 

0.03 

1.83  1 

-16.3 

0.14 

C17-H17 

1.91  i 

-18.5 

0.04 

1.80  I 

-16.7 

0.03 

Cl8-Hi8A 

1.82  i 

-113 

0.14  : 

1.95  : 

-19.3 

0.03 

C18-H18B 

1.63  1 

mm 

0.15 

1.85 

-16.5 

0.15 

C18-H18C 

1.81 

-13.7 

0.14 

1.93 

-18.5 

0.10 

9 


MODEL  2 
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Monopole  populations  from  the  multipole  refinement 
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Figure  2.  Model  multipole  deformation  electron  density  in  the  aromatic  ring  plane. 

Contour  interval  is  0.05  eA‘^. 
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Bond  critical  points  in  a-  estradiol 
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Crystal  structure  and  charge  density  analysis  of  3,16a,17|3-estriol 
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Charge  density  and  electrostatic  potential  study  on  estrogen  molecules: 

16a,  17p  -  Estriol 

Abstract 

The  measurements  of  the  electron  density  distribution  of  the  principal  members  of  the  estrogen 
family,  the  derivatives  estriol,  estradiol  and  estrone  are  now  in  progress.  As  a  part  of  our  study,  the 
electron  density  distribution  of  16a,17P-estriol  (C18H24O3)  was  determined  from  the  high-resolution 
single  crystal  X-ray  diffraction  measurements  at  100  K.  Estriol  crystallizes  with  two  molecules  in  the 
asymmetric  unit.  Intensity  data  were  collected  using  SMART  2K  CCD  area  detector  with  Ag  Ka 
radiation  up  to  resolution  (sin0/A,)max  =  1*29  A'*.  The  structure  was  solved,  and  a  conventional  spherical 
atom  refinement  was  carried  out.  A  multipole  pseudo  atom  refinement  was  then  performed  using  the 
program  XD.  The  electron  density  at  the  critical  point  of  aromatic  C-C  bonds  are  larger  than  the  other 
ring  systems  in  the  molecule.  The  electrostatic  potential  has  been  calculated  for  both  molecules  removed 
from  the  crystal  lattice.  Strong  electro  negative  regions  were  found  in  the  vicinity  of  the  oxygen  atoms 
of  the  16a,  17p  regions  of  both  molecules.  A  small  net  negative  electrostatic  potential  aroimd  3- 
hydroxyl  regions  in  both  molecules  may  be  related  to  low  binding  affinity  of  estriol.  Net  atomic  charges 
were  calculated  from  the  kappa  refinement.  A  group  charge  calculation  indicates  that  the  central  region 
is  close  to  neutral,  but  the  head  and  tail  regions  form  charged  units.  The  topological  analysis  of  the 
hydrogen  bonds  indicates  the  strength  of  the  intermolecular  interactions. 

Introduction 

In  the  recent  years  much  attention  has  been  paid  on  estrogen  research'.  The  major  role  of 
estrogen  molecules  in  the  body  influences  the  growth,  development,  and  behavior,  and  affect  many  other 
body  parts  such  as  tissues,  skin,  bone,  liver.  Importantly,  this  family  of  steroid  hormones  regulates  and 
sustains  female  sexual  development  and  reproductive  fimction^®.  In  other  words,  it  stimulates  tissue 
growth  by  promoting  cell  proliferation  (DNA  synthesis  and  cell  division)  in  the  female  sex  organs 
(breast,  uterus)^'’.  Subtle  changes  in  the  structure  of  the  estrogen  molecules  affect  their 
chemical/biological  behavior  and  can  result  in  the  development  of  growth  inhibitors  for  tumors.  For 
example,  moving  the  hydroxy  group  from  C(3)  (estradiol)  to  one  of  the  adjacent  C-atoms  can  change 
the  carcinostatic  potentials  from  agonistic  (estradiol)  to  inhibitory  (  2-and  4-hydroxy  estratrien  -17P- 
ol)"'’^.  In  this  way,  more  than  60%  of  all  breast  cancers  are  known  to  be  hormone  dependent,  i.e.. 
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initiation  and  progress  can  be  influenced  by  estrogens  and  related  compounds.  These  molecules  bind  as 
ligands  to  the  estrogen  receptor  and  hereby  initiate  a  series  of  events  resulting  in  the  activation  or 
repression  of  selective  genes.  While  most  changes  modify  the  affinity  of  the  ligand  for  the  estrogen 
receptor,  this  does  not  necessarily  correlate  with  the  stimulation  ability  of  the  transcription  of  estrogen 
responsive  genes.  Several  articles  focused  on  identification  and  characterization  of  receptor  ligand 
binding  affinities.  A  recent  comparative  QSAR  analysis  of  estrogen  receptor  ligands  generalizes  the 
effect  of  different  substituents  in  estrogen  molecules;  more  importantly,  the  substituents  which  increase 
the  electron  density  in  the  A-ring  appear  to  increase  the  binding  affinity,  and  the  steric  character  of 
substituent  and  polar  character  reduce  the  ligand  receptor  binding  affinity^.  An  extensive 
crystallographic  structural  study  reports  the  conformational  features  of  these  molecules  and  proposing 
the  possibilities  of  ligand  binding^’^.  The  knowledge  of  receptor  or  receptor  complex  model  is  essentially 
important  to  understand  the  binding  mechanism.  In  this  context  the  structural  analysis  of  ligand  binding 
domain  (LBD)  of  two  estrogen  receptors  (a  and  P)  complexed  with  various  ligands  has  been 
reported.  These  studies  emphasize  the  agonist  and  antagonist  effect  of  certain  hormones  and  drugs  with 
the  receptor.  Furthermore,  molecular  modeling  of  a  series  of  estrogen  molecules  has  been  conducted 
suggesting  that  conformational  flexibility  is  the  important  property  of  speeific  ligands  for  the  estrogen 
receptor*^.  However,  still  the  mechanism  by  which  the  ligands  regulate  the  gene  expression  is  currently 
unknown,  but  it  has  been  suggested'^  that  the  differences  in  eleetrostatic  potential  (ESP)  can  contribute 
to  the  variation  in  the  regulation  of  hormone  dependent  genes.  This  information  ean  be  obtained  from 
studies  of  the  electronic  properties  in  estrogen  derivatives. 

An  experimental  multipole  representation  of  the  eharge  distribution  has  been  becoming  very 
successful  method  to  determine  the  eleetrostatic  properties  and  to  investigate  the  intermolecular 
interactions.  An  interesting  study  on  the  determination  of  ESP  from  multipole  moments  shows  the 
importance  of  multipole  moments  in  the  intermolecular  interactions*'*’^^.  Currently,  various  quantum 
mechanical  models  are  available  to  explore  the  charge  distribution,  moleeular  ESP  and  its  topological 
characteristics  using  multipole  model*^'*®.  However,  it  is  hi^ly  desirable  for  the  model  to  include  the 
effect  of  intermolecular  interactions  to  understand  the  molecular  recognition  processes.  The  multipole 
model  experimental  electron  density  analysis  of  molecules  by  the  X-ray  diffraetion  methods  ineludes 
this  aspect.  In  the  recent  years,  experimental  multipole  model  results  have  been  extensively  analyzed  to 
elucidate  the  electronic  properties  of  atoms  in  molecules.  A  widely  aecepted  Bader’s  theory^**  of  “Atoms 
in  Molecules”  is  a  useful  tool  to  obtain  electronic  properties  such  as  electron  density  at  the  critical  point, 
Laplacian  of  charge  density  (charge  concentration  or  depletion),  bond  ellipticity,  eharges  of  atoms  and 
the  electrostatic  potential,  which,  in  turn,  aid  in  the  prediction  of  how  this  molecule  might  react  with  a 
given  environment.  Moreover,  the  hydrogen  bond  topology  studies  characterize  the  strength  of 
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intermolecular  interactions. 

We  have  begun  a  systematic  study  of  estrogen  compounds  by  investigating  crystals  of  the 
derivatives  of  the  principal  estrogens,  estriol,  estradiol  and  estrone.  The  first  study  in  this  series,  on 
estra-l,3,5(10)-triene-3,16a,17P-triol,  (see  Figure  1)  is  reported  here.  Estriol  is  the  one  of  the  active 
estrogens  found  in  the  body.  The  main  source  of  estriol  is  the  placenta,  which  produces  a  large  amount 
of  estriol  during  pregnancy^',  while  the  other  two  estrogens  are  produced  in  comparatively  small 
amounts.  Estriol  does  not  associated  with  the  cancer  activity  in  the  female  body  but  it  protects  this 
disease.  Estriol  has  much  less  stimulating  effect  on  the  breast  and  utrine  lining  than  estradiol  and 
estrone.  Receptor  binding  studies  have  indicated  that  estriol  has  low  binding  affinity  to  estrogen 
receptors  and,  in  general,  it  has  the  shortest  receptor  occupancy.  From  the  distinct  biological 
characteristics  of  estriol  over  the  other  two  estrogens  (estrone  and  estradiol),  it  implies  that  these 
molecules  must  differ  in  some  way  in  their  electronic  properties  and/or  physical  properties.  The 
knowledge  of  the  properties  of  these  estrogens  may  help  to  understand  the  detailed  mechanism  of 
estrogenic  activity  of  these  molecules.  The  commercially  available  estriol  has  been  crystallized  from 
acetonitrile.  We  have  carried  out  the  crystallographic  structure  and  electron  density  analysis  from  the 
high  resolution  X-ray  diffraction  measurement. 

Experimental  Section 


1.  Data  collection 

Estriol  has  been  crystallized  from  acetonitrile  by  slow  evaporation  at  room  temperature.  The 
crystals  were  found  to  be  colorless  and  prismatic.  In  order  to  find  a  sample  with  good  reflection  profiles, 
several  crystals  were  examined  on  the  diffractometer  for  the  low  and  high  angle  scattering.  A  well 
shaped  0.48  x  0.26  x  0.24  mm  size  crystal  was  selected  and  mounted  on  a  Bruker  SMART  2K  CCD 
diffractometer  using  graphite  monochromatized  Ag  Ka  (X  =  0.56089A)  radiation  (50V,  35mA).  As  a 
preliminary  check,  a  room  temperature  data  set  was  collected  for  a  hemisphere.  Refinement  of  the 
structure  from  this  data  gave  results  comparable  to  the  reported  crystal  structure^^.  Then,  the  crystal  was 
cooled  to  100  K  using  an  Oxford  Cryostream  N2  open  flow  cryostat.  After  cooling,  no  significant  unit 
cell  modification  was  observed  and  the  unit  cell  was  very  stable  over  the  period  of  the  entire  data 
collection. 

Initially  the  orientation  matrix  and  the  unit  cell  parameters  were  obtained  from  least-squares 
refinement  using  88  well  centered  reflections  measured  in  three  sets  of  matrix  runs  of  each  20  frames. 
The  low  temperature  cell  obtained  in  this  way  was  ~2  %  smaller  than  that  at  room  temperature.  A  full 
sphere  of  low  and  high  angle  reflection  intensities  were  measured  using  two  detector  settings  at  20  =  0® 
and  20  =  -50°,  using  the  co  scan  method  for  different  9  angles  (  0,  45,  90,  135,  180,  225,  270  and  315°). 
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The  CCD  detector  was  placed  at  5.1  cm  from  the  crystal  to  minimize  the  reflection  overlap.  A  total  of 
600  frames  were  collected  for  each  cp  angle  with  a  scan  width  of  0.3°.  The  frame  times  for  low  (40  sec) 
and  high  angles  (160  sec)  were  picked  after  the  careful  observation  of  the  most  intense  reflection, 
avoiding  detector  overflow,  but  utilizing  the  full  dynamic  range.  The  crystal  decay  was  monitored  by  re¬ 
measuring  the  initial  50  frames  at  the  end  of  each  data  collection,  and  it  was  found  to  be  negligible.  The 
data  collection  strategy  was  checked  with  the  program  ASTRO^^  incorporated  into  the  software  package 
SMART^^.  We  obtained  100%  coverage  for  a  resolution  of  (sin0/A,)niax  =1.2  A'^  with  good  redimdancy. 
A  total  of  9600  frames  were  collected  over  the  period  of  two  weeks. 

2.  Data  reduction  and  averaging 

The  best  unit  cell  and  the  orientation  matrix  were  determined  by  thresholding  [I/a(I)  >  20]  of 
reflections  from  each  run  of  the  data  set.  This  orientation  matrix  fits  well  with  the  entire  reciprocal  space 
observed.  The  data  has  been  integrated  and  corrected  for  the  background,  Lorentz  and  polarization 
effects  using  the  software's  advanced  integrating  procedures  SAINT^^.  The  method  of  integration  of  the 
collected  frames  was  as  described  by  Kabch^®  for  area  detectors.  During  the  integration  the  orientation 
matrix  and  xmit  cell  were  not  allowed  to  vary.  Appropriate  three  dimensional  box  sizes  were  used  in  low 
and  high  angle  integrations  to  account  for  the  Kai  and  Ka2  splitting  in  hi^  order  reflections.  No 
absorption  correction  was  applied  because  of  the  low  value  of  the  absorption  coefficient  (0.055  mm’’). 
All  the  corrected  reflection  intensities  were  scaled  and  averaged  based  on  the  point  group  symmetry 
using  the  program  SORTAV^^.  A  total  of  277952  reflections  were  measured  up  to  sinO/X,  =  1.29  A'\ 
43918  are  unique  reflections  with  an  average  redundancy  of  6.2.  Of  these,  989  unique  data  were 
measured  only  once,  3028  twice  and  39882  three  or  more  times.  We  note  that,  in  the  complete  set  of 
unique  reflections,  29528  reflections  were  found  to  be  greater  than  3a(I).  The  internal  agreement 
factors^^  based  on  the  averaging  of  symmetry  equivalent  reflections  are  R1  =  0.044,  R2  =  0.031,  wR  = 
0.11  and  S  =  1.4.  The  unit  cell  parameters,  data  collection,  averaging  and  other  crystallographic 
parameters  are  summarized  in  Table  1. 

3.  Conventional  refinement 

The  room  temperature  crystal  structure^^  of  16a,  17p-estriol  has  been  previously  reported.  We 
solved  the  100  K  low  temperature  structure  by  direct  methods  and  the  spherical  atom  least  squares 
refinement  was  carried  out  using  the  program  SHELXTL^®.  All  the  hydrogen  atoms,  except  one,  were 
located  from  a  difference  Fourier  synthesis.  In  the  least  squares  refinement,  all  non-hydrogen  atoms 
were  treated  anisotropically,  and  the  hydrogen  atoms  were  isotropic.  The  refinement  converged  with  all 
data  to  R  =  0.058,  wR  =  0.151  GOF  =  1.01.  Comparison  of  Fo  and  Fc  values  of  die  largest 
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structure  factors  gave  no  indication  of  the  effect  of  extinction  in  the  data.  A  thermal  ellipsoid  diagram 
(Figure  2)  shows  the  two  molecules  in  the  asymmetric  unit  and  the  standard  atom  numbering  scheme. 

4.  Multipole  refinement 

In  order  to  model  the  deformation  electron  density  and  explore  the  bond  topological 
characteristics  like  p(r),  V^p(r),  bond  ellipticity  (s),  molecular  electrostatic  potential  etc.,  a  multipole 
aspherical  atom  refinement  was  carried  out  using  the  Hansen-Coppens  multipole  formalism 
implemented  in  the  program  XD^®.  According  to  this  formalism,  the  total  electron  density  of  an 
aspherical  atom  is  parameterized  using  an  atom-centered  multipole  expansion 

3  ^max  3  ^ 

Patom(r)  =  Pc  pc(r)  +  PyK^  pv(Kr)  +  Z k'i  Ri  ( K'r)! Pim±  dim±  ( 9,  (p) 

1=0  m=0 

The  first  two  terms  pc(r)  and  Pvfnf)  are  the  spherically  averaged  Hartree-Fock  atomic  core  and  valence 
densities  normalized  for  one  electron;  the  third  term  describes  the  aspherical  part  of  the  pseudo  atom; 
Ri(r)  =  [ai  "‘'^^/(ni+2)!]r^‘  are  the  normalized  radial  distribution  functions;  oj  are  the  Hartree  Fock 
optimized  single  zeta  values;  k  and  id  are  radial  expansion/contraction  parameters  which  are  optimized 
in  the  least-square  refinement  along  with  Py  and  Pim,  Pim  are  the  multipole  population  coefficients;  dim± 
are  density  normalized  real  spherical  harmonic  angular  functions  and  r,  0  and  (p  are  user  defined  local 
atom  centered  coordinates  expressed  in  a  polar  coordinate  system.  In  the  multipole  refinement,  all 
carbon  and  oxygen  atoms  were  allowed  to  expand  up  to  the  hexadecapole  and  the  hydrogen  atoms  up  to 
the  dipole  level.  The  k  and  id  parameters  of  all  non-H  atoms  were  refined,  but  hydrogen  atoms  id  was 
fixed  at  1.2. 

In  order  to  obtain  accurate  positional  and  thermal  parameters  for  the  non-hydrogen  atoms,  a 
high-angle  refinement  was  carried  out  using  the  data  with  sin0/A,  >  0.7  A’\  All  C-H  distances  were  then 
reset  to  the  neutron  bond  lengths  [Csp2-H  =  1.083,  Csps-H  =  1.099(CH),  1.092(CH2),  1.059(CH3)  and 
0-H  =  0.967  A]  as  reported  in  the  International  Tables^  ^  This  high-angle  refinement  converged  without 
any  significant  change  of  scale  factor,  thus  indicating  the  quality  of  the  high  angle  data.  The  multipole 
refinement  was  carried  out  on  the  25520  reflections  [I  >  3ct(I)]  with  the  resolution  (sin0/A,)max  =  1-1  A''. 
The  refinement  strategy  was  as  follows: 

In  the  first  step,  Pv,  k  and  scale  factor  were  refined;  second  step  P/^,  k';  third  step  xyz,  Ujj,  scale 
factor.  Finally,  xyz,  and  all  electronic  parameters  were  refined  together.  Since  there  are  two 
molecules  in  the  asymmetric  unit,  their  multipole  parameters  were  constrained  to  be  equal  at  the 
beginning  of  the  refinement,  except  the  0-H  group  atoms  as  these  are  involved  in  the  strong  hydrogen 
bonding.  In  addition,  to  avoid  the  phase  problem^^,  all  chemically  equivalent  atoms  were  initially 
constrained  to  be  equal,  but  no  site  symmetry  was  imposed  on  any  atoms  in  the  molecules.  At  the  later 
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stages  of  the  refinement,  the  chemical  atom  constraints  within  the  molecules  were  released,  but  not 
between  the  molecules.  The  residual  index  from  the  final  refinement  is  given  in  Table  1.  An  attempt  was 
made  to  refine  the  two  molecules  independently  after  the  final  refinement.  This  refinement  did  not  give 
a  meaningful  result  so  we  did  not  proceed  further.  The  net  atomic  charges  were  calculated  from  a  k- 
refinement^^  (Pv,  k,  and  scale).  During  the  refinements,  the  unit  cell  was  constrained  to  be  neutral  and  no 
charge  transfer  was  allowed  between  the  molecules.  The  featureless  residual  density  maps  of  the  two 
molecules  (Figure  4)  from  the  Fourier  summation  (Fo-Fmui)  indicates  the  adequacy  of  the  model 
refinement. 

It  is  essential  to  test  the  physical  significance  and  the  correctness  of  the  thermal  parameters 
obtained  from  the  model  refinement.  To  accomplish  this,  Hirshfeld's^"^  "rigid  bond"  test  was  performed 
for  both  the  molecules  leading  to  a  satisfactory  result.  The  maximum  difference  Aa,b  between  the  atomic 
mean  square  displacement  amplitudes  along  the  interatomic  vectors  was  0.0013  A. 

Result  and  Discussion 


1.  Structural  aspects 

The  report  of  the  room  temperature  crystal  structure^^  indicated  significant  conformational 
differences  between  the  two  molecules  in  the  asymmetric  unit.  We  observe  some  notable  differences 
between  the  reported  room  temperature  and  current  low  temperature  structure,  so  it  is  worth  revisiting 
the  structural  analysis.  Figure  2  shows  the  thermal  ellipsoid  representation  of  estriol  with  two 
independent  molecules  in  the  asymmetric  unit  obtained  from  the  low  temperature  (100  K)  measurement. 
The  geometrical  comparison  of  the  two  molecules  explicitly  indicates  how  the  two  molecules  are 
conformationally  different.  On  comparing  the  bond  distances  between  the  two  molecules,  no  significant 
difference  were  observed,  in  contrast  to  the  previous  report^^.  The  maximum  deviation  in  the  bond 
lengths  of  the  ring  (A-D)  systems  between  the  two  molecules  is  0.01  A  (Table  2).  There  is  no  evidence 
for  C(2)-C(3)  bond  shortening  in  the  A-ring  systems  as  the  room  temperature  structure  reports^^.  The 
average  Car-Car  bond  distances  of  A-ring  atoms  of  molecule  1  and  2  are  1.3995  and  1.3990  A 
respectively.  The  CH3  group  of  the  second  molecule  lies  above  the  aromatic  ring  of  the  first  producing 
distortion  in  the  aromatic  ring.  This  effect  is  reproduced  here  giving  a  folding  in  the  A-ring  of  molecule 
1  with  a  dihedral  angle  of  2.99°.  There  is  no  such  effect  in  the  second  molecule  where  the  A-ring 
system  is  planar  (0.52°).  There  are  differences  in  bond  angles  between  the  molecules  1  and  2  around  the 
atoms  C(9),  C(13)  and  C(14)  [Table  2].  However,  we  observe  a  maximum  deviation  of  2.5°  compared  to 
3.9°  as  previously  reported.  The  torsion  angles  in  the  Table  2  indicate  the  distortion  of  the  B  and  C  ring 
systems.  The  steric  hindrance  between  the  C-ring  equatorial  hydrogen  atom  at  C(ll)  and  the  hydrogen 
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atom  at  C(l)  produce  different  conformations  of  the  B  and  C  ring  systems  in  the  two  molecules. 
Comparison  of  torsion  angles  between  the  molecules  shows  the  major  distortion  appears  at  the  head  and 
tail  part  of  the  molecules  (see  Table  2).  This  is  attributed  to  the  different  environment  of  intermolecular 
interactions  of  the  molecules  (1  and  2)  with  their  neighbors.  The  relative  strengths  of  the  hydrogen 
bonds  to  0(3)  in  the  two  molecules  was  previously  discussed  on  the  basis  of  0.04  A  difference  in  the 
C(16)-0(3)  bond  length.  We  observe  less  than  0.01  A  difference  in  this  distance  between  molecules  1 
and  2.  Hydrogen  bonding  does,  however,  have  an  effect  on  the  geometry  of  the  hydroxyl  groups.  In  the 
second  molecule  the  16a  and  17P-hydroxyl  C-O-H  bond  angles  are  systematically  less  compared  with 
the  first,  the  maximum  difference  being  8.8°.  There  is  less  of  a  differentiation  in  the  3-hydroxyl  bonds, 
the  deviation  being  2.1°.  The  structure  is  characterized  by  strong  hydrogen  bonding.  Figure  3  depicts 
the  hydrogen-bonding  environment  of  the  two  molecules  in  the  asymmetric  unit.  Each  hydroxyl  group  is 
involved  in,  at  least,  two  hydrogen  bonds.  The  hydrogen  bonds  0-H— O  and  O—H-0  type  of 
interactions  in  the  molecules  indicates  the  0-atoms  behave  as  donor  as  well  as  acceptor.  Moreover,  0(T) 
has  one  additional  hydrogen  bond  resulting  in  a  tetrahedral  coordination  of  the  0-atom.  The  head  to  tail 
hydrogen  bonds  link  the  molecules  into  chains.  The  hydrogen  bond  geometrical  parameters  are  given  in 
Table  3. 

2.  Charge  density  distribution 

The  model  static  deformation  density  map  in  Figure  5  shows  the  bonding  features  of  the  different 
atoms  in  the  estriol  molecule.  All  chemical  bonds  show  an  increase  in  the  electron  density  with  respect 
to  neutral  spherical  atoms  as  well  as  density  associated  with  the  lone  pair  regions  of  the  oxygen  atoms. 
The  featureless  residual  electron  density  map  was  computed  using  the  program  XDFOUR^^  after  the 
multipole  refinement.  The  featureless  residual  maps  (Figures  4a  and  b)  of  the  two  molecules  in  the 
asymmetric  unit  confirm  the  validity  and  the  adequacy  of  the  model  refinement  and  the  quality  of  the 
data  set.  The  minimum  and  maximum  residual  density  calculated  for  the  molecules  are  0.21,  -0.29 
(moleculel)  and  0.22,  -0.27eA'^  (molecule  2)  respectively.  As  specified  earlier,  in  the  multipole 
refinement,  the  respective  chemically  equivalent  ring  atoms  of  the  two  molecules  were  constrained  pair 
wise.  No  constraint  was  imposed  on  the  hydroxyl  groups,  since  they  involved  in  strong  hydrogen 
bonding.  All  further  discussions  are  based  on  this  constrained  model. 

In  order  to  further  characterize  the  interaction  of  atoms  in  a  molecule,  a  topological  analysis  was 
performed  based  on  Bader’s  theory^*^  of  “Atoms  in  Molecules”  implemented  in  XDPROP^^.  This 
analysis  includes  the  determination  of  electron  density,  p(rc),  at  the  bond  critical  point  (CP),  where 
Vp(rc)  =  0,  between  the  bonded  atoms,  its  Laplacian,  V  p(rc),  and  bond  ellipticity  [e  =  (>.1/^.2)-!].  Also, 
the  eigenvalues  of  the  second  derivatives  of  p(rc)  was  calculated  to  characterize  the  nature  of  the 
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bond  critical  point.  The  bond  path  analysis^®’^^  gives  information  about  the  amount  of  deviation  of  the 
bond  critical  point  from  the  intemuclear  axis,  thus  obtaining  an  additional  measure  of  a  strain.  All  these 
results,  including  the  length  of  bond  paths  and  intemuclear  distances,  are  listed  in  Table  4. 

A-ring  system:  The  topological  analysis  results  reported  in  Table  4,  such  as  Xi,  X2,  ^3,  p(rc)»  V  p(rc),  8 
reveals  the  electronic  insight  into  the  estriol  molecules.  The  magnitude  of  p(rc)  at  the  bond  critical  point 
is  a  clear  indicator  of  the  chemical  bond  strength.  The  electron  density,  p(rc),  at  the  critical  points  of  the 
aromatic  ring  (Csp2-Csp2 )  bonds  range  from  2.04  to  2.19  eA’^.  The  average  value  (2.1  eA'^)  is  a  typical 
value  for  aromatic  bonds,  but  these  values  are  0.1  eA'^  smaller  than  in  the  reported  p-nitrophenol^*.  The 
Laplacian,  V^p(rc),  is  a  very  sensitive  parameter  which  identifies  the  areas  where  charge  is  concentrated 
(negative  value)  or  depleted  (positive  value)  in  the  bonding  region.  In  this  respect,  we  observe  the 
maximum  negative  Laplacian  in  the  A-ring  bonds  being  -22.1(1)  eA'^,  and  this  value  appears  less  in  the 
bonds  linked  to  B  ring  atoms.  However,  the  average  value  is  -18.52  eA'^,  which  is  very  close  to  an 
earlier  observation  for  the  benzene  molecule^^.  The  dispersion  of  bond  ellipticity  is  attributed  to  the 
preferential  charge  accumulation  in  the  bond.  Bader’s  theory  on  ellipticity^*^’^*^  emphasis  the  concept  of 
the  CT-7I  character  of  the  bond.  Explicitly,  it  is  the  measure  of  deviation  of  charge  density  from  spherical 
symmetry  (at  the  critical  point).  In  this  aspect  the  average  value  (0.18)  of  ellipticity  of  Csp2-Csp2  bonds 
indicates  the  significant  7t  (aromatic)  character.  This  value  is  smaller  than  reported  for  the  benzene  ring 
bonds  but  completely  agree  with  those  in  tris-annealed  benzene"*  ^  The  p(rc)’s  in  the  Csp2-H  bonds  range 
from  1.78(4)  to  1.90(1)  eA'^,  and  this  density  is  consistent  with  one  would  expect  for  the  single  C-H 
bonds.  The  Laplacian  of  the  C(l)-H(l)  bond  is  small  compared  with  the  other  two  C-H  bonds,  perhaps 
indicating  weakening  of  the  bond  due  to  steric  interactions  with  hydrogen  atoms  at  C(ll).  The  small 
ellipticity  of  Csp2-H  bonds  is  consistent  with  only  a  ct  contribution  in  the  single  bonds. 

B  and  C  ring  systems:  The  density  of  chemically  equivalent  Csp2-Csp3  bonds  in  the  B-ring  are  equal  and 
their  Laplacian’s  are  very  close  (see  Table  4).  In  both  ring  systems,  the  average  values  (1.70  and  1.67 
eA‘^)  of  p(rc)’s  in  the  Csp3-Csp3  bonds  are  not  very  different  and  a  similar  trend  appears  for  the 
Laplacian.  The  p(rc)’s  and  Laplacian  for  these  bonds  are  significantly  less  than  the  aromatic  ring  system 
in  agreement  with  their  single  bond  character.  In  the  C-ring  the  ellipticity  of  the  bonds  C(1 1)-C(12)  and 
C(12)-C(13)  are  unexpectedly  larger  compared  to  the  similar  bonds  having  small  values  (see  Table  4), 
indeed  the  values  approach  those  of  the  aromatic  ring.  In  all  Csp3-H  (CH2)  bonds  the  density  ranges 
from  1.64  to  1.86  eA'^  and  Csp3-H  (CH)  bonds  1.75  to  1.81  eA'^.  The  average  effect  is  same  in  both 
types  of  bonds.  The  Laplacian  of  these  bonds  is  consistently  smaller  than  for  the  Csp2-H  bonds.  The 
maximum  and  minimum  values  are  -8.9(2)  and  -18.6(1)  eA'^ . 

D-ring  system:  Having  5  members  in  the  ring  and  due  to  the  different  hydrogen  bonding 
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environments,  this  ring  is  highly  strained  in  both  molecules  (1  and  2).  Nevertheless,  no  significant 
difference  was  observed  in  the  p(re)’s  at  the  bond  critical  points,  except  for  C(14)-C(15)  =  1.57  eA'^ 
implying  a  weakening  of  this  bond.  The  average  value  of  the  Laplacian  is  -12.88  eA'^.  We  observe  the 
ellipticity  of  this  five-member  system  being  close  to  the  B-ring  bonds,  and  it  is  relatively  small, 
consistent  with  the  single  bonds.  The  electron  density  (1.94  eA'^ )  and  the  Laplacian  of  tiie  C(16)-H(16) 
bond  is  surprisingly  higher  than  the  similar  C(17)-H(17)  bond  in  this  ring  (see  Table  4). 

C-O  and  O-H  groups:  In  general,  diese  groups  play  a  significant  role  in  the  whole  geometry  of 
estrogen  molecules.  Especially  the  OH  group  perturbs  the  electron  density  of  aromatic  system  and  its 
electron  density  distribution  is  important  in  the  ligand  binding  domain^^.  The  Pc(i')’s  of  Csp2-0  [A-rings] 
bonds  of  both  molecules  (1  and  2)  are  larger  than  pc(r)’s  of  the  Csp3-0  bonds  of  D-ring,  in  agreement 
with  conjugation  of  0(1)  to  the  aromatic  ring.  The  average  values  are  2.11  and  1.89  eA‘^.  The  Laplacians 
of  the  Csp2-0  bonds  are  not  equal  to  Csp3-0  bonds,  which  indicates  a  strengthening  of  the  bond  to  the 
aromatic  ring.  Whereas  the  ellipticity  of  the  C(3)-0(l)  bond  in  molecule  1  is  a  strong  indicator  of  a  ti 
component  of  the  bond,  that  in  molecule  2  is  closer  to  the  value  observed  for  the  D  ring.  We  attribute 
this  to  the  difference  in  the  hydrogen  bonding  in  these  two  molecules.  The  small  values  of  the  ellipticity 
(0.05  -  0.1)  in  the  Csp3-0  for  both  molecules  indicate  an  unperturbed  pure  a  bond. 

Bond  path:  In  Table  4,  Ad  describes  the  deviation  of  bond  critical  point  (CP)  fi'om  an  intemuclear  axis 
of  a  pair  of  bonded  atoms  in  molecule.  Essentially,  the  bond  path^°’^^  is  a  line  which  links  two  atoms  by 
a  maximum  charge  density.  The  bond  path  length  significantly  differs  fi-om  intemuclear  distance  and  the 
amount  of  deviation  indicates  the  extent  at  which  the  bond  is  strained.  The  distances  di  and  da  (see  Table 
4)  are  the  bond  path  lengths  between  CP  and  the  respective  bonded  atoms.  Their  inequality  implies  the 
degree  of  polarization  of  atom  in  the  bond^^’^^ .  On  inspecting  all  homonuclear  atoms  of  the  ring  systems 
A  through  C,  the  minimum  and  maximum  deviations  are  0.004  and  0.03  A.  These  deviations  are  smaller 
compared  with  the  D-ring  homonuclear  bonds  showing  maximum  deviation  0.044  A.  It  confirms  that  the 
D-ring  bonds  are  strained.  However,  this  deviation  is  smaller  than  the  highly  strained  C-C  bonds  of 
cyclopropane‘s^,  bullvalence'S^  (0.12  A)  and  cubane'S^  (0.17  A)  stmctures.  The  hetronuclear  C-O  bonds 
show  marked  differences  in  both  molecules.  The  CP  of  Csp2-0  bond  in  the  first  molecule  deviates  by 
0.03  A  fi'om  the  intemuclear  axis,  while  in  molecule  2  it  lies  very  close  (0.01  A).  In  molecule  1,  the 
deviation  of  BCP  in  the  bonds  C(17)-0(2)  and  C(16)-0(3)  are  0.03  and  0.01  A  respectively.  But  these 
deviations  are  exactly  opposite  in  the  bonds  C(17’)-0(2')  [0.01  A]  and  C(16')-0(3')  [0.04  A]  of 
molecule  2. 

Topological  properties  of  hydrogen  bonds:  Table  5  reports  the  hydrogen  bond  (O-H— O)  topological 
properties  of  molecules  1  and  2.  A  (3,-1)  hydrogen  bond  critical  point’s  (HCP)  which  are  located  in 
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between  H  and  O  atoms  of  neighboring  molecules.  As  we  mentioned  earlier,  there  are  three  hydrogen 
bonds  exist  in  each  molecule.  Of  the  three  hydrogen  bonds  from  molecule  1,  the  two  bonds 
[0(l)-H(lA)—0(3)  and  0(3)-H(3A)— 0(1')]  display  approximately  linear  hydrogen  bonding,  171.7, 
170.5°  respectively.  The  corresponding  electron  densities  of  these  bonds  at  the  HOP  are  similar  with 
values  of  0.23  and  0.20  eA’^,  respectively.  Both  of  these  bonds  have  a  positive  Laplacian  V^p(rc)  at  the 
HOP.  This  is  characteristic  of  closed  shell  interactions^®  which  are  governed  by  the  contraction  of  the 
electron  density  towards  each  of  the  interacting  atoms.  In  molecule  2,  we  observe  the  similar  electron 
densities  and  Laplacians  for  the  bonds  0(r)-H(TA)—  0(3')  and  0(2')-H(2'A)-0(2').  The  hydrogen 
bond  0(2)-H(2A)—0(l)  in  molecule  1  shows  a  moderate  electron  density  at  the  HCP  (0.11  eA’^)  with  a 
positive  Laplacian.  Although  the  hydrogen  bond  0(2')-H(2'A)— 0(1')  is  approximately  linear  (169.2°), 
the  electron  density  is  very  small  at  the  CP  (0.05  eA'^).  This  may  be  attributed  to  the  longer 
H(2'A)—0(TA)  distance  and  the  fact  that  O(rA)  accepts  one  additional  hydrogen  bond.  Further  more, 
we  have  calculated  bond  path  of  H-0  bonds.  In  H(1A)— 0(3)  bond  the  CP  is  not  much  deviated  (0.056 
A)  from  intemuclear  axis  and  it  is  close  (0.045  A)  to  similar  bond  H(l'A)—0(3')  in  molecule  2.  But  in 
the  16a,  17P  region  of  molecule  1,  the  H(2A)— 0(1)  deviates  significantly  (0.26  A)  compared  with 
H(3A)-0(r)  showing  0.05  A.  Similarly  in  molecule  2,  the  HCP’s  both  for  the  bonds  H(2'A)— O(l')  and 
H(3'A)  "0(2)  are  deviated  by  0.36  and  0.1  A  respectively.  This  strange  deviation  in  17p  regions  of  both 
molecules  indicates  that  these  are  highly  strained  hydrogen  bonds.  Wide-angle  deviations  are  larger  than 
the  reported  value  (0.12  A)  which  was  obtained  from  extensive  electron  density  study  on  0-H— 0 
hydrogen  bonds'*®. 

3.  Electrostatic  potential 

An  additional  important  application  of  the  multipole  model  is  that  it  provides  an  opportxmity  to 
understand  the  chemical  reactivity  of  molecules.  That  is  to  identify  the  places  where  nucleophilic  or 
electrophilic  attack  may  take  place.  The  regions  of  the  positive  potential  will  attract  nucleophilic 
reagents  and  the  electronegative  regions  of  the  molecule  will  determine  the  approach  of  the  electrophilic 
reagents.  The  electrostatic  potential  not  only  plays  an  important  role  in  directing  chemical  attack,  but 
also  in  the  molecular  recognition'*’’'*^  (hormone  receptor,  drug  receptor,  enzyme-substrate  interaction). 
This  is  because  electrostatic  forces  are  long-range  interactions  which  determine  the  most  stable 
geometrical  orientation  of  two  molecules.  To  investigate  the  aspect  of  the  estrogen  hormone  receptor 
interaction,  we  have  calculated  the  electrostatic  potential  around  the  molecules  by  the  method  of  Su  and 
Coppens^®  based  on  the  populations  obtained  from  the  multipole  refinement.  The  advantage  of  this 
method  is  that  it  allows  to  compute  the  ESP  at  points  inside  and  outside  the  Van  der  Waals  surface  of  the 
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charge  distribution^®  and  explore  the  topology  of  V(r)^®’^^  Generally,  intermolecular  potential  involves 
the  effect  of  electrostatic,  induction,  short  range  and  dispersion  interactions^^.  Although  all  these  kinds 
of  interactions  alter  the  intermolecular  interaction,  the  first  one  contributes  significantly^®.  Thus,  the 
molecules  extracted  fi:om  crystal  lattice  account  for  this  many  body  effect.  We  have  mapped  the  ESP  of 
the  isolated  estriol  single  molecule  and  the  fragments  of  the  molecules  which  are  involved  in  hydrogen 
bonding. 

Figures  6a  and  6b  show  contour  plots  (in  the  molecular  mean  plane)  of  the  electrostatic  potential 
of  molecules  1  and  2  isolated  from  the  crystal  lattice.  The  view  of  positive  contours  around  the  aromatic 
rings  are  more  diffuse  than  those  of  the  saturated  rings.  In  molecule  1,  the  positive  contours  are  very 
compact  around  the  atoms  of  the  rings  B,  C,  D  compared  with  molecule  2.  In  Figure  6a,  an  empty  (no 
ESP  contours  inside  the  ring)  view  of  A-ring  of  molecule  1  indicates  that  the  net  electrostatic  effect  is 
zero,  in  comparison  with  the  other  ring  systems  showing  the  moderate  potential  values.  The  negative 
contours  around  the  16a,  ivp  oxygen  atoms  are  much  more  pronounced  than  for  the  A-ring.  The 
potential  aroimd  molecule  1  is  more  symmetric  than  for  molecule  2. 

We  have  computed  a  3D  isosurface  electrostatic  potential  of  the  pseudomolecule  extracted  from 
the  crystal.  In  Figure  7,  the  anti-parallel  view  of  two  molecules  reveals  the  difference  of  their  ESP  in  the 
asymmetric  unit.  The  positive  equipotential  (+0.5  eA"’)  surface  shows  the  extension  of  ESP  bound  to  the 
molecules  1  and  2.  Especially,  the  bulk  surface  bound  by  aromatic  ring  of  molecule  1  confirms  the 
previous  2D  observation.  The  negative  equipotential  surface  indicates  the  amoimt  of  extended  region  for 
the  value  -0.15  eA'*.  More  importantly,  both  the  molecules  have  strong  extended  electronegative  regions 
at  the  16a,  17P-hydroxy  groups,  in  contrast  to  the  less  electronegative  regions  observed  at  the  3-hydroxy 
groups.  The  negative  region  of  molecule  1  exhibit  spherically  bounded  ESP  surface  than  in  molecule  2. 
These  distinct  features  can  be  attributed,  at  least  to  some  extent,  to  the  difference  in  polarization  of  such 
as  weak  and  strong  hydrogen  bonded  atoms.  This  kind  of  isosurface  ESP  analysis  is  becoming  common 
in  the  biological  and  medicinal  important  molecules®"^'®®.  Finally,  the  view  of  the  rich  electronegative 
16a,  17P  hydroxyl  region  of  two  molecules  together  implies  the  diagonal  orientation  of  ESP  field 
vector®^'®®  in  the  molecular  geometry. 

Furthermore,  we  have  added  the  fi'agments  of  both  molecules  involved  in  strong  hydrogen 
bonding  to  visualize  the  polarization.  Figures  8a  and  b  show  the  hydrogen  bond  environment  of  the  O- 
atoms  in  molecules  1  and  2.  Since  ESP  is  attributed  to  the  long-range  electrostatic  interaction,  the 
hydrogen  bonded  atoms  affect  the  existing  potential  around  the  molecule.  The  amount  of  variation  of  the 
ESP  of  the  unperturbed  molecules  depends  upon  the  direction  and  charge  distribution  of  the  two  or  more 
interacting  molecules  or  fi'agments.  The  hydrogen  bonded  fragments  of  the  neighboring  molecules  show 
the  visual  representation  of  polarization  of  hydrogen  bonded  atoms  and  the  net  effect  of  ESP  in  the  inter 
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molecular  region.  Similar  effects  were  observed  in  phosphorylethanolamine^^  crystal  structure.  Figure  8 
emphasizes  the  net  effect  of  ESP  around  each  0-H  groups  of  both  molecules,  after  adding  the  hydrogen 
bonding  fragments. 

4.  Net  atomic  charges 

The  net  atomic  charges  (Table  6a)  have  been  calculated  from  a  K-refinement  (k,  Py,  scale).  In  this 
refinement  the  chemically  equivalent  atoms  were  constrained  as  described  in  the  multipole  refinement, 
and  no  charge  transfer  has  been  allowed  between  the  molecule  1  and  2.  Both  molecules  were  kept 
neutral.  All  oxygen  atoms  carry  significant  negative  charges  in  the  range  -0.63  to  -0.87e.  Similarly  all 
hydroxyl  hydrogen  atoms  carry  high  positive  charges,  approximately  double  the  amoimt  of  charge  of  the 
other  hydrogen  atoms  in  the  molecules.  The  charges  of  chemically  equivalent  aromatic  atoms  C(2)  and 
C(3)  are  equal  and  they  are  much  smaller  than  C(l),  as  expected.  The  equivalent  Csps  atoms  C(6),  C(7), 
C(ll),  and  C(12)  are  carrying  moderate  charges  and  they  are  almost  equal.  A  similar  Csps  (CH)  kind  of 
atoms  C(8),  C(9)  and  C(14)  carrying  small  amoimt  of  charge  except  C(9).  These  charges  are  effectively 
larger  than  for  the  similar  D-ring  atoms  C(16)  and  C(17)  (see  Table  6),  this  anomaly  may  be  due  to  the 
different  environment.  The  overall  charge  distribution  of  C-  atoms  in  the  ring  systems  A  through  D 
indicates  that  the  aromatic  A-ring  and  D-ring  are  holding  lesser  amount  of  charge  compared  to  B  and  C 
ring  S3^tems.  A  group  charge  calculation  (Table  6b)  on  both  molecules  indicates  that  the  central  regions 
are  close  to  neutral,  whereas  the  head  and  tails  form  charged  units.  Surprisingly,  the  3-hydroxy 
(molecule  1)  and  16a  (molecule  2)  regions  are  holding  less  negative  charges  (see  Table  6b), 

Conclusion 

The  present  study  on  estriol  gives  an  information  on  both  the  structural  and  electronic  properties 
of  the  molecule.  The  differences  in  hydrogen  bonding,  molecular  packing  environment,  and 
interamolecular  steric  effects  are  responsible  for  the  small  differences  in  the  two  molecules  in  the 
asymmetric  unit.  The  significant  differences  reported  for  the  room  temperature  structure  were  not 
confirmed  here.  The  comparison  of  torsion  angles  between  the  two  molecules  shows  that  the  major 
difference  at  the  head  and  tail  portions  of  the  molecules  may  be  attributed  to  the  different  hydrogen 
bonding  environment,  both  molecules  being  involved  in  the  strong  hydrogen  bonding. 

The  electronic  information  of  A-ring  system  of  the  estrogen  molecules  is  very  important  for  the 
ligand  receptor  binding^^’®.  The  electron  densities  at  the  CP  of  the  aromatic  C-C  bonds  are  larger  than 
for  the  Csp3-Csp3  bonds  of  the  B,  C  and  D  ring  systems.  The  additional  electron  density  accumulation  in 
the  aromatic  bonds  along  with  the  ellipticity  indicates  a  n  contribution.  The  unconstrained  refinement  of 
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the  hydroxyl  groups  between  the  two  molecules  shows  differences  in  both  C-0  and  0-H  bond  electron 
densities  at  the  CP.  The  average  Csp3-0  bond  density  is  significantly  less  than  Csp2-0  density.  The 
charge  accumulation  in  the  0-H  (16a,17P  position)  bond  is  slightly  larger  than  in  the  3 -hydroxy  0-H 
bond.  Strong  electronegative  regions  were  found  in  the  electrostatic  potential  maps  at  the  vicinity  of  the 
oxygen  atoms  of  the  16a,17P  groups  of  both  molecules.  The  extension  of  positive  and  negative  contours 
explicitly  indicates  the  electrostatic  interaction  with  the  neighboring  molecules. 

The  hydrogen  bonding  O-H—O  and  O— H-0  interaction  at  the  head  and  tail  portion  of  the 
molecule  shows  the  possibility  of  having  the  positive  and  negative  regions  of  the  ligand  binding  domain 
of  the  receptor.  The  favored  direction  of  alignment  of  the  non-identical  molecules  together  in  the  crystal 
lattice  is  along  the  direction  of  maximum  negative  ESP.  The  large  electron  density  at  the  hydrogen  bond 
CP’s  encourages  the  very  strong  nature  of  the  hydrogen  bonding  in  the  estriol  crystal. 

Net  atomic  charges  computed  fi'om  the  K-refinement  gave  a  reasonable  distribution  of  charges.  In 
the  16a,17P  position,  the  oxygen  atoms  0(2),  0(3)  (molecule  1)  and  0(2'),  0(3')  (molecule  2)  carry 
more  negative  charges  than  3-hydroxy  atoms  0(1)  and  0(T). 

To  relate  the  results  obtained  fi"om  the  electron  density  analysis  and  ESP  with  biological 
characteristics  of  estriol,  the  small  negative  ESP  at  the  3-hydroxy  positions  and  the  broad  shape  of 
negative  ESP  at  the  16a,17P  around  0-atoms  overlapping  with  the  neighboring  groups,  and  the 
conformational  flexibility  (i.e.,  two  molecules  arises  in  the  asymmetric  unit  with  different 
conformations)  reduces  the  effective  binding  of  molecules  with  the  receptor.  The  first  two  of  these 
results  support  the  conclusion  derived  in  isopotential  study' ^  in  which  the  strong  extended  ESP  at  the  3- 
hydroxyl  region  is  expected. 

Thus,  the  way  of  representing  and  analyzing  the  experimental  ESP  brings  out  the  approximate 
electrostatic  characteristics  of  a  molecule  in  the  nonbonding  environment.  And  the  ESP  extracted  from 
the  crystal  lattice  is  suitable  for  the  simulation  of  hormone-receptor  interactions  at  different 
physiological  situations. 
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Table  1.  Crystal  data  and  Experimental  conditions^ 

Crystal  Data 


Empirical  formula 

C,8H2403 

Formula  weight 

288.37 

Crystal  system 

Monoclinic 

Space  group 

P2, 

a  (A) 

7.5077(4) 

b(A) 

23.0809(12) 

c  (A) 

9.1632(5) 

P  (deg) 

110.957(2) 

V(A') 

1482.8 

Z 

4 

F(OOO) 

624 

Dx(Mg/m^) 

1.292 

p,  (mm"') 

0.055 

Data  collection 

Diffractometer 

Brukers  SMART  2K  CCD 

Radiation,  X  (A) 

Ag  Ka,  0.56087 

Scan  method,  width  (deg) 

(O-scan,  0.3 

No.  of  reflections  measured 

278184 

sin(0/X,)tnax  (A  ) 

1.29 

Rint 

0.042 

No.  of  unique  data 

43918 

Multipole  refinement 

No.  of  relections  (F  >  3cr,  sm(0/X)max=  1-1  A"') 

25520 

R(F),  wR(F) 

0.032,  0.035 

GOF 

0.98 

Nobs/Npar 

26.6 

‘Rint=  Ill-Imeanl/B;  R(F)=I(  I  lF„l-lFel  I /SIf„I);  wR(F)=(S(F„-Fe)7lwFo2)’''^;w=l/a\Fo) 
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Table  2.  Selected  geometrical  parameters 
of  molecule  1  and 

Bond  lengths  (A) _ 

bond  ml  m2 


C(l)-C(2)  1.3938(8)  1.3962(7) 

C(2)-C(3)  1.3952(8)  1.3936(7) 

C(3)-C(4)  1.3939(7)  1.3921(7) 

C(4)-C(5)  1.3991(7)  1.3978(6) 

C(5)-C(10)  1.4084(7)  1.4111(7) 

C(l)-C(10)  1.4064(7)  1.4034(7) 

0(1)-C(3)  1.3693(8)  1.3785(7) 

C(5)-C(6)  1.5146(7)  1.5154(7) 

C(6)-C(7)  1.5255(7)  1.5250(7) 

C(7)-C(8)  1.5293(7)  1.5264(7) 

C(8)-C(9)  1.5444(7)  1.5424(6) 

C(9)-C(10)  1.5278(7)  1.5253(6) 

C(9)-C(ll)  1.5388(7)  1.5408(7) 

C(ll)-C(12)  1.5409(7)  1.5441(7) 
C(12)-C(13)  1.5325(7)  1.5313  (7) 
C(13)-C(14)  1.5508(7)  1.5415(7) 
C(13)-C(17)  1.5360(7)  1.5462(7) 
C(13)-C(18)  1.5368(8)  1.5334(7) 
C(14)-C(15)  1.5400(7)  1.5360(7) 
C(14)-C(8)  1.5235(7)  1.5172(7) 

C(15)-C(16)  1.5529(7)  1.5499(7) 
C(16)-(17)  1.5413(7)  1.5463(7) 

C(16)-0(3)  1.4387(9)  1.4278(8) 


C(17)-0(2)  1.4293(7)  1.4168(8) 


*/M7-molecule  1  and  w2-molecule  2 


Bond  angles  (deg) 


Bond 

ml 

m2 

C(l)-C(2)-C(3) 

119.3(1) 

199.1(1) 

0(1)-C(3)-C(2) 

118.4(1) 

118.1(1) 

0(1)-C(3)-C(4) 

121.9(1) 

121.8(1) 

C(2)-C(3)-C(4) 

119.6(1) 

120.0(1) 

C(3)-C(4)-C(5) 

121.0(1) 

120.6(1) 

C(4)-C(5)-C(6) 

118.1(1) 

117.2(1) 

C(4)-C(5)-C(10) 

120.3(1) 

120.4(1) 

C(6)-C(5)-C(10) 

121.6(1) 

122.4(1) 

C(5)-C(6)-C(7) 

113.1(1) 

113.6(1) 

C(6)-C(7)-C(8) 

110.2(1) 

110.7(1) 

C(7)-C(8)-C(9) 

109.0(1) 

108.4(1) 

C(7)-C(8)-C(14) 

112.7(1) 

111.9(1) 

C(9)-C(8)-C(14) 

107.8(1) 

109.4(1) 

C(8)-C(9)-C(10) 

112.2(1) 

109.9(1) 

C(8)-C(9)-C(ll) 

111.3(1) 

112.5(1) 

C(10)-C(9)-C(ll) 

114.8(1) 

114.3(1) 

C(l)-C(10)-C(5) 

117.6(1) 

117.7(1) 

C(l)-C(10)-C(9) 

121.1(1) 

121.8(1) 

C(5)-C(10)-C(9) 

121.0(1) 

120.5(1) 

C(9)-C(ll)-C(12) 

111.9(1) 

112.7(1) 

C(ll)-C(12)-C(13) 

111.4(1) 

110.6(1) 

C(12)-C(13)-C(14) 

108.5(1) 

108.4(1) 

C(12)-C(13)-C(17) 

115.4(1) 

117.1(1) 

C(12)-C(13)-C(18) 

110.6(1) 

109.8(1) 

C(14)-C(13)-C(17) 

97.6(1) 

98.6(1) 

C(14)-C(13)-C(18) 

113.6(1) 

113.5(1) 

C(17)-C(13)-C(18) 

110.7(1) 

109.1(1) 

C(8)-C(14)-C(13) 

113.7(1) 

113.8(1) 

C(8)-C(14)-C(15) 

120.6(1) 

118.1(1) 

C(13)-C(14)-C(15) 

103.5(1) 

103.8(1) 

C(14)-C(15)-C(16) 

104.0(1) 

104.2(1) 

0(3)-C(16)-C(15) 

113.9(1) 

113.1(1) 

0(3)-C(16)-C(17) 

111.6(1) 

113.0(1) 

C(15)-C(16)-C(17) 

104.8(1) 

105.6(1) 

0(2)-C(17)-C(13) 

119.7(1) 

116.6(1) 

0(2)-C(17)-C(16) 

112.0(1) 

110.2(1) 

C(13  )-C(17)-C(16) 

104.6(1) 

103.2(1) 

C(3)-0(1)-H(1A) 

110.4(2) 

108.4(2) 

C(17)-0(2)  -H(2A) 

112.7(2) 

103.8(2) 

Ca6VOf31-H13A') 

116.9fn 

mssmrn 

21 


Torsion  angles  (deg) 


Bond 

ml 

m2 

A-ring 

C(10)-C(l)-C(2)-C(3) 

-0.5(2) 

-0.5(2) 

C(l)-C(2)-C(3)-C(4) 

2.1(2) 

0.8(2) 

C(2)-C(3)-C(4)-C(5) 

-2.1(2) 

-0.3(2) 

C(3)-C(4)-C(5)-C(10) 

0.5(2) 

-0.4(2) 

C(4)-C(5)-C(10)-C(l) 

1.1(2) 

0.7(2) 

C(5)-C(10)-C(l)-C(2) 

-1.1(2) 

-0.3(2) 

B-ring 

C(10)-C(5)-C(6)-C(7) 

-19.3(2) 

-10.5(2) 

C(5)-C(6)-C(7)-C(8) 

48.4(2) 

41.1(2) 

C(6)-C(7)-C(8)-C(9) 

-65.0(2) 

-66.0(2) 

C(7)-C(8)-C(9)-C(10) 

50.7(2) 

57.8(2) 

C(8)-C(9)-C(10)-C(5) 

-22.2(2) 

-27.9(2) 

C(9)-C(10)-C(5)-C(6) 

6.2(2) 

4.3(2) 

C-ring 

C(8)-C(9)-C(ll)-C(12) 

56.2(2) 

51.8(2) 

C(9)-C(ll)-C(12)-C(13) 

-55.2(2) 

-54.7(2) 

C(ll)-C(12)-C(13)-C(14) 

54.4(2) 

56.8(2) 

C(12)-C(13)-C(14)-C(8) 

-58.7(2) 

-60.2(2) 

C(13)-C(14)-C(8)-C(9) 

59.3(2) 

56.9(2) 

C(14)-C(8)-C(9)-C(ll) 

-56.6(2) 

-51.4(2) 

D-ring 

C(13)-C(14)-C(15)-C(16) 

-33.3(2) 

-31.2(2) 

C(14)-C(15)-C(16)-C(17) 

4.2(2) 

2.0(2) 

C(15)-C(16)-C(17)-C(13) 

26.7(2) 

27.6(2) 

C(16)-C(17)-C(13)-C(14) 

-45.9(2) 

-45.8(2) 

C(17)-C(13)-C(14)-C(15) 

48.6(2) 

47.7(2) 

C(13)-C(17)-C(16)-0(3) 

150.4(2) 

151.8(2) 

C(15)-C(16)-C(17)-0(2) 

157.7(2) 

152.9(2) 

C(13)-C(17)-0(2)-H(2A) 

-4.7(2) 

-78.1(2) 

C(15)-C(16)-0(3)-H(3A) 

84.4(2) 

50.8(2) 

C(2)-C(3)-0(1)-H(1A) 

169.7(2) 

-176.3(2) 

C(4)-C(3)-0(1)-H(1A) 

-8.5(2) 

3.6(2) 
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Table  3.  Hydrogen  bonds* 


bond 

0-H  (A) 

0 -  O  (A) 

H  -O  (A) 

O-H-0  (deg) 

ml 

0(l)-H(lA)-0(3)* 

0.97 

2.646(2) 

1.682(2) 

171.7(1) 

0(2)-H(2A)-0(l)^ 

0.97 

1.851(1) 

158.6(1) 

0(3)-H(3A)  -OClT 
m2 

0(l’)-H(rA)-0(3y 

0.97 

2.890(1) 

1.930(1) 

170.5(1) 

0.97 

2.640(2) 

1.684(1) 

169.9(1) 

0(2')-H(2'A)-0(iy 

0.97 

2.884(1) 

1.927(1) 

169.2(1) 

0(3')-H(3’A)-0(2/ 

0.97 

2.736(1) 

1.857(1) 

149.9(1) 

“Symmetry  codes:  (a)  -x+2,  y-!4, -z+1;  (b)  -x+1,  y+Vi, -z+1;  (c)  x,  y,  z-1; 
(d)-x+l,y+'/4, -z+2;  (e) -x+2,  y-'/2,  -z+2;  (f)  -x+1,  y-Vi, -z+1 
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Table  4.  Analysis  of  the  properties  of  the  bond  critical  points  for  the  molecule  1 


(first  row)  and  2  (second  row)* 


bond 

>.1 

X2 

X3 

p(rc)(eA'^) 

V'p(rc)(eA-*) 

Rij(A) 

d.(A) 

da  (A) 

8 

Ad  (A) 

A-ring 

C(l)-C(2) 

-15.07 

-12.99 

9.95 

2.14(1) 

-18.1(1) 

1.3939 

0.7040 

0.6899 

0.16 

0.01 

C(2>-C(3) 

-16.61 

-14.18 

8.67 

2.19(2) 

-22.1(1) 

1.3967 

0.6846 

0.7121 

0.17 

0.03 

C(3)-C(4) 

-15.88 

-13.36 

8.29 

2.15(1) 

-21.0(1) 

1.3947 

0.7607 

0.6340 

0.19 

0.02 

C(4)-C(5) 

-14.72 

-12.60 

9.51 

2.04(1) 

-17.8(1) 

1.3993 

0.6975 

0.7017 

0.17 

0.01 

C{5)-C(10) 

-14.47 

-12,27 

9.74 

2.05(1) 

-17.0(1) 

1.4087 

0.7253 

0.6834 

0.18 

0.01 

C(10)-C(l) 

-13.89 

-11.35 

10.12 

2.04(1) 

-15.1(1) 

1.4065 

0.7349 

0.6716 

0.22 

0.01 

C(l)-H(l) 

-15.65 

-14.68 

16.05 

1.80(1) 

-14.3(1) 

1.0817 

0.7067 

0.3750 

0.07 

0.01 

C(2)-H(2) 

-17.65 

-15.99 

11.59 

1.90(1) 

-22.1(1) 

1.0832 

0.6889 

0.3943 

0.10 

0.02 

C(4)-H(4) 

-16.16 

-15.03 

11.97 

1.78(4) 

-19.2(1) 

1.0849 

0.6950 

0.3900 

0.07 

0.002 

C(3)-0(l) 

-18.17 

-14.42 

11.28 

2.05(3) 

-21.3(1) 

1.3708 

0.8343 

0.5365 

0.26 

0.03 

-17.65 

-16.46 

12.08 

2.16(2) 

-22.0(1) 

1.3787 

0.8627 

0.5159 

0.07 

0.01 

0(1)-H(1A) 

-43.38 

-42.37 

27.52 

2.08(8) 

-58.2(8) 

0.9700 

0.8205 

0.1495 

0.02 

0.004 

-40.29 

-39.52 

31.58 

2.18(8) 

-48.2(7) 

0.9661 

0.7967 

0.1694 

0.02 

0.004 

B-ring 

C(5>-C(6) 

-11.20 

-10.55 

10.64 

1.71(1) 

-11.1(1) 

1.5148 

0.7741 

0.7408 

0.06 

0.01 

C(6)-C(7) 

-11.48 

-10.65 

10.97 

1.74(1) 

-11.2(1) 

1.5258 

0.7633 

0.7625 

0.08 

0.02 

C(7)-C(8) 

-11.20 

-10.26 

10.40 

1.67(1) 

-11.1(1) 

1.5293 

0.7806 

0.7487 

0.09 

0.004 

C(8)-C(9) 

-11.94 

-11.08 

9.82 

1.71(1) 

-13.2(1) 

1.5453 

0.7704 

0.7749 

0.08 

0.03 

C(9)-C(10) 

-11.36 

-10.41 

10.20 

1.69(1) 

-11.6(1) 

1.5284 

0.7911 

0.7373 

0.09 

0.02 

C(6)-H(6A) 

-15.34 

-13.61 

12.99 

1.82(4) 

-16.0(1) 

1.0944 

0.6649 

0.4295 

0.13 

0.01 

C(6)-H(6B) 

-13.83 

-13.29 

18.21 

1.64(4) 

-8.9(2) 

1.0919 

0.7474 

0.3446 

0.04 

0.02 

C(7)-H(7A) 

-17.50 

-16.33 

15.29 

1.86(4) 

-18.6(1) 

1.0933 

0.7508 

0.3425 

0.07 

0.004 

C(7)-H(7B) 

-16.31 

-15.14 

14.93 

1.83(4) 

-16.5(1) 

1.0922 

0.7205 

0.3717 

0.08 

0.01 

C(8)-H(8) 

-16.79 

-15.76 

15.69 

1.76(4) 

-16.9(2) 

1.0990 

0.7745 

0.3245 

0.07 

0,01 

C(9)-H(9) 

-16.14 

-15.65 

14.21 

1.81(4) 

-17.6(1) 

1.0984 

0.7295 

0.3688 

0.0 

0.01 

C-ring 

C(9>-C(ll) 

-10.75 

-10.54 

10.25 

1.65(2) 

-11.0(1) 

1.5389 

0.7555 

0.7834 

0.02 

0.01 

C(11)~C(12) 

-10.95 

-9.56 

10.61 

1.63(1) 

-9.9(1) 

1.5410 

0.7595 

0.7816 

0.14 

0.01 

C(12)-C(13) 

-12.33 

-10.71 

9.98 

1.73(1) 

-13.1(1) 

1.5328 

0.7586 

0.7742 

0.15 

0.02 

C(13)-C(14) 

-11.29 

-10.42 

9.46 

1.62(1) 

-12.3(1) 

1.5510 

0.7749 

0.7761 

0.08 

0.01 

C(14)-C(8) 

-12.98 

-11.80 

8.96 

1.73(1) 

-15.8(1) 

1.5236 

0.7627 

0.7609 

0.10 

O.Ol 

C(11)-H(11A) 

-14.98 

-14.54 

15.54 

1.75(4) 

-14.0(1) 

1,0938 

0,7228 

0.3710 

0.03 

0.02 

C(11)-H(11B) 

-16.06 

-15.38 

14.43 

1.80(4) 

-17.0(1) 

1.0909 

0.7187 

0,3722 

0.04 

0.01 

C(12>-H(12A) 

-15.92 

-15.31 

17.02 

1.77(4) 

-14.2(2) 

1.0917 

0.7602 

0.3315 

0.04 

0.004 

C(12)-H(12B) 

-15.35 

-14.45 

13.31 

1.79(4) 

-16.5(1) 

1.0901 

0.6871 

0.4031 

0.06 

0.01 

C(14)-H(14) 

-16.50 

-16.13 

15.31 

1.75(4) 

-17.3(2) 

1.1011 

0.7792 

0.3219 

0.02 

0.02 

D-ring 

C(14)-C(15) 

-10.27 

-10.19 

10.12 

1.57(1) 

-10.4(1) 

1.5409 

0.7816 

0.7593 

0.01 

0.03 

C(15)-C(16) 

-12.10 

-11.18 

9.38 

1.71(1) 

-13.9(1) 

1.5544 

0.7375 

0.8169 

0.08 

0.035 

C(16)-C(17) 

-12.71 

-12.07 

9.51 

1.79(1) 

-15.3(1) 

1.5416 

0.7916 

0.7499 

0.05 

0.014 

C(17>-C(13) 

-11.13 

-10.28 

9.56 

1.65(1) 

-11.9(1) 

1.5386 

0.7199 

0.8187 

0.08 

0.044 

C(15)-H(15A) 

-17.40 

-15.51 

12.54 

1.82(4) 

-20.4(1) 

1.0927 

0.7242 

0,3684 

0.12 

0.01 

C(15>-H(15B) 

-16.14 

-15.14 

16.09 

1.68(5) 

-15.2(2) 

1.0948 

0.7847 

0.3101 

0.07 

0.02 

C(16>-H(16) 

-19.35 

-18.60 

13.23 

1.94(5) 

-24.7(2) 

1.0968 

0.7741 

0.3227 

0.04 

0.01 

C(17)-H(17) 

-17.41 

-16.91 

15.27 

1.78(4) 

-19.1(2) 

1.1009 

0.7915 

0.3094 

0.03 

0.014 

C(17)-0(2) 

-14.85 

-14.05 

13,75 

1.96(2) 

-15.2(1) 

1.4308 

0.8536 

0.5772 

0.06 

0.03 

-13.95 

-12.67 

14.94 

1.91(3) 

-11.7(1) 

1.4170 

0.8319 

0.5852 

0.10 

0.01 

0(2)-H(2A) 

-47.97 

-45.24 

24.28 

2.26(6) 

-68.9(7) 

0.9665 

0.8117 

0.1548 

0.06 

0.01 

-46.44 

-43.87 

23.04 

2.34(8) 

-67.3(7) 

0.9683 

0.8000 

0.1683 

0.06 

0.005 

C(16)-0(3) 

-14.32 

-13.65 

14.11 

1.87(2) 

-13.9(6) 

1.4389 

0.8588 

0.5800 

0.05 

0.01 

-13.64 

-12.59 

14.89 

1.83(2) 

-11.3(1) 

1.4301 

0.8364 

0.5937 

0.08 

0.04 

0(3)-H(3A) 

-35.31 

-34.33 

40.97 

2.34(7) 

-28.7(5) 

0.9694 

0.7599 

0.2095 

0.03 

0.01 

-41.18 

-39.53 

34.28 

2.22(7) 

-46.4(6) 

0.9680 

0.7983 

0.1698 

0.04 

0.02 

C(13)-C(18) 

-11.30 

-10.73 

10.44 

1.71(1) 

-11.6(1) 

1.5369 

0.7892 

0.7477 

0.05 

0.01 

C(18)-H(18A) 

-16.71 

-15.95 

15.99 

1.89(5) 

-16.7(2) 

1.0629 

0.7085 

0.3544 

0,05 

0,004 

C(18)-H(18B) 

-16.28 

-15.68 

13.80 

1.88(5) 

-18.2(2) 

1.0549 

0.6779 

0.3771 

0.04 

0.01 

C(18^H(18C) 

-16.05 

-14.81 

15.84 

1.85(5) 

-15.0(1) 

1.0577 

0.6869 

0.3708 

0.08 

0.015 

*  Xu  X2  and  X3  are  the  eigenvalues  of  Hessian  matix.  p  (Tc)  is  the  electron  density;  V^p  (fc)  is  the  Laplacian  at  the  CP;  di  and 
d2  are  the  distances  from  the  critical  point;  Ry  is  the  sum  of  di  &  d2 ;  e  is  the  bond  ellipticity.  Ad  is  the  deviation  of  CP 
from  the  interauclear  axis 
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Table  5.  Topological  Characteristics  of  the  electron  density  at  the  hydrogen 
bondcritical  points* _ 


bond 

p(rc)(eA‘^) 

V^p(rc)(eA'^) 

Rii  (A) 

d,(A) 

da  (A) 

ml 

0(1)-H(1A)  -OCS)® 

0.23(5) 

5.07(7) 

1.6869 

1.1526 

0.5343 

0(2)-H(2A)-0(l)’’ 

0.11(2) 

3.00(1) 

1.9287 

1.2551 

0.6737 

0(3)-H(3A)-0(lT 

0.20(3) 

3.00(2) 

1.9328 

1.2381 

0.6947 

m2 

0(l')-H(l'A)-0(3')‘' 

0.24(5) 

5.02(5) 

1.6863 

1,1299 

0.5664 

0(2')-H(2'A)-0(lT 

0.05(2) 

2.31(1) 

2.0748 

1.3544 

0.7204 

0(2')-H(2'A)-0(2)^ 

0.20(3) 

3.41(2) 

1.8682 

1.2186 

0.6496 

*  Symmetry  codes:  (a)  -x+2,  y-Vi,  -z+1;  (b)  -x+1,  y+Vi,  -z+1;  (c)  x,  y,  z-1 ;  (d)  -x+1,  y+!4,  -z+2; 

(e)  -x+2,  y-Vi,  -z+2;  (f)  -x+1,  y-14,  -z+1 
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Table  6a.  Net  atomic  charges  (e)  for  the 
molecule  1  (first  row)  and  2  (second  row) 


Atom  charge  Atom  charge 


C(l) 

-0.21(3) 

H(6A) 

0.20(2) 

C(2) 

-0.06(3) 

H(6B) 

0.12(2) 

C(3) 

0.16(3) 

H(7A) 

0.20(2) 

C(4) 

-0.06(3) 

H(7B) 

0.24(2) 

C(5) 

0.00 

H(8) 

0.16(2) 

C(10) 

-0.04(2) 

H(9) 

0.20(2) 

C(6) 

-0.36(3) 

H(llA) 

0.20(2) 

C(7) 

-0.48(3) 

H(llB) 

0.22(2) 

C(8) 

-0.12(3) 

H(12A) 

0.08(2) 

C(9) 

-0.29(3) 

H(12B) 

0.26(2) 

C(ll) 

-0.31(3) 

H(14) 

0.17(2) 

C(12) 

-0.47(3) 

H(15A) 

0.30(2) 

C(13) 

0.00 

H(15B) 

0.21(2) 

C(14) 

-0.07(3) 

H(16) 

0.21(2) 

C(15) 

-0.24(3) 

H(17) 

0.16(2) 

C(16) 

0.05(3) 

H(18A) 

0.20(2) 

C(17) 

-0.08(3) 

H(18B) 

0.27(2) 

C(18) 

-0.61(3) 

H(18C) 

0.19(2) 

0(1) 

-0.63(3) 

H(1A) 

0.43(3) 

-0.79(3) 

0.35(3) 

0(2) 

-0.86(3) 

H(2A) 

0.48(2) 

-0.87(3) 

0.46(2) 

0(3) 

-0.84(3) 

H(3A) 

0.25(2) 

-0.77(3) 

0.46(2) 

H(l) 

0.16(2) 

H(2) 

0.29(2) 

H(4) 

0.31(2) 
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Table  6b.  Group  charge(e)  calculation  for  the 
molecule  1  (first  row)  and  2  (second  row) 


Atoms 

charge(e) 

C(3),  0(1),H(1A) 

-0.04 

-0.20 

C(17),H(17),  0(2A),H(2A) 

-0.30 

-0.33 

C(16),  H(16),  0(3A),  H(3A) 

-0.33 

-0.05 

C(1),H(1),C(10) 

-0.05 

C(2),  H(2) 

0.23 

C(4),  H(4) 

0.28 

C(6),  C(7),  H(6A),  H(6B) 

-0.08 

H(7A),  H(7B) 

C(ll),  C(12),H(11A),H(11B) 

-0.02 

H(12A),  H(12B) 

C(14),  H(14),  C(15),  H(15A),  H(15B) 

0.37 

j  2.  ORTEP  view  of  two  molecules  in  the  as)mimetric  unit  showing 
75%  probability  ellipsoids  of  the  atoms  and  atomic  labeling. 


Figure  4a.  Residual  density  maps  for  molecule  1  in  the  plane  of  C(l),  C(5)  and 
C(16) .  Contour  interval  0.05  eA'^.  Positive,  negative  and  zero  contours  are  drawn  as 

solid,  dotted  and  broken  lines  respectively. 
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Figure  4b.  Residual  density  maps  for  molecule  2  in  the  0(1')  C(5')  C(16')  plane. 
Contour  interval  0.05  eA'^.  Positive,  negative  and  zero  contours  are  drawn  as  solid, 

dotted  and  broken  lines  respectively. 
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Figure  5.  Static  deformation  density  map  of  molecule  1.  Contour  intervals  0.05  eA' 
Positive  solid  line,  negative  dotted  line  and  zero  contour  broken  line. 


Figure  6a.  Electrostatic  potential  contour  map  of  molecule  1  in  the  plane  of  0(1), 
C(5)  and  C(1 6).  Contour  interval  0.05  eA'^;  Positive  contour  is  solid  line, 
negative  contour  is  dotted  line  and  zero  contour  is  broken  line. 


Figure  6b.  Electrostatic  potential  contour  map  of  molecule  2  in  the  plane  of  0(1'), 
C(5')  and  C(16').  Contour  interval  0.05  eA*^;  Positive  contour  is  solid  line, 
negative  contour  is  dotted  line  and  zero  contour  is  broken  line. 
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Figure  7.  Isosurface  representation  of  electrostatic  potential  of  molecule  1  (right) 
and  2  (left).  Surfaces  are  drawn  at  constant  positive  and  negative  potentials  of 
0.5  /e/A  (red)  and  -0.1  e/A  (green).  Grey  dots  represent  zero  surface. 
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(b) 


Figure  8.  The  electrostatic  potential  in  the  hydrogen  bonding  environment  of  oxygen 
atoms  of  (a)  molecule  1  and  (b)  molecule  2.  Contour  interval  0.05  eA'^; 
Positive  solid  lines,  negative  dotted  lines  and  zero  contours  are  broken  lines. 
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Supplementary  Tables 
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Table  SI.  Fractional  coordinates,  equivalent  displacement  parameters  (A^) 
after  multipole  refinement* 


Atom 

X 

y 

z 

Ue, 

Molecule  1 
0(1) 

0.51768(10) 

0.85352(11) 

0.52120(11) 

0.018 

0(2) 

0.79921(9) 

1.32161(11) 

0.40944(9) 

0.015 

0(3) 

1.19065(9) 

1.28089(11) 

0.40369(9) 

0.015 

C(l) 

0.45782(8) 

1.00729(11) 

0.42503(8) 

0.014 

C(2) 

0.41442(8) 

0.94969(11) 

0.44370(8) 

0.016 

C(3) 

0.56255(8) 

0.91059(11) 

0.51231(7) 

0.012 

C(4) 

0.75074(7) 

0.93004(11) 

0.56560(6) 

0.011 

C(5) 

0.79388(7) 

0.98789(11) 

0.54672(6) 

0.010 

C(6) 

1.00170(8) 

1.00579(11) 

0.60531(7) 

0.013 

0(7) 

1.03656(7) 

1.06192(11) 

0.53178(7) 

0.013 

0(8) 

0.89489(7) 

1.10844(11) 

0.53900(6) 

0.010 

0(9) 

0.69131(7) 

1.08931(11) 

0.43663(6) 

0.010 

0(10) 

0.64664(7) 

1.02772(11) 

0.47414(6) 

0.010 

0(11) 

0.54265(8) 

1.13507(11) 

0.43759(7) 

0.013 

0(12) 

0.58776(7) 

1.19493(11) 

0.38487(7) 

0.012 

0(13) 

0.78969(7) 

1.21485(11) 

0.48383(6) 

0.010 

0(14) 

0.93243(7) 

1.16700(11) 

0.47888(6) 

0.010 

0(15) 

1.12818(8) 

1.19665(11) 

0.55149(8) 

0.015 

0(16) 

1.08802(7) 

1.26090(11) 

0.49990(7) 

0.012 

0(17) 

0.87053(7) 

1.26401(11) 

0.41297(6) 

0.010 

0(18) 

0.80161(10) 

1.23048(11) 

0.65010(7) 

0.017 

H(1A) 

0.63209 

0.82990 

0.55330 

0.034 

H(2A) 

0.69428 

1.32373 

0.44686 

0.012 

H(3A) 

1.13395 

1.27295 

0.29273 

0.025 

H(l) 

0.34450 

1.03712 

0.36414 

0.029 

H(2) 

0.26604 

0.93672 

0.39812 

0.044 

H(4) 

0.86510 

0.89941 

0.61941 

0.021 

H(6A) 

1.09283 

0.97171 

0.58700 

0.037 

H(6B) 

1.04093 

1.01387 

0.73024 

0.020 

H(7A) 

1.01476 

1.05396 

0.40896 

0.015 

H(7B) 

1.18207 

1.07730 

0.59177 

0.021 

H(8) 

0.89840 

1.11151 

0.65964 

0.017 

H(9) 

0.68554 

1.08694 

0.31532 

0.016 

H(11A) 

0.40304 

1.12061 

0.35796 

0.024 

H(llB) 

0.53292 

1.13757 

0.55336 

0.020 

H(12A) 

0.57626 

1.19302 

0.26278 

0.018 

H(12B) 

0.48167 

1.22649 

0.38902 

0.019 

H(14) 

0.90581 

1.15747 

0.35520 

0.014 

H(15A) 

1.23477 

1.17769 

0.51002 

0.022 

H(15B) 

1.18913 

1.19417 

0.67908 

0.024 

H(16) 

1.12808 

1.29047 

0.60016 

0.025 

H(17) 

0.84636 

1.25096 

0.29187 

0.011 

H(18A) 

0.76328 

1.19583 

0.70939 

0.033 

H(18B) 

0.70167 

1.26365 

0.64136 

0.026 

H(18C) 

0.94099 

1.24472 

0.71651 

0.025 

Molecule  2 
O(l’) 

0.98786(9) 

1.26560(11) 

1.07236(8) 

0.013 

0(2’) 

0.74949(10) 

0.79556(11) 

1.08083(9) 

0.016 
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0(3') 

0.31379(9) 

0.81942(11) 

0.90942(8) 

0.015 

C(l’) 

1.04605(7) 

1.10997(11) 

1.03189(7) 

0.013 

C(2') 

1.08893(7) 

1.16877(11) 

1.05868(7) 

0.013 

C(3') 

0.94251(7) 

1.20768(11) 

1.04700(6) 

0.010 

C(4') 

0.75670(7) 

1.18771(11) 

1.01090(6) 

0.010 

C(5') 

0.71459(7) 

1.12874(11) 

0.98527(6) 

0.009 

C(6’) 

0.51022(7) 

1.11043(11) 

0.95313(8) 

0.013 

C(7') 

0.47065(7) 

1.04728(11) 

0.90310(7) 

0.012 

C(8') 

0.63563(7) 

1.00865(11) 

0.99989(6) 

0.009 

C(9') 

0.81329(7) 

1.02453(11) 

0.96109(6) 

0.009 

C(IO') 

0.86001(7) 

1.08860(11) 

0.99460(6) 

0.010 

cxii*) 

0.98227(7) 

0.98340(11) 

1.03926(7) 

0.012 

C(12’) 

0.92732(7) 

0.91892(11) 

1.00723(7) 

0.012 

C(13') 

0.75705(7) 

0.90425(11) 

1.05525(6) 

0.009 

C(14') 

0.59061(7) 

0.94487(11) 

0.96631(6) 

0.009 

C(15') 

0.41989(8) 

0.91948(11) 

0.99963(8) 

0.014 

C(16-) 

0.46119(8) 

0.85355(11) 

1.01718(7) 

0.012 

C(17’) 

0.65666(8) 

0.84557(11) 

0.99892(6) 

0.011 

C(18’) 

0.81622(9) 

0.90814(11) 

1.23340(6) 

0.015 

H(l'A) 

0.87426 

1.28653 

1.06698 

0.027 

H(2'A) 

0.84956 

0.78829 

1.03968 

0.022 

H(3A) 

0.28464 

0.83420 

0.80471 

0.020 

H(l’) 

1.15901 

1.07987 

1.03660 

0.026 

H(2') 

1.23032 

1.18590 

1.08056 

0.026 

H(4') 

0.64298 

1.21784 

1.00111 

0.021 

H(6A) 

0.41147 

1.14015 

0.87102 

0.027 

H(6'B) 

0.48205 

1.11487 

1.06167 

0.014 

H(7A) 

0.44868 

1.04194 

0.77938 

0.024 

H(7'B) 

0.33470 

1.03427 

0.91168 

0.021 

H(8’) 

0.66573 

1.01546 

1.12513 

0.012 

H(9') 

0.77340 

1.01922 

0.83423 

0.011 

H(llA) 

1.09205 

0.99035 

0.98873 

0.020 

H(ll’B) 

1.04447 

0.99134 

1.16493 

0.017 

H(12A) 

0.89389 

0.91066 

0.88280 

0.018 

H(12’B) 

1.05438 

0.89402 

1.07392 

0.029 

H(14’) 

0.56308 

0.93793 

0.84142 

0.012 

H(15A) 

0.28167 

0.92827 

0.90891 

0.023 

H(15A) 

0.40421 

0.93935 

1.10220 

0.016 

H(16') 

0.48007 

0.83935 

1.13620 

0.015 

H(17') 

0.62869 

0.84147 

0.87314 

0.014 

H(18A) 

0.86807 

0.94995 

1.27305 

0.017 

H(18'B) 

0.92207 

0.87636 

1.28493 

0.022 

H(18’C) 

0.70244 

0.89942 

1.27329 

0;032 

"Ueq=/2(SiSjUijai*a/ai.aj) 


Table  S2.  Anisotropic  thermal  parameters  of  non-Hydrogen  atoms  (A^)  and 
estimated  standard  deviations. 


Atom 

u„ 

U22 

U33 

U,2 

U,3 

U23 

Molecule  1 

0(1) 

0.0125(2) 

0.0088(2) 

0.0312(3) 

-0.0001(2) 

0.0108(2) 

0.0023(2) 

0(2) 

0.0146(2) 

0.0087(2) 

0.0217(3) 

0.0028(1) 

0.0102(2) 

0.0011(2) 

0(3) 

0.0113(2) 

0.0137(2) 

0.0187(2) 

-0.0034(2) 

0.0066(2) 

-0.0017(2) 

C(l) 

0.0083(2) 

0.0111(2) 

0.0211(2) 

0.0005(1) 

0.0039(2) 

0.0019(2) 

C(2) 

0.0094(2) 

0.0109(2) 

0.0250(3) 

-0.0003(1) 

0.0062(2) 

0.0016(2) 

C(3) 

0.0109(2) 

0.0089(1) 

0.0164(2) 

0.0002(1) 

0.0071(1) 

0.0008(1) 

C(4) 

0.0103(2) 

0.0096(1) 

0.0120(2) 

0.0005(1) 

0.0048(1) 

0.0009(1) 

C(5) 

0.0087(1) 

0.0086(1) 

0.0102(2) 

0.0005(1) 

0.0029(1) 

0.0007(1) 

C(6) 

0.0085(2) 

0.0104(2) 

0.0167(2) 

0.0010(1) 

0.0015(1) 

0.0032(1) 

C(7) 

0.0083(2) 

0.0104(2) 

0.0170(2) 

0.0018(1) 

0.0041(1) 

0.0024(1) 

C(8) 

0.0080(1) 

0.0089(1) 

0.0104(2) 

0.0008(1) 

0.0024(1) 

0.0010(1) 

C(9) 

0.0076(1) 

0.0087(1) 

0.0108(2) 

0.0011(1) 

0.0022(1) 

0.0008(1) 

C(10) 

0.0080(1) 

0.0090(1) 

0.0112(2) 

0.0004(1) 

0.0026(1) 

0.0007(1) 

C(ll) 

0.0085(2) 

0.0099(2) 

0.0194(2) 

0.0014(1) 

0.0054(1) 

0.0017(1) 

C(12) 

0.0086(2) 

0.0094(2) 

0.0163(2) 

0.0014(1) 

0.0034(1) 

0.0016(1) 

C(13) 

0.0092(1) 

0.0085(1) 

0.0106(2) 

0.0005(1) 

0.0042(1) 

-0.0005(1) 

C(14) 

0.0078(1) 

0.0086(1) 

0.0111(2) 

0.0010(1) 

0.0029(1) 

0.0007(1) 

C(15) 

0.0081(2) 

0.0115(2) 

0.0214(2) 

0.0006(1) 

0.0020(2) 

0.0023(2) 

C(16) 

0.0097(2) 

0.0096(2) 

0.0147(2) 

-0.0005(1) 

0.0041(1) 

-0.0009(1) 

C(17) 

0.0097(2) 

0.0080(1) 

0.0123(2) 

0.0010(1) 

0.0048(1) 

0.0001(1) 

qi8) 

0.0216(2) 

0.0145(2) 

0.0123(2) 

-0.0001(2) 

0.0089(2) 

-0.0018(2) 

Molecule  2 

0(1') 

0.0118(2) 

0.0071(1) 

0.0182(2) 

-0.0011(1) 

0.0073(2) 

-0.0007(1) 

0(2-) 

0.0171(2) 

0.0093(2) 

0.0201(3) 

0.0029(2) 

0.0096(2) 

0.0024(2) 

0(3') 

0.0145(2) 

0.0111(2) 

0.0167(2) 

-0.0051(2) 

0.0065(2) 

-0.0028(2) 

C(l') 

0.0082(1) 

0.0087(1) 

0.0217(2) 

0.0000(1) 

0.0066(1) 

-0.0009(2) 

C(2') 

0.0092(2) 

0.0083(2) 

0.0211(2) 

-0.0007(1) 

0.0067(2 

-0.0005(1) 

C(3’) 

0.0097(1) 

0.0073(1) 

0.0128(2) 

-0.0004(1) 

0.0053(1) 

-0.0003(1) 

C(4') 

0.0089(1) 

0.0069(1) 

0.0134(2) 

0.0002(1) 

0.0041(1) 

0.0002(1) 

C(5') 

0.0074(1) 

0.0070(1) 

0.0125(2) 

0.0004(1) 

0.0032(1) 

0.0001(1) 

C(6') 

0.0071(1) 

0.0089(1) 

0.0212(2) 

0.0007(1) 

0.0034(1) 

-0.0006(1) 

C(7’) 

0.0075(1) 

0.0087(2) 

0.0175(2) 

0.0003(1) 

0.0008(1) 

-0.0020(1) 

C(8') 

0.0072(1) 

0.0073(1) 

0.0115(2) 

-0.0001(1) 

0.0034(1) 

-0.0016(1) 

C(9’) 

0.0084(1) 

0.0067(1) 

0.0114(2) 

0.0003(1) 

0.0041(1) 

-0.0006(1) 

C(IO') 

0.0078(1) 

0.0076(1) 

0.0123(2) 

0.0004(1) 

0.0040(1) 

0.0000(1) 

C(ll') 

0.0080(1) 

0.0085(1) 

0.0175(2) 

0.0008(1) 

0.0047(1) 

0.0005(1) 

C(12') 

0.0103(2) 

0.0084(1) 

0.0157(2) 

0.0014(1) 

0.0069(1) 

0.0002(1) 

C(13') 

0.0090(1) 

0.0083(1) 

0.0089(1) 

0.0000(1) 

0.0040(1) 

-0.0007(1) 

C(14') 

0.0080(1) 

0.0073(1) 

0.0113(2) 

-0.0003(1) 

0.0032(1) 

-0.0020(1) 

C(15') 

0.0099(2) 

0.0091(2) 

0.0220(2) 

-0.0015(1) 

0.0078(2) 

-0.0035(2) 

C(16') 

0.0119(2) 

0.0080(1) 

0.0144(2) 

-0.0017(1) 

0.0066(1) 

-0.0019(1) 

C(17') 

0.0117(2) 

0.0076(1) 

0.0121(2) 

-0.0003(1) 

0.0060(1) 

-0.0013(1) 

C(18') 

0.0175(2) 

0.0149(2) 

0.0095(2) 

-0.0022(2) 

0.0044(2) 

-0.0009(1) 

The  thermal  parameters  are  listed  in  the  form  of  exp{-27i^(Uiih^a*^+2Ui2hka*b*+  — )} 
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Table  S3.  Bond  angles  (deg) 


bond 


m  1  m2 


C(2)-C(l)-H(l) 

C(10)-C(l)-H(l) 

C(l)-C(2)-H(2) 

C(3)-C(2)-H(2) 

C(3)-C(4)-H(4) 

C(5)-C(4)-H(4) 

C(5)-C(6)-H(6A) 

C(5)-C(6)-H(6B) 

C(7)-C(6)-H(6A) 

C(7)-C(6)-H(6B) 

H(6A)-C(6)-H(6B) 

C(6)-C(7)-H(7A) 

C(6)-C(7)-H(7B) 

C(8)-C(7)-H(7A) 

C(8)-C(7)-H(7B) 

H(7A)-C(7)-H(7B) 

C(7)-C(8)-H(8) 

C(9)-C(8)-H(8) 

C(14)-C(8)-H(8) 

C(8)-C(9)-H(9) 

C(10)-C(9)-H(9) 

C(ll)-C(9)-H(9) 

C(9)-C(11)-H(11A) 

C(9)-C(11)-H(11B) 

C(12)-C(11)-H(11A) 

C(12)-C(11)-H(11B) 

H(11A)-C(11)-H(11B) 

C(11)-C(12)-H(12A) 

C(11)-C(12)-H(12B) 

C(13)-C(12)-H(12A) 

C(13)-C(12)-H(12B) 

H(12A)-C(12)-H(12B) 

C(8)-C(14)-H(14) 

C(13)-C(14)-H(14) 

C(15)-C(14)-H(14) 

C(14)-C(15)-H(15A) 

C(14)-C(15)-H(15B) 

C(16)-C(15)-H(15A) 

C(16)-C(15)-H(15B) 

H(15A)-C(15)-H(15B) 

0(3)-C(16)-H(16) 


119.7(2) 

119.2(2) 

117.9(2) 

118.7(2) 

118.0(2) 

122.5(2) 

122.6(2) 

118.3(2) 

119.2(2) 

120.1(2) 

119.8(2) 

119.2(2) 

111.8(1) 

110.5(2) 

106.0(1) 

108.0(1) 

109.1(1) 

112.5(2) 

107.1(2) 

106.4(2) 

109.6(2) 

105.3(2) 

109.0(2) 

111.1(2) 

110.6(2) 

109.1(2) 

108.2(1) 

109.1(1) 

109.6(2) 

112.5(2) 

109.1(1) 

104.2(1) 

108.8(1) 

110.1(1) 

107.0(1) 

109.6(1) 

111.4(2) 

107.4(2) 

107.5(1) 

106.6(1) 

104.4(2) 

106.1(2) 

105.9(1) 

106.9(1) 

108.1(2) 

109.5(2) 

110.1(2) 

110.9(2) 

109.3(2) 

104.6(2) 

110.9(2) 

110.3(2) 

106.5(1) 

108.2(1) 

109.8(2) 

108.3(2) 

110.1(2) 

106.3(2) 

108.5(1) 

110.9(1) 

111.1(1) 

112.0(1) 

105.8(2) 

108.7(2) 

101.9(2) 

107.3(2) 

107.4(1) 

106.1(1) 

109.1(1) 

107.0(1) 

112.1(1) 

114.3(1) 

114.1(1) 

111.5(1) 

111.3(2) 

111.2(2) 

109.7(2) 

113.8(2) 

105.8(1) 

102.2(2) 

105.5(1) 

108.4(2) 
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C(15)-C(16)-H(16) 

C(17)-C(16)-H(16) 

0(2)-C(17)-H(17) 

C(13)-C(17)-H(17) 

C(16)-C(17)-H(17) 

C(13)-C(18)-H(18A) 

C(13)-C(18)-H(18B) 

C(13)-C(18)-H(18C) 

H(18A)-C(18)-H(18B) 

H(18A)-C(18)-H(18C) 

H(18B)-C(18)-H(18C) 


120.0(2) 

110.2(2) 

109.0(1) 

106.4(1) 

107.8(2) 

111.0(2) 

106.1(2) 

108.2(2) 

105.6(1) 

106.9(1) 

113.5(1) 

110.5(1) 

107.9(1) 

108.6(1) 

109.8(1) 

113.1(1) 

105.6(1) 

110.7(1) 

109.4(2) 

107.3(2) 

110.4(3) 

106.4(2) 
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Table  S4.  Multipole  population  coeffieients  for  the  atoms  of  molecule  1  (first  row) 
and  2  (second  row)“. 


All  the  ring  (A,  B,  C,  D)  atoms  of  the  molecule  1  has  been  constrained  equal  to  molecule  2. 
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The  Experimental  Charge  Density  Study  of  Estrogens:  17 /5-Estradiol  •Urea 
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In  order  to  relate  molecular  electrostatic  potential  to  biological  activities  of  estrogens,  a 
comparative  charge  density  study  of  different  derivatives  has  been  initiated.  The  first  completed 
charge  density  analysis  of  this  series  for  17>^estradiol«urea  is  presented  here.  It  has  been 
demonstrated  that  it  is  possible  to  perform  a  quality  experimental  charge  density  study  for  the  big 
organic  system  with  56  atoms  in  the  non-centrosymmetric  space  group.  The  special  tools  such  as 
the  optimal  coordinate  system  and  slow  initially  constrained  refinement  have  been  used  to 
accomplish  this  study. 

Our  results  for  the  urea  molecule  reasonably  agree  with  the  previous  experimental  and 
theoretical  results  showing  the  quality  of  our  data.  In  the  17yff-estradiol  molecule,  oxygen  atoms 
appeared  to  be  sp^  in  shape,  exhibiting  two  consistent,  distinct  lone  pairs  despite  of  different 
chemical  environments.  No  interaction  of  hydroxy  group  oxygen  lone  pairs  with  the  orbitals  of 
the  aromatic  ring  has  been  observed.  The  analysis  of  the  electrostatic  potential  revealed  that  the 
negative  potential  in  the  lone  pair  region  of  the  oxygen  atoms  is  quite  different.  The  Bader's 
topological  analysis  has  been  performed,  and  atomic  charges  have  been  estimated.  The  results 
are  compared  and  discussed. 


Estrogens  are  known  to  be  responsible  for  the  development  of  secondary  sexual  characteristics, 
as  well  as  they  effect  growth,  differentiation  and  function  of  wide  range  of  tissues  [1].  It  has  also 
been  shown  [2]  that  some  estrogens  are  responsible  for  the  initiation  and  progression  of  certain 
types  of  breast  cancer.  These  molecules  have  the  ability  to  bind  as  ligands  to  the  estrogen 
receptor  in  the  first  of  many  steps  which  could  result  in  the  activation  (agonistic  effect)  or 
repression  (antagonistic  effect)  of  genes  critical  in  the  mechanism  of  tumor  growth.  While  the 
mechanism  by  which  the  estrogens  influence  cancer  is  currently  unknown,  subtle  changes  in  the 
chemical  structure  of  estradiol  and  the  other  estrogens  are  known  to  elicit  different  biological 
responses  in  the  development  of  cancer.  It  has  been  suggested  [3]  that  the  agonistic/antagonistic 
responses  of  the  different  estrogens  can  be  related  to  such  physical  properties  as  their 
electrostatic  potential  and  charge  density  distribution. 

A  comparative  charge  density  study  of  a  series  of  estrogen  derivatives  has  been  initiated  in  our 
group,  with  the  primary  goal  of  relating  the  electronic  and  physical  structure  of  the  molecule  to 
biological  action.  In  addition  to  molecular  properties,  charge  density  studies  provide  significant 
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fundamental  information  about  the  functional  groups  and  the  atoms  which  make  them  up, 
helping  to  investigate  basic  chemical  principles  of  small  organic  molecules  in  general.  This  paper 
will  pay  particular  attention  to  the  hydroxy  groups  of  the  estrogen  molecule.  The  hydroxy  groups 
are  some  of  the  most  chemically  interesting  portion  of  the  estrogen  molecule,  and  their 
interaction  with  the  receptor  is  thought  to  be  mostly  responsible  for  their  activity. 

Described  here  is  the  charge  density  study  of  17y0-estradiol»urea.  This  particular  crystal  system 
was  chosen  to  be  the  first  in  the  series  for  two  reasons.  First,  ITy^estradiol  (Ej)  is  the  most 
common  naturally  occurring  estrogen  and  is  typically  the  standard  by  which  all  other  estrogenic 
activities  are  related.  Ea  demonstrates  an  agonistic  response  in  in  vitro  tests  which  typically 
utilize  the  estrogen  receptor  in  either  competitive  yeast  or  human  cell  assays  [4-5].  Secondly,  it 
contained  a  solvated  urea  molecule  which  has  been  extensively  studied.  Both  theoretical  and 
experimental  charge  density  studies  have  been  performed  on  this  molecule  [6-7].  This  provides  a 
direct  means  of  comparison  to  test  our  methods  against  previous  results. 

Experimental 

Data  Collection  and  Reduction 

Crystals  were  grown  by  slow  evaporation  as  described  by  Duax  in  the  original  publication  of 
the  crystal  structure  [8].  The  clear,  colorless  crystal  of  approximate  dimensions  0.35  x  0.37  x 
0.40  mm  was  attached  to  a  30  /an  carbon  fiber  using  a  small  amount  of  epoxy  resin.  Data  were 
collected  at  100.0(1)  K  on  a  Bruker  Platform  Diffractometer  equipped  with  a  SMART  6000 
CCD  area  detector  located  at  7.12  cm  from  the  goniometer  center.  The  Oxford  Cryostream 
system  has  been  used  to  cool  the  crystal.  The  crystals  were  irradiated  with  graphite 
monochromated  Mo  Ka  radiation  (k  =  0.7107  A).  0.3“  omega  scans  were  carried  out  at  four  phi 
settings  and  three  detector  positions  in  20,  resulting  in  more  than  7400  frames.  A  frame  time  of 
60  seconds  was  used  for  the  low  angle  setting,  and  180  seconds  for  the  medium  and  high  angle 
settings.  A  complete  experimental  protocol  is  deposited. 

The  intensities  were  integrated  using  the  program  SAINT  [9].  After  several  integrations  and 
preliminary  refinements,  it  was  found  that  different  integration  parameters  were  required  for 
each  detector  setting  to  obtain  the  best  possible  structure  factors.  The  box  sizes,  profile  fitting 
parameters,  and  the  simple  sum  perimeter  limits  were  adjusted  for  each  detector  setting.  The  unit 
cell  parameters  have  been  finally  refined  in  SAINT. 

Lorentz-polarization  and  A/2  corrections  were  applied  to  the  integrated  intensities  with  a  A/2 
correction  factor  of  0.0012,  this  value  being  taken  from  the  earlier  study  [10].  The  program 
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SORTAV  [11]  was  used  for  outlier  determination,  averaging,  and  applying  the  absorption 
correction.  Only  reflections  with  I  >  3  0(7)  and  below  1.00  A  *  were  used  during  the  multipole 
refinements. 

Table  1.  The  crystal  information  and  experimental  details. 


Formula 

C19H28N2O3 

Formula  Weight 

332.43 

Crystal  System 

Orthorhombic 

Space  Group 

P2i2i2i 

a 

7.9022(9)  A 

b 

9.2228(10)  A 

c 

24.5890(28)  A 

II 

II 

90.0" 

Volume 

1792.06  A^ 

Z 

4 

T 

100.0  K 

/I 

0.7107  A 

(sin  ^^)max 

1.180  A* 

Reflections  Collected 

110999 

Rejected  Outliers 

779 

Unique  Reflections 

13187 

Included  in  the  Refinement  [/>  3o(7)] 

7626 

Data  completeness 

98.6  % 

Average  redundancy 

8.4 

Rint 

0.0411 

R]  (spherical) 

0.065 

wR2  (spherical) 

0.145 

R(F2)  [7  >30(7)1  (multipole) 

0.029 

R(F^)  all  data  (multipole) 

0.036 

wR(F^) 

0.045 

Weighting  scheme 

1/a^ 

Goodness  of  Fit 

1.393 

n\ 


Least-Squares  ReHneinents 

The  crystal  structure  was  re-solved,  and  the  spherical  atom  model  was  refined  using  the 
program  SHELXL-97  [12].  17P-estradiol  and  urea  molecules  are  shown  in  Figure  1.  All 
hydrogen  atoms  were  located  from  the  difference  Fourier  map.  The  results  of  this  refinement 
served  as  the  starting  point  for  the  aspherical  atom  refinements  that  were  performed  using  the  XD 
software  package  [13]. 


CHS) 


0(31 


Figure  1.  17P-estradiol  and  urea  molecules  showing  50%  probability  thermal  ellipsoids  as  at 
lOOK.  Hydrogen  atoms  are  omitted  for  clarity. 


The  aspherical  atom  model  [14],  used  in  further  refinements  consists  of  three  components.  The 
first  and  second  describe  spherically  averaged  core  and  valence  densities  normalized  to  one 
electron,  while  the  third  one  describes  the  aspherical  valence  density.  p(r)  is  the  electron  density. 
Pc,  Pv,  and  Pin,  are  population  coefficients,  the  core  population  (Pc)  was  fixed  in  this  study. 
Asphericity  is  taken  into  account  by  expanding  of  the  electron  density  over  real  combinations  of 
spherical  harmonics,  yim,  which  satisfy  local  symmetry.  The  variable  radial  parameters,  /rand  /r', 
allow  for  expansion  or  contraction  of  the  valence  terms.  Ri  are  normalized  Slater-type  radial 
functions. 

Patomic(r)  =  PcPcore  +  PvK  ^  P  valence  (KT) 

4  1 

+  ]^K'^Ri(K'r)  ^PimyimCr/r) 

1=1  m=“l 
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Because  the  model  is  based  on  atom  centered  multipoles,  a  coordinate  system  is  required 
for  each  atom  in  the  system.  The  core  structure  of  the  estrogen  molecule  was  set  up  using  the 
standard  estrogen  coordinate  system  to  take  advantage  of  the  shape  of  certain  multipoles  to 
increase  the  stability  of  the  initial  multipole  refinements  [15].  The  vectors  for  the  oxygen  atoms 
of  the  estrogen  molecule  were  directed  to  bound  atoms.  The  coordinate  system  setup  for  the  urea 
molecule  is  deposited. 

Before  the  refinements  were  begun,  the  hydrogen  atom  positions  were  corrected  by  extending 
the  X  -  H  distances  to  their  neutron  bond  lengths  (Csp2H  =  1.08  A,  CspsH  =  1.10  A,  Csp3H2  =  1.09 
A,  CspsHa  =  1.06  A,  OH  =  0.97  A,  Nsp2H2  =  1.01  A).  These  values  were  fixed  in  the  refinement. 
A  well  known  correlation  exists  between  the  thermal  parameter  and  the  kappa  parameters  of  a 
hydrogen  atom,  only  allowing  the  kappa  parameters  to  be  refined  when  the  thermal  parameters 
of  the  hydrogen  atom  are  known  [16].  To  avoid  this  problem,  fixed  kappa  values  (/t//  =  1.20  and 
k'h  =  1.20)  were  used,  while  the  thermal  parameters  were  allowed  to  refine  [17-18]. 

Initially,  to  reduce  the  number  of  parameters  in  the  least  squares  process,  chemical  constraints 
and  only  the  most  significant  multipoles  were  considered.  By  the  end  of  the  refinement,  all  non¬ 
hydrogen  atoms  were  treated  up  to  the  octupole  level  with  no  chemical  constraints.  The 
hydrogen  atoms  were  refined  with  only  two  multipoles  having  cylindrical  symmetry  along  the  X 
-  H  bond.  Hydrogen  atoms  bound  to  the  same  carbon  have  been  constrained  to  be  equal.  The 
electroneutrality  constraint  has  been  imposed  on  both  of  two  molecules.  Complete  tables  giving 
the  final  structural  information,  such  as  atomic  coordinates,  Uy’s,  Uiso’s,  bond  lengths,  bond 
angles  and  multipole  parameters  are  deposited. 

Topological  Analysis,  Atomic  Charges  and  Electrostatic  Potential 

The  topological  analysis  was  performed  based  on  the  Atoms  in  Molecules  (AIM)  Quantum 
Theory  [19].  The  XD  property  assessment  program,  XDPROP,  was  used  to  perform  a 
quantitative  critical  point*  analysis  on  all  covalent  and  hydrogen  bonding  interactions,  as  well  as 
oxygen  lone  pair  densities.  The  atomic  basin  boundaries  were  defined  and  the  atomic  charges 
(Bader's  charges)  were  integrated  within  the  basins  with  the  program  TOPXD  [20].  The 
integrated  atomic  charge  q(ii)  over  the  atomic  basin  is  defined  as  a  difference  between  nuclear 
Zq  and  electronic  N{Q.)  charges: 


*  Critical  points  are  special  points  in  the  electron  density,  where  Vp(rcp)— 0.  The  saddle  (3,-1)  point  reflects  the 
chemical  bonding.  In  the  Laplacian,  the  (3,+3)  minimum  point  corresponds  to  the  maximum  curvature  of  the 
electron  density  in  the  charge  concentration  regions.  See  ref.  [191. 
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q(Q)  =  7^ -N{Q.X  N(a)=^j  p(r)dT. 

£2 

This  method  allows  one  an  exhaustive  partitioning  of  the  molecules  into  the  respective  atoms, 
leading  to  the  additivity  and  transferability  of  it's  properties  [21]. 

The  calculation  of  the  electrostatic  potential  (ESP)  was  performed  using  the  program 
XDPROP.  It  can  be  calculated  by  the  evaluation  of 


rP(r’)dr’ 
-r| 


where  Za  is  the  charge  on  nucleus  A  located  at  Ra.  The  construction  assembles  positive  charges 
at  the  nuclei  and  a  continuous  distribution  of  electrons.  The  calculations  were  performed  based 
on  the  isolated  from  the  crystal  P-estradiol  molecule. 


Results  and  Discussion 

Packing 

The  molecular  structure  and  hydrogen  bonding  pattern  were  found  to  be  the  same  as  described 
by  Duax  [8].  The  packing  of  this  structure  results  in  a  formation  of  six  hydrogen  bonds  ranging 
in  donor-acceptor  length  from  2.64  -  3.10  A,  these  are  listed  in  Table  2.  The  0(1)  atom  acts  as 
an  acceptor  in  two  hydrogen  bonds  with  hydrogens  of  two  independent  urea  molecules  (H(INB) 
and  H(2NB)).  The  H  -0  vectors  of  these  hydrogen  bonds  are  pointed  almost  directly  at  the 
positions  of  the  lone  pairs  of  0(1).  The  0(2)  atom  has  a  very  similar  hydrogen  bonding  structure 
to  0(1)  in  which  it  acts  as  an  acceptor  in  two  hydrogen  bonds  and  a  donator  in  one.  The 
hydrogen  atoms  from  the  0(1)  atom  of  another  independent  E2  molecule  and  H(2NB)  of  a  urea 
molecule  line  up  almost  directly  with  the  lone  pair  positions  of  0(2).  Topological  analysis 
yielded  in  a  (3,-1)  critical  point  and  a  bond  path^  for  each  hydrogen  -  acceptor  interaction, 
verifying  the  true  bonding  interaction  in  these  hydrogen  bonds  [22].  The  complete  table  of 
results  from  hydrogen  bond  topology  is  deposited. 


^  Bond  path  is  a  line  along  which  the  electron  density  decreases  for  any  lateral  displacement. 
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Table  2.  The  geometrical  information  for  the  hydrogen  bonds. 


D-H-A 

d  (D-A),  A 

ZDHA,® 

0(2)-H(20)-  -0(3) 

2.6395(7) 

169.1(16) 

0(1)-H(10)*  -0(2) 

2.6698(7) 

171.0(20) 

N(1)-H(1NA)-  -0(3) 

2.8231(10) 

168.6(18) 

N(1>-H(1NB)-  *0(1) 

2.9783(9) 

135.2(16) 

N(2)-H(2NB)-  -0(1) 

3.0086(9) 

161.1(17) 

N(2)-H(2NA)-  -0(2) 

3.1023(8) 

169.1(16) 

Kappa  Parameters 

The  coefficients  of  expansion/contraction  of  the  spherical  valence  and  deformation 
densities,  K  and  k’  parameters  respectively,  are  shown  in  Table  3.  The  grouping  of  individual 
atoms  were  based  on  the  atom  type,  hybridization  and  chemical  environment.  A  total  of  twelve 
kappa  sets  were  utilized  in  order  to  allow  the  necessary  flexibility  while  attempting  to  maintain  a 
minimum  number  of  parameters.  As  mentioned  earlier,  the  hydrogen  atom  kappa  values  were 
fixed  at  1.20.  All  other  kappa  values  were  allowed  to  refine  during  the  least  squares 
minimization. 

Table  3.  The  expansion/contraction  (k  and  K')parameters  for  the  17P-estradiol*urea  crystal 


Atoms 

Kappa 

K 

K’ 

0(1),  0(2) 

1 

0.971 

0.854 

C(3) 

2 

1.009 

0.879 

C(17) 

3 

1.012 

0.933 

C(1),C(2),C(4) 

4 

0.976 

0.896 

C(5),C(10) 

5 

0.986 

0.838 

C(6),C(7),C(8),C(9), 

C(ll),  C(12),  C(13),  C(14), 

6 

0.985 

0.878 

C(15),  C(16),C(17),C(18) 

all  C  -  H  hydrogen  atoms 

7 

1.200 

1.200 

8 


H(10),  H(20) 

8 

1.200 

1.200 

C(19) 

9 

0.986 

0.818 

0(3) 

10 

0.991 

1.080 

N(1),N(2) 

11 

0.996 

1.004 

all  N  -  H  hydrogen  atoms 

12 

1.200 

1.200 

Atomic  Charges 

The  net  atomic  charges  were  derived  directly  from  the  monopole  populations  of  the  multipole 
model,  as  well  as  from  integrating  of  the  electron  density  within  the  atomic  basins  as  described 
by  the  AIM  theory.  A  complete  list  of  atomic  charges  determined  by  both  methods  are  listed  in 
Table  4. 

As  mentioned  earlier,  the  urea  molecule  gives  us  the  opportunity  to  compare  our  results  with 
previous  theoretical  and  experimental  results  [6-7].  The  atomic  charges  derived  directly  from  the 
monopole  populations  of  the  multipole  model  compared  favorably  to  the  experimental  values 
reported*  by  Zavodnik  et  al.  [6],  with  values  matching  within  one  e.s.d.  Bader’s  charges  are  in  a 
fair  agreement  with  those  reported  by  Gatti  et  al.  [7].  Generally,  Bader's  charges  are  significantly 
bigger  than  the  charges  derived  from  the  monopole  populations;  this  is  not  surprising  due  to  the 
different  ways  of  calculations.  Upon  first  inspection,  a  charge  of  +2.13  e  for  the  central  carbon 
seems  tremendously  high.  At  the  same  time,  there  is  evidence  in  the  literature,  from  both 
theoretical  and  experimental  data,  that  this  is  not  unreasonable  for  a  carbon  surrounded  by 
electronegative  oxygen  and  nitrogen  atoms  [7, 23]. 

The  net  atomic  charges  for  the  atoms  of  the  estrogen  molecule  derived  from  the  monopole 
populations  of  the  multipole  model  were  as  expected  for  the  given  atom  types  and  their  chemical 
neighbors.  In  addition,  like  atoms  with  similar  chemical  environments  exhibited  very  consistent 
monopole  populations.  Again,  the  charge  integrations  over  atomic  basins  yielded  much  greater 
charges  for  the  oxygen  atoms  and  the  atoms  directly  bound  to  them,  maintaining  the  consistency 
between  chemically  equivalent  atoms. 

Table  4.  Net  atomic  charges  for  the  atoms  of  17P-estradiol  •  urea 


*  See  the  deposition  material  of  Zavodnik  et  al.  (1999)  Acta  Crystallogr.  B55, 45-54. 
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Atom 

q(Pv) 

q(^) 

Atom 

q(Pv) 

q(i2) 

0(1) 

-0.534(9) 

-1.24 

H(10) 

0.401(9) 

0.68 

0(2) 

-0.527(9) 

-1.11 

H(20) 

0.391(9) 

0.67 

C(l) 

-0.233(16) 

-0.19 

H(l) 

0.206(10) 

0.14 

C(2) 

-0.227(16) 

-0.12 

H(2) 

0.248(9) 

0.20 

C(3) 

0.148(14) 

0.49 

H(4) 

0.249(9) 

0.21 

C(4) 

-0.267(15) 

-0.21 

H(6A) 

0.147(7) 

0.15 

C(5) 

-0.110(15) 

-0.09 

H(6B) 

0.147(7) 

0.16 

C(6) 

-0.265(16) 

-0.28 

H(7A) 

0.146(7) 

0.08 

C(7) 

-0.279(16) 

-0.15 

H(7B) 

0.146(7) 

0.10 

C(8) 

-0.121(14) 

-0.12 

H(8) 

0.180(9) 

0.16 

C(9) 

-0.131(15) 

-0.07 

H(9) 

0.180(9) 

0.09 

C(10) 

-0.106(15) 

-0.05 

H(llA) 

0.147(7) 

0.08 

C(ll) 

-0.277(15) 

-0.16 

H(llB) 

0.147(7) 

0.09 

C(12) 

-0.273(15) 

-0.14 

H(12A) 

0.141(7) 

0.07 

C(13) 

-0.158(16) 

-0.08 

H(12B) 

0.141(7) 

0.06 

C(14) 

-0.111(14) 

-0.04 

H(14) 

0.176(8) 

0.11 

C(15) 

-0.294(15) 

-0.23 

H(15A) 

0.141(7) 

0.10 

C(16) 

-0.273(16) 

-0.18 

H(15B) 

0.141(7) 

0.13 

C(17) 

0.168(13) 

0.36 

H(16A) 

0.146(7) 

0.11 

C(18) 

-0.402(17) 

-0.24 

H(16B) 

0.146(7) 

0.14 

0(3) 

-0.234(9) 

-1.27 

H(17) 

0.120(9) 

0.14 

N(l) 

-0.282(12) 

-1.41 

H(18A) 

0.122(6) 

0.06 

N(2) 

-0.286(11) 

-1.40 

H(18B) 

0.122(6) 

0.06 

C(19) 

-0.030(14) 

2.13 

H(18C) 

0.122(6) 

0.06 

H(INA) 

0.205(10) 

0.49 

H(INB) 

0.209(11) 

0.49 

H(2NA) 

0.207(10) 

0.47 

H(2NB) 

0.208(10) 

0.49 

q(Pv)  -  net  atomic  charges  derived  directly  from  the  monopole  populations  of  the  multipole  model; 
q(Q)  -  Bader’s  net  atomic  charges 
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a)  urea  molecule 

Overall  values  of  the  electron  density  at  the  bond  critical  points  are  in  a  fair  agreement  with 
those  reported  by  Zavodnik  et  al.  [6]  and  Gatti  et  al.  [7],  with  the  exception  of  the  N  -  H  bonds 
(Table  5).  Zavodnik  et  al.  [6]  reported  significant  differences  (-0.5  eA'^)  between  the  theoretical 
and  experimental  electron  densities  at  the  (3,-1)  critical  point  for  the  N  -  H  bond  of  the  urea 
molecule.  The  p(x)  values  in  this  study  are  much  closer  to  the  theoretical  values.  The  difference 
in  the  Laplacian  values  at  (3,-1)  critical  points  between  experimental  and  theoretical  values  is 
much  bigger,  showing  the  sensitivity  of  the  electron  density  curvature  to  the  different  approaches 
of  it's  estimation,  including  the  deficiency  of  the  multipole  model  [24]. 


Table  5.  Topological  analysis  for  the  urea  molecule. 


this  work 

Zavodnik,  et  al.  [6] 

Gatti,  et  al.  [7] 

experiment 

experiment 

theory 

theory 

bond 

Ar) 

fKr) 

V"Ar) 

Ar) 

P(X) 

C-0 

2.79 

-33.15 

2.54 

-18.86 

2.57 

-7.92 

2.57 

-7.95 

(C“N)average 

2.26 

-24.20 

2.54 

-37.66 

2.35 

-27.61 

2.36 

-27.71 

(N-H)average 

2.22 

-28.91 

1.82 

-31.91 

2.34 

-47.06 

2.34 

-47.23 

The  units  for  the  electron  density,  /0(r),  are  eA‘*,  and  for  the  Laplacian,  V^/0(r),  they  are  eA‘ 


A  critical  point  search  was  also  performed  in  the  Laplacian  to  find  the  charge  density 
concentration  associated  with  the  lone  pairs  on  the  urea  oxygen.  The  critical  points  are  (3, +3)  in 
type  indicating  that  V^/0(r)  is  a  local  minimum  at  rc.  The  critical  point  information,  as  well  as  a 
plot  of  the  Laplacian  in  the  plane  of  the  lone  pairs  of  the  0(3)  atom  can  be  seen  Table  6  and  Fig. 
2  respectively.  The  angle  between  the  lone  pairs,  being  greater  that  140®,  is  much  larger  than  the 
ideal  trigonal  planar  geometry  in  agreement  with  VESPR  theory.  The  inconsistency  between  the 
C(3)  -  0(1)  -  LP  angles  are  most  likely  due  to  the  shallow  curved  shape  of  the  Laplacian  in  the 
lone  pair  region.  An  extremely  small  change  in  V^/7(r)  could  easily  move  the  critical  point  by  5  - 
10°.  Gatti  et  al.  [7]  also  found  two  (3,+3)  critical  points  in  the  Laplacian,  corresponding  to  the 
urea  oxygen  lone  pairs.  Although  one  of  their  points  is  very  similar  to  those  in  Table  6,  having 
the  VVc(r)  s  -137  eA'^  and  the  0(1)-LP  distance  of  -0.339  A,  another  (3,+3)  point  is 
significantly  away  from  the  oxygen  nucleus  (-0.383  A)  having  much  lower  value  of  the 
Laplacian  (~  -48  eA'^).  As  it  has  been  shown  before,  the  Laplacian  values  (and  it's  curvatures) 
are  very  'method-sensitive'. 
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Table  6.  Lone  pair  Laplacian  critical  point  information  for  the  0(3)  atom. 


Atom 

Lone 

Pair 

vVW 

(eA-^) 

CP 

Type 

0(1)-LP 

(A) 

LP1-0(1)-LP2 

n 

C3-0(1)-LP 

(°) 

0(3) 

LPl 

-149 

(3, +3) 

0.337 

144 

118 

LP2 

-148 

(3, +3) 

0.336 

144 

92 

Figure  2.  The  negative  part  of  the  Laplacian  in  the  plane  of  the  0(3)  atom  lone  pairs.  Contour 
interval  is  -5  eA'^  starting  at  -100  eA'^. 

b)  E2  molecule 

The  critical  point  analysis  of  covalent  bonds  of  the  estrogen  molecule  demonstrated  the 
characteristic  (3,-1)  critical  points,  with  expected  yO(rc),  V^/j(rc),  and  ellipticity  values  given  the 
hybridization  and  connectivity.  The  complete  list  of  all  bond  critical  point  properties  is 
deposited.  The  average  /X^c)  values  are  1.91  eA'^  for  the  C  -  O  bond,  2.17  eA  ^  for  Car  -  Car,  1.66 
eA'^  for  C  -  C,  1.91  eA'^  for  C  -  H  and  2.07  eA'^  for  the  O  -  H  bond.  The  critical  points  are  well 
centered  in  the  homonuclear  C  —  C  bonds,  and  proportionally  displaced  away  from  the  more 
electronegative  atom  in  the  heteronuclear  bonds. 

Critical  points  in  the  Laplacian  corresponding  to  the  lone  pairs  of  both  hydroxy  oxygen  atoms 
are  of  the  special  interest.  Table  7  indicates  the  positions  of  the  lone  pair  critical  points  in  the 
Laplacian.  The  deformation  electron  density  map  of  the  lone  pairs  of  the  oxygen  0(1),  bound  to 
the  aromatic  ring  is  presented  in  Fig.  3,  while  plot  of  the  Laplacian  in  the  plane  of  the  lone  pairs 
of  the  0(2)  atom  is  shown  in  Fig.  4.  The  finding  of  two  distinct  lone  pair  regions  after  the 
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multipole  refinement  for  0(1)  was  somewhat  surprising.  We  expected  that  the  aromatic  hydroxy 
group,  in  which  hydrogen  lies  in  the  plane  of  the  aromatic  ring,  would  be  of  sp  type  geometry 
having  some  interaction  with  the  7t-  orbitals  of  the  aromatic  ring  (the  same  as  in  phenol), 
however,  the  results  of  this  study  indicate  that  this  is  not  the  ground  state  configuration  in  the 
estrogen  molecule.  The  distances  of  the  lone  pairs  from  the  oxygen  atoms  are  very  consistent. 
This  is  believed  to  be  at  least  partially  to  due  to  the  fact  that  the  contraction/expansion 
coefficients  for  both  oxygen  atoms  are  constrained  to  be  equal.  Same  as  in  the  urea  molecule, 
the  angles  between  the  lone  pairs  agree  with  VESPR  theory  in  that  they  are  much  larger  than  the 
ideal  tetrahedral  value. 

Table  7.  Lone  pair  Laplacian  critical  point  information  for  the  0(1)  and  0(2)  atoms. 


Atom 

Lone 

Pair 

vVrc) 

(eA-=) 

0-LP 

(A) 

LP1-0-LP2 

(°) 

C-O-LP 

(°) 

H-O-LP 

(°) 

0(1) 

LP3 

-133 

0.346 

146 

104 

103 

LP4 

-119 

0.347 

146 

90 

100 

0(2) 

LP5 

-129 

0.346 

235 

108 

115 

LP6 

-120 

0.346 

235 

88 

101 

Figure  3.  Dynamic  model  deformation  density  map  in  the  plane  of  the  lone  pairs  of  the  0(1) 
atom.  Contour  intervals  are  0.05  eA'^  with  solid  lines  positive,  dashed  lines  zero,  and  dotted  lines 
negative. 
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Figure  4.  The  negative  part  of  the  Laplacian  in  the  plane  of  the  lone  pairs  of  the  0(2)  atom. 
Contour  intervals  are  -5  eA'^  starting  at  -80  eA'^. 


Electrostatic  Potential 

The  core  structure  of  the  estrogen  molecule  generally  is  hydrophobic  and  relatively  predictable 
in  terms  of  the  electrostatic  potential.  The  ESP  plotted  on  the  van  der  Waals  molecular  surface 
for  the  E2  molecule  with  the  program  package  MOLMOL  [25]  is  shown  in  Fig.  6.  Since  the 
isolated  molecule  (taken  from  the  crystal)  is  calculated  in  this  case,  the  ESP  does  not  reflect  the 
interactions  in  the  condensed  phase,  such  as  hydrogen  bonding.  The  negative  potential  areas 
(red)  are  produced  by  the  aromatic  ring,  oxygens  and  some  of  carbon  atoms  (Fig.  6),  while  the 
hydrogen  atoms  involved  in  the  hydrogen  bonding  produce  highly  positive  areas  (blue)  on  the 
molecular  surface. 


Figure  6.  The  E2  electrostatic  potential  plotted  on  the  van  der  Waals  molecular  surface.  Negative 
values  are  in  red,  positive  are  in  blue. 

The  hydroxy  groups  are  of  particular  interest  in  the  E2  molecule,  offering  varying 
environments  being  bound  to  sp^  and  sp^  type  carbons  and  being  involved  in  several  hydrogen 
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bonds.  The  consistency  in  the  lone  pair  density,  in  shape,  size,  position,  and  hydrogen  bonding 
geometry  suggested  that  the  electrostatic  potential  surrounding  the  oxygen  atoms  might  be 
relatively  consistent  in  shape  and  magnitude.  This  however  was  not  the  case.  Fig.  7  demonstrates 
the  electrostatic  potential  of  the  0(1)  hydroxy  group  lone  pair  region.  The  significant  portion  of 
the  negative  ESP  extends  approximately  180°  around  the  oxygen.  A  maximum  ESP  of  ~  -0.20 
eA  *  was  observed  for  this  hydroxy  group.  The  additional  negative  regions  displayed  represent 
the  potential  created  by  the  n  density  of  the  aromatic  ring.  The  oxygen  atom  forms  a  continuum 
with  the  negative  potential  above  and  below  the  aromatic  ring. 


Figure  7.  Electrostatic  potential  isosurface  around  the  0(1)  atom.  The  blue  surface  represents 
+1.0  eA  ’,  while  the  red  surfaces  represent  -0.15  eA  *. 

Before  discussing  the  0(2)  hydroxy  group,  it  is  important  to  recognize  the  orientation  of  the 
hydrogen  atom  H(20).  A  hydroxy  group  attached  to  an  aromatic  ring  almost  exclusively  has  the 
hydrogen  atom  in  the  plane  of  the  ring.  Obviously,  a  hydroxy  group  bound  to  an  sp  carbon  can 
freely  rotate.  In  the  solid  state,  the  orientation  of  the  OH  group  will  almost  always  be  defined  by 
the  hydrogen  bonding  scheme.  Fig.  8  displays  the  D-ring^  of  the  E2  molecule  and  the  torsion 
angle  (y)  of  49°  between  H(20)  and  H(17). 


5  The  C(13)-C(14)-C(15)-C(16)-C(17)  ring. 
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Figure  8.  The  0(2)  hydroxy  group  geometry  showing  torsion  angle  (y)  between  the  H(20)  and 
H(17)  atoms.  The  0(2)  atom  is  shown  in  red. 

The  negative  electrostatic  potential  around  0(2)  has  much  greater  maximum  value  than  that  of 
the  0(1)  hydroxy  group,  reaching  a  value  of  ~  -0.35  eA’’.  In  addition,  the  ESP  does  not  evenly 
cover  the  lone  pair  region  as  on  the  0(1)  hydroxy  group.  A  shift  toward  the  lone  pair  opposite 
H(17)  appears  to  have  taken  place  (Fig.  9).  The  coverage  area  of  the  negative  potential  seems 
smaller,  covering  near  120°,  versus  the  approximate  180°  coverage  on  0(1). 


Figure  9.  Electrostatic  potential  isosurface  around  the  0(2)  atom.  The  blue  surface  represents 
+1.0  eA  ’,  while  the  red  surfaces  represent  -0.25  eA  '. 

It  is  believed  that  the  agonistic  or  the  antagonistic  responses  of  different  estrogens  can  be 
related  to  their  electrostatic  potential  [3].  We  intend  to  construct  a  database  with  ESP  calculated 
for  estrogens  with  different  biological  activities.  Once  enough  data  are  available,  we  hope  to 
identify  specific  topological  features  in  these  estrogens  that  distinguish  binding  ability  of  these 
molecules  to  the  estrogen  receptor. 
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Conclusions 

This  is  one  of  the  first  charge  density  studies  accomplished  for  the  big  (56  atoms)  organic 
molecules  in  the  non-centrosymmetric  unit  cell.  Our  study  showed  that  it  is  possible  to  get  a 
quality  data  for  such  kind  of  system;  the  optimal  coordinate  system  and  slow,  initially 
constrained  refinement  are  the  necessary  tools  for  this  type  of  study.  Since  most  estrogens  are 
identical  in  the  core  structure  of  the  molecule,  this  study  will  provide  a  better  starting  point  for 
future  multipole  refinements  of  other  estrogen  derivatives.  Overall  our  results  pertaining  to  the 
urea  molecule  reasonably  agree  with  the  previous  literature  results  showing  the  quality  of  our 
data. 

For  the  E2  molecule,  the  lone  pair  densities  of  the  oxygen  atoms,  which  are  an  extremely  small 
features  of  the  system,  were  found  to  be  quite  similar.  This  indicates  that  the  oxygen  atom  bound 
at  the  C(3)  position  do  not  have  a  strong  interaction  with  the  ;rorbitals  of  the  aromatic  ring  in  the 
ground  state  configuration.  Despite  this  consistency,  the  electrostatic  potentials  in  the  region  of 
the  hydroxy  groups  differ  greatly  in  both  shape  and  magnitude.  Future  studies  may  reveal  how 
the  H(17)  -  C(17)  -  0(2)  -  H(20)  torsion  angle  affects  the  electrostatic  potential  of  the  C(17) 
hydroxy  group.  Obviously  many  more  estrogen  systems  must  be  completed  before  the 
comparative  portion  of  the  study  can  take  place  and  real  conclusions  can  be  made  on  the 
relationship  between  the  electrostatic  potential  and  estrogenic  action. 
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and  Prof  P.  Coppens  for  making  TOPXD  program  available  for  us.  We  would  also  like  to  thank 
the  Department  of  Defense  for  their  financial  support  (Grant  #  DAMD 17-00- 1-0468  and  # 
DAMD17-99-1-9408). 
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Fig.  Dl.  Local  multipole  coordinate  system  for  the  urea  molecule 
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Table  S2.  Fractional  atomic  coordinates  for  E2*urea  after  final  refinement. 
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Table  S8.  Multipole  populations  of  Carbon  atoms  of  E2*urea.  Populations  less  than  one  estimated  standard  deviation  were  set  to 
zero. 
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Table  S9.  Multipole  populations  of  Hydrogen  atoms  of  E2*urea.  Populations  less  than  one  estimated  standard  deviation  were  set 
to  zero. 
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Table  S 10.  Topological  properties  of  bond  critical  points  in  E2*urea. 
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Table  S 1 1 .  Topological  properties  of  bond  critical  points  in  E2«urea  continued. 
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Appendix  H. 


Crystal  structure  and  charge  density  analysis  of  17P-estradiol*V2  MeOH 
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C(180 


Figure  1.  Thermal  ellipsoid  plot  of  17;5-estradiol«*/2methanol  where  ellipsoids  represent 
50%  probability  electron  density  of  the  atom.  Hydrogen  atoms  are  omitted  for 
clarity. 


Run 

le 

(O 

Scan 
Width  (°) 

#of 

Frames 

Frame 
Times  (sec) 

1 

-40 

-46 

22 

-0.15 

1100 

16 

2 

-40 

-46 

112 

-0.15 

1100 

16 

3 

-40 

-46 

202 

-0.15 

1100 

16 

4 

-40 

-46 

292 

-0.15 

1100 

16 

5 

-80 

-86 

67 

-0.15 

1100 

32 

6 

-80 

-86 

157 

-0.15 

1100 

32 

7 

-80 

-86 

247 

-0.15 

1100 

32 

8 

-80 

-86 

337 

-0.15 

1100 

32 

9 

-102 

-108 

22 

-0.15 

1100 

64 

10 

-102 

-108 

112 

-0.15 

1100 

64 

11 

-102 

-108 

202 

-0.15 

1100 

64 

12 

-102 

-108 

292 

-0.15 

1100 

64 

Table  1.  Data  collection  parameters  for  17>0-estradiol»‘/2methanol. 
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Crystal  Data 

Chemical  Formula 

C37H52O5 

Temperature 

100.0(1)  K 

Crystal  Dimensions 

0.22  X  0.26  X  0.42  mm 

Space  Group 

PI 

A 

7.2910(1)  A 

B 

9.2768(1)  A 

C 

12.3873(2)  A 

a 

89.4704(6) 

e 

87.8577(6) 

Y 

70.7607(7) 

Volume 

790.489(33) 

Z  (Crystallographic) 

2 

Integration  Par 

ameters 

Box  Size  (°) 

Profile 
Fitting  (Vd) 

Simple  Sum 
Perimeter 
Limit 

Low  Angle 

1.5  X  1.5  X  1.0 

20  20 

0.02 

Medium  Angle 

1.2  X  1.2  X  0.8 

20  20 

0.02 

High  Angle 

1.0  X  1.0  X  0.6 

10  10 

0.02 

Reflection  Statistics  (from  SORTAV) 

Total  Reflections 

86369 

Rejected  Outliers 

33 

Unique  Reflections 

29051 

Average  Redundancy 

3.0 

Resolution 

1.329  A'' 

Completeness 

91.9  % 

Ri 

5.77  % 

R2 

5.34  % 

Rw 

15.25  % 

Z  (Refinement) 

1.219 

Table  2. 


Selected  crystal,  integration,  and  reflection  data  for  17;0-estradiol»Vimethanol. 
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Figure  2.  Coordinate  system  for  the  methanol  molecule. 
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n 

m 

<n> 

Ri 

Rz 

/?W 

z 

y 

0<-4 

0 

0 

0.0 

0.0000 

0.0000 

0.0000 

0.000 

0.000 

-4<0<-3 

0 

0 

0.0 

0.0000 

0.0000 

0.0000 

0.000 

0.000 

-3  <  Q  <  -2 

8 

4 

2.0 

0.3321 

0.3844 

0.3018 

1.034 

0.263 

-2  <  Q  <  -1 

173 

74 

2.3 

0.5492 

0.6949 

0.3870 

0.761 

0.455 

■EEBBSi 

2508 

941 

2.7 

1.1243 

1.1755 

1.0245 

1.183 

2.153 

0<Q<1 

7959 

2832 

2.8 

1.0254 

1.0485 

0.9274 

1.329 

■reiTM 

1<Q<2 

7845 

2724 

2.9 

0.5759 

0.6902 

0.4618 

1.252 

esei 

2<Q<3 

6736 

2209 

3.0 

0.3661 

0.4556 

0.2863 

1.194 

3<Q<4 

5725 

1806 

3.2 

0.2602 

0.3270 

0.2132 

1.220 

0.235 

4<Q<6 

10058 

2909 

3.5 

0.1800 

0.2248 

0.1548 

1.231 

E|j^ 

6<Q<8 

8497 

2264 

3.8 

0.1266 

0.1558 

0.1169 

1.254 

8<Q<10 

6921 

1723 

4.0 

0.0969 

0.1194 

0.0927 

1.242 

10  <  Q  <  20 

18395 

4459 

4.1 

0.0596 

0.0679 

0.0645 

1.224 

20  <  Q  <  30 

4702 

1214 

3.9 

0.0389 

0.0677 

0.0370 

1.193 

30  <  0  <  50 

1174 

361 

3.3 

0.0232 

0.0273 

0.0241 

1.141 

0.025 

50<o<i6o 

212 

75 

2.8 

0.0115 

0.0136 

0.0142 

1.028 

0.011 

100<Q 

0 

0 

0.0 

0.0000 

0.0000 

0.0000 

0.000 

0.000 

Table  3.  Intensity-Significance  Intervals  where  n  is  the  number  of  reflections,  m  is  the 
number  of  unique  reflections,  <n>  is  the  average  measurement  multiplicity,  and  Q=I/Max 
((Tint/tText)  respectively  for  17y^estradiol»‘/imethanol. 


n 

m 

<n> 

Ri 

/?W 

z 

y 

D>  1.016 

4805 

1570 

3.1 

.0345 

.0529 

.0721 

1.568 

.035 

1.016  >D>  0.806 

6705 

1582 

mm 

.0512 

.0579 

.0943 

1.358 

.057 

0.806  >  D  >  0.705 

6464 

1594 

4.1 

.0554 

.0578 

.1004 

1.337 

.061 

0.705  >D >0.640 

7744 

1576 

4.9 

.0711 

.0719 

.1065 

1.271 

.075 

0.640  >D  >0.594 

8049 

1558 

5.2 

.0792 

.0777 

.1080 

1.250 

.081 

0.594  >D  >0.559 

7545 

1532 

4.9 

.1092 

.1085 

.1249 

1.189 

.105 

0.559  >D  >0.531 

6152 

1521 

4.0 

.1373 

.1371 

.1424 

1.155 

.124 

.531  >D  >0.508 

4376 

1441 

3.0 

.1334 

.1241 

.1677 

1.274 

.135 

.508  >D  >0.488 

4098 

1401 

2.9 

.1508 

.1428 

.1755 

1.245 

.148 

.488  >D  >0.472 

4118 

1436 

2.9 

.1593 

.1445 

.1895 

1.280 

.155 

.472  >D  >0.457 

3747 

1354 

2.8 

.1854 

.1681 

.1988 

1.230 

.177 

.457  >D  >0.444 

3746 

1387 

.2289 

.2205 

.2122 

1.192 

.444  >D  >0.432 

3456 

1331 

2.6 

.3003 

.3007 

.2497 

1.183 

.432  >D>  0.422 

3298 

1303 

2.5 

.3737 

.3823 

.2719 

1.167 

.333 

.422  >D>  0.412 

2886 

1194 

.3703 

.3773 

.2801 

1.166 

.321 

.412  >D  >0.403 

1376 

641 

2.1 

.3296 

.3234 

.3198 

1.396 

.309 

.403  >D  >0.395 

986 

493 

2.0 

.3404 

.3040 

.3670 

1.502 

.352 

.395  >D  >0.388 

754 

377 

2.0 

.4069 

.3446 

.4136 

1.574 

.431 

.388  >D  >0.381 

500 

250 

2.0 

.4569 

.4011 

.4303 

1.430 

.486 

.381  >D  >0.374 

108 

54 

2.0 

.5333 

.4996 

.5168 

1.417 

.584 

Table  4.  Equal-Volume  Resolution  Shells  where  n  is  the  number  of  reflections,  m  is  the 
number  of  unique  reflections,  <n>  is  the  average  measurement  multiplicity,  and 
S=sin6'/l  (A  ’)  respectively  for  17>^estradiol*‘/2methanol. 
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Monopole 

HIO 

0.40 

H20 

0.38 

HI 

0.23 

H2 

0.22 

H4 

0.26 

H6x 

0.20 

H7x 

0.17 

H8 

0.20 

H9 

0.16 

mix 

0.17 

H12x 

0.16 

H14 

0.19 

H15x 

0.16 

H16x 

0.18 

H17 

0.13 

H18x 

0.18 

Monooole 

o«2 

_ SP _ 

sp^ 

20 

33+ 

32- 

01 

-0.50 

02 

-0.49 

Cl 

-0.30 

-0.22 

0.34 

C2 

-0.38 

-0.19 

0.37 

C3 

0.27 

-0.21 

0.38 

C4 

-0.33 

-0.17 

0.36 

C5 

-0.18 

-0.22 

0.33 

C6 

-0.26 

0.31 

C7 

-0.31 

0.34 

C8 

-0.21 

0.39 

C9 

-0.17 

0.31 

CIO 

-0.25 

-0.18 

0.37 

Cll 

-0.31 

0.35 

C12 

-0.28 

0.31 

C13 

-0.16 

0.38 

C14 

-0.20 

0.38 

C15 

-0.26 

0.33 

C16 

-0.35 

0.42 

C17 

0.20 

0.38 

C18 

-0.32 

0.27 

Atoms 

Kappa 

K 

k’ 

01, 02, 03 

1 

0.97 

1.16 

C3 

2 

1.01 

0.92 

C17 

3 

1.02 

0.95 

Cl,  C2,  C4 

4 

0.97 

0.92 

C5,  CIO 

5 

0.98 

0.87 

C6,  C7,  C8,  C9, 
C11,C12,C13,C14, 
C15,C16,C17,C18,C19 

6 

0.98 

0.95 

all  C-H  hydrogen  atoms 

7 

1.20 

1.29 

HIO,  H20,  H30 

8 

1.20 

1.29 

Table  5.  Starting  values  entered  into  the  model  for  the  multipole  refinement  for  \lfi- 
estradiol»Vimethanol.  Units  for  multipole  populations  are  e’. 
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N 

0.9652(12) 

0.0311(12)1 

0.6260(11)1 

0.8125(11)1 

0.8194(10) 

0.6651(11)1 

0.5919(11)1 

0.5233(11) 

0.4609(11)1 

0.4195(11)1 

0.5028(10) 

0.4670(11)1 

0.4192(11) 

0.3270(11)1 

0.2685(11)1 

0.3477(11)1 

0.2797(12)  1 

0.2507(12) 

0.1119(11)1 

0.0791(12)1 

0.1875(11) 

0.1658(11) 

0.2845(11) 

0.1751(12) 

0.2510(18) 

0.0361(18) 

0.1644(16) 

0.2411(16) 

(91)91110 

-0.0323(15) 

0.1370(16) 

-0.1709(16) 

-0.0687(15) 

0.1410(16) 

-0.0908(15) 

0.1042(16) 

0.2269(16) 

-0.0845(16) 

0.0952(15) 

(91)68910- 

-0.1869(17) 

0.0144(17) 

-0.2360(16) 

-0.0401(17) 

-0.1975(15) 

0.2267(16) 

0.2577(16) 

0.1811(17) 

X 

-0.0503(22) 

0.1284(24) 

-0.2009(20) 

-0.2242(21) 

0.3944(20) 

0.5897(19) 

0.5127(20) 

0.4844(20) 

0.6269(20) 

0.3190(21) 

0.1338(19) 

-0.1825(20) 

-0.0617(20) 

-0.0363(20) 

-0.1555(20) 

0.3164(20) 

0.6267(22) 

0.5532(22) 

0.5101(21) 

0.4559(22) 

0.1966(20) 

0.0354(21) 

0.1336(20) 

0.2817(22) 

Atom 

0 

a 

0 

a 

HI 

H2 

H4 

V9H 

H6B 

H7A 

H7B 

00 

a 

a 

HllA 

aim 

H12A 

H12B 

H14 

H15A 

H15B 

a 

'O 

a 

H17 

H18A 

H18B 

H18C 

N 

0.936405(8) 

0.06557(7) 

00 

3 

SO 

d 

0.77285(7)  1 

0.82894(7) 

0.77447(7) 

0.66531(7)  1 

0.60930(8) 

0.50539(7)  1 

0.43507(7)  1 

0.49496(7) 

0.60917(7)  1 

0.42644(7) 

0.31502(7)  1 

0.25505(6)  1 

0.32756(7) 

0.25121(7)  1 

0.14061(7)  1 

0.15990(8)  1 

00 

0 

00 

f-H 

d 

0.21162(10) 

-0.05126(10) 

0.15097(10) 

0.19541(10) 

0.17723(11) 

0.12231(11) 

0.08148(10) 

0.03222(10) 

-0.05741(10) 

0.02383(10) 

0.02870(10) 

0.09174(11) 

0.10675(10) 

0.03067(11) 

0.02225(10) 

-0.05715(10) 

-0.08784(11) 

-0.11963(10) 

(01)16060  0- 

0.18231(10) 

X 

0.08240(13) 

0.10594(13) 

-0.06613(13) 

-0.08242(13) 

0.08629(14) 

0.26644(14) 

0.28180(13) 

0.48051(13) 

0.49100(12) 

0.32007(13) 

0.12862(13) 

0.11375(14) 

-0.04837(13) 

-0.03113(14) 

0.15891(13) 

0.33043(13) 

0.51074(13) 

0.43778(12) 

0.21795(13) 

0.15271(14) 

^H 

)R 

r- 

00 

0 

2 

in 

s 

00 

Atom 

0 

0 

u 

u 

U 

U 

U 

U 

U 

U 
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Table  7.  Fractional  atomic  coordinates  for  molecule  2  and  methanol  of  ITy^stradioFVimethanol. 


Atom 

u“ 

- iF - 

- iF - 

Ul3 

- JP - 

01 

0.02315(12) 

0.01838(11) 

0.01172(8) 

-0.00821(9) 

0.00038(8) 

-0.00143(7) 

02 

0.02742(14) 

0.01660(10) 

0.01237(8) 

-0.00827(10) 

-0.00257(8) 

-0.00055(7) 

Cl 

0.01183(10) 

0.02079(13) 

0.01240(9) 

■0.00540(9) 

0.00016(8) 

-0.00046(9) 

C2 

0.01489(11) 

0.01999(13) 

0.01241(10) 

-0.00567(10) 

0.00116(8) 

-0.00063(9) 

C3 

0.01715(12) 

0.01344(10) 

0.01130(9) 

-0.00612(9) 

-0.00016(8) 

0.00036(8) 

C4 

0.01523(11) 

0.01497(11) 

0.01220(9) 

-0.00595(9) 

-0.00180(8) 

0.00039(8) 

C5 

0.01212(10) 

0.01365(10) 

0.01193(9) 

-0.00438(8) 

-0.00169(8) 

0.00091(8) 

C6 

0.01108(10) 

0.02189(14) 

0.01563(11) 

-0.00391(10) 

-0.00226(8) 

-0.00091(10) 

C7 

0.01153(10) 

0.01874(13) 

0.01690(11) 

-0.00068(9) 

-0.00173(9) 

-0.00146(10) 

C8 

0.01108(9) 

0.01379(10) 

0.01239(9) 

-0.00381(8) 

-0.00023(7) 

-0.00045(8) 

C9 

0.01150(10) 

0.01434(11) 

0.01169(9) 

-0.00501(8) 

-0.00072(7) 

0.00044(8) 

CIO 

0.01091(9) 

0.01456(11) 

0.01085(9) 

-0.00475(8) 

-0.00064(7) 

0.00068(8) 

Cll 

0.01118(10) 

0.02348(15) 

0.01289(10) 

-0.00428(10) 

-0.00113(8) 

-0.00183(9) 

C12 

0.01433(11) 

0.02257(15) 

0.01333(10) 

-0.00799(10) 

-0.00128(8) 

-0.00153(9) 

C13 

0.01540(11) 

0.01279(10) 

0.01145(9) 

-0.00557(9) 

-0.00003(8) 

-0.00079(7) 

C14 

0.01351(10) 

0.01445(11) 

0.01285(9) 

-0.00400(9) 

0.00044(8) 

-0.00125(8) 

C15 

0.01524(13) 

0.02956(19) 

0.01735(13) 

-0.00556(13) 

0.00290(10) 

-0.00514(12) 

C16 

0.02131(15) 

0.02492(17) 

0.01594(12) 

-0.00703(13) 

0.00430(11) 

-0.00569(11) 

C17 

0.02088(13) 

0.01382(11) 

0.01227(10) 

-0.00648(15) 

-0.00006(11) 

-0.00162(13) 

C18 

0.02731(16) 

0.01397(12) 

0.01504(11) 

-0.00815(11) 

-0.00065(11) 

0.00034(9) 

or 

0.02050(12) 

0.02548(14) 

0.01264(9) 

-0.00112(10) 

0.00182(8) 

-0.00266(9) 

02’ 

0.02908(14) 

0.01675(10) 

0.01253(8) 

-0.00712(10) 

-0.00434(9) 

-0.00013(7) 

cr 

0.01201(11) 

0.02698(17) 

0.01402(11) 

0.00007(11) 

-0.00084(13) 

-0.00137(10) 

C2’ 

0.01366(12) 

0.02877(18) 

0.01463(11) 

-0.00027(11) 

0.00023(9) 

-0.00244(11) 

C3’ 

0.01544(12) 

0.01728(12) 

0.01231(10) 

-0.00199(10) 

0.00038(9) 

-0.00223(9) 

C4’ 

0.01434(11) 

0.01596(11) 

0.01173(9) 

-0.00340(9) 

-0.00136(8) 

-0.00047(8) 

C5’ 

0.01218(10) 

0.01417(10) 

0.01179(9) 

-0.00328(8) 

-0.00185(8) 

0.00034(8) 

C6’ 

0.01206(11) 

0.02589(16) 

0.01333(10) 

-0.00459(10) 

-0.00254(8) 

0.00307(10) 

C7’ 

0.01301(11) 

0.02461(15) 

0.01363(10) 

-0.00787(10) 

-0.00229(8) 

0.00273(10) 

C8’ 

0.01167(9) 

0.01277(10) 

0.01186(9) 

-0.00365(8) 

-0.00136(7) 

0.00172(7) 

C9’ 

0.01227(10) 

0.01343(10) 

0.01241(9) 

-0.00207(8) 

-0.00242(8) 

0.00063(8) 

CIO’ 

0.01146(10) 

0.01537(11) 

0.01218(9) 

-0.00169(8) 

-0.00187(8) 

-0.00036(8) 

Cll’ 

0.01272(11) 

0.02953(18) 

0.01384(11) 

-0.00590(12) 

-0.00282(9) 

0.00010(11) 

C12’ 

0.01563(12) 

0.02041(14) 

0.01340(10) 

-0.00364(11) 

-0.00403(10) 

0.00092(9) 

C13’ 

0.01718(11) 

0.01090(9) 

0.01195(9) 

-0.00572(9) 

-0.00183(8) 

0.00114(7) 

C14’ 

0.01539(11) 

0.01184('9) 

-0.00583(9) 

-0.00107(8)  1 

0.00163(8) 

C15’ 

0.01570(13) 

0.03738(24) 

0.01515(12) 

-0.00321(14) 

0.00057(10) 

0.00073(13) 

C16’ 

0.02051(15) 

0.03456(22) 

0.01457(12) 

-0.00670(15) 

0.00243(11) 

-0.00056(13) 

C17’ 

0.02185(14) 

0.01486(11) 

0.01165(9) 

-0.00722(10) 

-0.00124(9) 

0.00096(8) 

C18’ 

0.03705(22) 

0.01437(12) 

0.01709(12) 

-0.01337(14) 

-0.00299(13) 

0.00180(10) 

03 

0.02682(15) 

0.02440(15) 

0.02598(14) 

-0.00745(12) 

0.01130(12) 

-0.00554(12) 

C19 

0.02730(20) 

0.03127(23) 

0.02636(19) 

-0.01196(18) 

0.00339(15) 

0.00142(16) 

Table  8.  Anisotropic  thermal  parameters  of  non-H  atoms  for  17y5-estradiol*‘/2methanol. 
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Atom 

Ujso 

HIO’ 

0.0326(27) 

H20’ 

0.0316(26) 

HI’ 

0.0470(24) 

H2’ 

0.0495(26) 

H4’ 

0.0393(22) 

H6C 

0.0478(24) 

H6D 

0.0580(29) 

H7C 

0.0512(26) 

H7D 

0.0428(22) 

H8’ 

0.0422(23) 

H9’ 

0.0479(25) 

HI  1C 

0.0550(28) 

HI  ID 

0.0519(26) 

H12C 

0.0515(26) 

H12D 

0.0505(25) 

H14’ 

0.0450(23) 

H15C 

0.0578(29) 

H15D 

0.0600(30) 

H16C 

0.0636(31) 

H16D 

0.0657(31) 

H17’ 

0.0553(24) 

H18D 

0.0590(28) 

H18E 

0.0611(29) 

H18F 

0.0749(36) 

Atom 

Uiso 

HIO 

0.0320(26) 

H20 

0.0353(29) 

HI 

0.0453(25) 

H2 

0.0460(26) 

H4 

0.0429(24) 

H6A 

0.0404(21) 

H6B 

0.0491(26) 

H7A 

0.0495(25) 

H7B 

0.0475(24) 

H8 

0.0457(24) 

H9 

0.0386(21) 

HllA 

0.0488(25) 

HUB 

0.0493(25) 

H12A 

0.0519(26) 

H12B 

0.0462(24) 

H14 

0.0414(22) 

H15A 

0.0601(29) 

H15B 

0.0558(28) 

H16A 

0.0592(28) 

H16B 

0.0602(29) 

H17 

0.0505(23) 

H18A 

0.0572(27) 

H18B 

0.0554(26) 

H18C 

0.0658(31) 

H30 

0.0402(30) 

H19A 

0.0770(36) 

H19B 

0.0837(40) 

H19C 

0.0903(43) 

Table  9. 


Isotropic  thermal  parameters  of  H  atoms  for  17/^estradiol»'/2methanol. 


Atoms 

Bond  Length  (A) 

01’-C3’ 

1.3687(4) 

U 

1 

o 

1.4284(5) 

Cl’-C2’ 

1.3930(5) 

cr-cio’ 

1.4029(5) 

C2’-C3’ 

1.3944(5) 

C3’-C4’ 

1.3938(5) 

C4’-C5’ 

1.4007(4) 

C5’-C6’ 

1.5124(5) 

C5’-C10’ 

1.4073(4) 

C6’-C7’ 

1.5247(5) 

C7’-C8’ 

1.5270(5) 

C8’-C9’ 

1.5447(4) 

C8’  -  C14’ 

1.5225(4) 

C9’-C10’ 

1.5254(4) 

C9’-Cir 

1.5404(5) 

Cll’-C12’ 

1.5405(5) 

C12’-C13’ 

1.5304(5) 

C13’-C14’ 

1.5408(5) 

C13’-C17’ 

1.5421(5) 

C13’-C18’ 

1.5319(5) 

C14’-C15’ 

1.5351(5) 

C15’-C16’ 

1.5499(6) 

C16’-C17’ 

1.5511(6) 

Atoms 

Bond  Length  (A) 

01 -C3 

1.3688(4) 

02-C17 

1.4259(5) 

C1-C2 

1.3926(5) 

Cl -CIO 

1.4009(4) 

C2-C3 

1.3961(5) 

C3-C4 

1.3925(5) 

C4-C5 

1.3981(5) 

C5-C6 

1.5111(5) 

C5-C10 

1.4069(4) 

C6-C7 

1.5247(5) 

C7-C8 

1.5264(5) 

C8-C9 

1.5435(4) 

C8-C14 

1.5222(5) 

C9-C10 

1.5223(5) 

C9-C11 

1.5377(5) 

C11-C12 

1.5381(5) 

C12-C13 

1.5266(5) 

C13-C14 

1.5421(5) 

C13-C17 

1.5374(5) 

C13-C18 

1.5366(5) 

C14-C15 

1.5392(5) 

C15-C16 

1.5532(6) 

C16-C17 

1.5453(6) 

03-C19 


1.4234(8) 


Table  10.  Bond  distances  of  non-H  atoms  of  17y^estradiol«Vimethanol. 


11 


VO  ovoS  ^  vovo  vo^  os^  ^ 

Vw/  V--'  's*^  'W'  V-^  'w'  'w'  ^ ''w^  'w-'  'W'  'w^  'W'  'w' 

5|T^;mo^oo^ooo^OTt^ovoa^^cNo^cn^cnp^r^'^^^^po^vqTtpr^^ 

<J|rn  drowSuSvdrna\r^^<N^^r'5^G5ooo5r^^3;^r:'asoo^oja\t::;i^ 

^l-H  (NOOOOOOO-Hr-H:i:o^-H^oo-Ht^cDdppq--|-HoqJ:^ 


r>in2^a2SiSSiS"§';5SSS§2;22s--^§2 

CGoSaaussssaossSaSooOBasSxa! 

c^  I  I  i  I  I  i  I  I  I  I  I  >  I  I  I  I  i  I  t  I  I  I  ’  I  t  I  i 
^  I  {  I  I  I  I  ^  I  I  I  I  *  *  I  I  1  I  1  i  I  I  I  I 

ggggggg^. 


S  i2,  s  ij  >2 
u  u  u  u  g 


uuuug  ^ 


u  u  u 


<— i'OVO'0^<-<'-<'0'0'0-H-Hi— I«-IVOVOVO 

'w''  'w'  'W^  'w '  'W^  'W''  'w^  'w^  'w'  'w' 

Clr^  r-iinON^'^'opc^'^t^csosprar^cN 
<3lt^  ONOsOOCSrH^t^d^t^^pcj^t^OO 
p^l0000’-Hi-H^000*-^r4f<l^»^00 

^1^  ^t-HT-Hi-Ht-Hi-Hi-H»-Hi-Hi-H»-Hi-H»— I'pHi-Hi-H 


/-*s 

>«-n' 

00 

o 

8 

i™H 

t«H 

VOVO'00^^»^’-HO^  1-H 

Sw/  V,^  Sw-'  .  'W*' 

voc^c^io^<N<N^tn  *0 
odo\dodoN*od|^cn  d 


Pn  P'  MO  O'  P' 


88S088SS0055S 


^pnCpaSt^^-^w-^eQ^Ti-  00  00  00 

^-HC^pffioSSSSffiOUUUU  u 
Kffi® ^7  I ® ® I  I  1  I  I  I  I 

IIT..  .111111  I 


I  •  •  _  I  i 

On  r7  ^  00  00 

u  u  u  u 


U  u  ^  '-' 


m'->uuuuuQp'P'S2J:J 


K  «  ®  ffi 

I  I  I  I 
CM  ri  <s  (M 


^  t-h  cn  m 

T— (  1-H  1-H  ^ 

o  u  u  u 


<;  (S  (S  o. 

^  i-M  fH  f-H  i-H  1-H  ip-H 

2  u  u  u  u  u  u 

a 


00  ^  SO  VO  ^  NO  NO 

W-  'w/  'w/  s^  'w' 

lo  Tt  ^  Tt  1-H  Os  Os 
Qs  cvj  vd  Os  Os  d 


^^^voso^^»— '1— ivo 

'w<'  'w/  's^ 

^Tt<Nosov«npu^Tros 

o6osi-Ht^c5o6d^cnos 

t-Hi-HCS^CV|^CS<S^O 


^OS^SOSONOSOOO^  t-h 

w  'W'  Nw^  W'  W'  w  S--/ 

cvi^»nr^»nr^T};pco  co 
dvdddosr^ososov 
^Ot-H»— tOOOOO  ^ 


ligssoaasssssssOgGlIllismsEs 


sSgSGGOOoCOSSSOgoSsSSSSDc  d|6|d|9 


u  u  u  u 

I  I  I  I 


QgGS§Guuoouuuu55S^CCuu|^uuuu^^ 


Table  1 1.  Bond  angles  for  molecule  1  of  17;ff-estradiol*‘/6methanol. 
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C7’-C8’-C9’ _ 108.7(1)  C8’ -C14’ -C13’ _ 113.3(1)  03-C19-H19C _ 111.3(5) 

C7’-C8’-C14’ _ 112.9(1)  C8’-C14’-C15’ _ 119.6(1)  H19A  -  C19  -  H19B  107.4(8) 

C9’-C8’-C14’ _ 108.7(1)  C13’ -  C14’ -  C15’ _ 103.7(1)  H19A  -  C19  -  H19C  111.3(8) 

C7’-C8’-H8’  I  106.6(6)  ~1  C8’ -  C14’ -H14’  |  106.6(6)  |  H19B  -  C19  -  H19C  |  106.1(8) 

Table  12.  Bond  angles  for  molecule  2  and  methanol  of  17y^estradiol»‘/2methanol. 


Atom 

Monopole 
Population  (Po.o) 

or 

6.539(21) 

02’ 

6.543(21) 

cr 

4.194(36) 

C2’ 

4.256(35) 

C3’ 

3.875(34) 

C4’ 

4.295(34) 

C5’ 

4.123(34) 

C6’ 

4.222(35) 

CT 

4.199(34) 

C8’ 

4.115(35) 

C9’ 

4.115(33) 

CIO’ 

4.090(35) 

cir 

4.215(36) 

C12’ 

4.215(36) 

C13’ 

4.220(34) 

C14’ 

4.120(35) 

C15’ 

4.309(37) 

C16’ 

4.344(37) 

C17’ 

3.852(33) 

C18’ 

4.365(36) 

Atom 

Monopole 
Population  (Po.o) 

oi 

6.536(20) 

02 

6.533(20) 

Cl 

4.200(34) 

C2 

4.241(36) 

C3 

3.857(34) 

C4 

4.279(36) 

C5 

4.107(32) 

C6 

4.216(37) 

C7 

4.211(37) 

C8 

4.124(36) 

C9 

4.101(35) 

CIO 

4.123(36) 

Cll 

4.205(36) 

C12 

4.187(37) 

C13 

4.232(32) 

C14 

4.131(33) 

C15 

4.313(37) 

C16 

4.366(36) 

C17 

3.874(34) 

C18 

4.355(37) 

03 

6.491(20) 

C19 

4.264(35) 

Table  13. 


Monopole  populations  (e')  of  non-H  atoms  of  17/^estradiol»Vimethanol. 


Atom 

Monopole 
Population  {Po.o) 

HIO’ 

0.638(18) 

H20’ 

0.651(18) 

HI’ 

0.796(19) 

H2’ 

0.788(21) 

H4’ 

0.790(19) 

H6C 

0.832(15) 

H6D 

0.832(15) 

H7C 

0.842(14) 

H7D 

0.842(14) 

H8’ 

0.803(18) 

H9’ 

0.829(20) 

HllC 

0.827(15) 

HI  ID 

0.827(15) 

H12C 

0.841(14) 

H12D 

0.841(14) 

H14’ 

0.821(19) 

H15C 

0.841(15) 

H15D 

0.841(15) 

H16C 

0.870(15) 

H16D 

0.870(15) 

H17’ 

0.939(20) 

H18D 

0.879(13) 

H18E 

0.879(13) 

H18F 

0.879(13) 

Atom 

Monopole 
Population  {Po.o) 

HIO 

0.665(19) 

H20 

0.636(19) 

HI 

0.783(19) 

H2 

0.769(19) 

H4 

0.776(19) 

H6A 

0.826(14) 

H6B 

0.826(14) 

H7A 

0.842(14) 

H7B 

0.842(14) 

H8 

0.819(19) 

H9 

0.833(17) 

HllA 

0.833(14) 

HUB 

0.833(14) 

H12A 

0.837(14) 

H12B 

0.837(14) 

H14 

0.820(18) 

H15A 

0.857(15) 

H15B 

0.857(15) 

H16A 

0.872(14) 

H16B 

0.872(14) 

H17 

0.916(20) 

HI  8  A 

0.887(13) 

H18B 

0.887(13) 

H18C 

0.887(13) 

H30 

0.696(20) 

H19A 

0.850(12) 

H19B 

0.850(12) 

H19C 

0.850(12) 

Table  14.  Monopole  populations  (e')  of  H  atoms  of  17y^estradiol«14methanol. 
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Multipoles 

01 

or 

02 

02’ 

03 

Pi.^i 

-0.025(17) 

-0.033(17) 

-0.043(17) 

-0.028(18) 

-0.046(20) 

Pi.-i 

0.030(19) 

0.0 

0.0 

0.0 

0.0 

Pi.o 

0.0 

0.0 

0.0 

0.0 

-0.068(19) 

P2.0 

0.117(12) 

0.127(13) 

0.083(12) 

0.089(12) 

0.063(14) 

Pz^I 

-0.033(11) 

-0.016(12) 

-0.021(11) 

-0.018(11) 

-0.038(14) 

P2.-I 

-0.047(11) 

-0.027(12) 

0.0 

-0.021(12) 

-0.014(13) 

P2.+2 

-0.042(11) 

-0.035(11) 

-0.077(11) 

-0.046(11) 

0.0 

P2..2 

-0.017(11) 

0.0 

0.0 

0.0 

0.0 

P3.0 

-0.029(20) 

0.0 

0.0 

0.050(22) 

-0.040(25) 

P3M 

-0.029(19) 

0.0 

0.0 

0.0 

-0.026(24) 

P3..I 

-0.028(20) 

0.0 

-0.024(20) 

0.0 

0.027(21) 

P3.^2 

0.0 

0.0 

0.0 

0.0 

0.0 

P3..2 

0.0 

0.034(21) 

0.0 

0.0 

0.0 

P3.+3 

0.083(18) 

0.095(18) 

0.075(19) 

0.102(19) 

0.073(20) 

P3.-3 

-0.019(19) 

0.0 

0.0 

-0.038(19) 

-0.037(21) 

0.037(19) 

0.057(22) 

0.0 

0.025(20) 

0.0 

P4.+1 

0.0 

-0.035(20) 

0.0 

0.0 

0.0 

P4..1 

0.053(18) 

0.0 

0.025(18) 

0.040(18) 

0.0 

P4.+2 

0.0 

0.0 

0.0 

0.023(19) 

0.023(22) 

P4..2 

0.0 

0.0 

0.0 

0.021(19) 

0.067(23) 

P4.+3 

0.0 

0.0 

-0.039(18) 

0.0 

0.0 

P4..3 

0.0 

-0.022(18) 

0.020(18) 

0.0 

0.0 

P4.+4 

0.024(16) 

0.038(16) 

0.0 

0.0 

0.020(18) 

P4.-4 

0.0 

0.036(17) 

0.0 

-0.020(16) 

-0.030(20) 

Table  15.  Multipole  populations  (e‘)  of  Oxygen  atoms  of  17;ff-estradiol»l^methanol. 
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Table  16.  Multipole  populations  (e  )  of  Carbon  atoms  of  17;^stradiol*Vimethanol. 
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u 
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00 
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00 

00 
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-0.098(26) 

00 
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Multipoles 
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cC 
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C19 
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00 
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00 
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0.044(26) 

C18 

-0.028(25) 

0.034(24) 

-0.034(24) 

-0.017(17) 

-0.039(17) 

(91)9100- 

0.054(18) 

0.023(17) 

-0.028(29) 

00 

00 

0.039(27) 

0.212(26) 

00 

00 

u 

0.039(19) 

00 

00 

0.028(14) 

0.014(13) 

00 

o 

d 

-0.031(13) 

-0.066(21) 

-0.087(21) 

00 

00 

0.224(21) 

o 

d 

-0.021(20) 

C17 

0.028(20) 

00 

-0.067(18) 

-0.018(13) 

0.052(13) 

00 

00 

-0.078(13) 

-0.074(22) 

-0.080(21) 

O 

d 

0.074(22) 

0.212(21) 

-0.076(21) 
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d 

Co 

u 

-0.034(26) 

00 

-0.049(24) 

o 

d 

-0.025(19) 

-0.050(20) 

O 

d 

-0.030(18) 

0.059(30) 

O 

d 

0.022(29) 

-0.070(29) 

0.205(28) 

-0.116(28) 

-0.039(27) 

C16 

o 

d 

-0.031(25) 

-0.061(26) 

0.033(20) 

00 

o 

d 

00 

-0.019(18) 

(oe)iwo 

-0.091(29) 

0.034(28) 

O 

d 

0.274(28) 

-0.135(25) 

0.061(26) 

Multipoles 
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cC 

1 

1 

1 

1 

1 

1 

1 

1 

1 

Table  17.  Multipole  populations  (e')  of  Carbon  atoms  of  17;^stradiol*l/2methanol  continued. 
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o 

0.127(19) 

0.129(20) 

0.191(23) 

0.150(23) 

0.112(22) 

0.138(16) 

0.138(16) 
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0.127(24) 

0.139(16) 

0.139(16) 

0.140(16) 

0.140(16) 

0.176(24) 

0.075(17) 

0.075(17) 

0.094(17) 

0.094(17) 

0.219(23) 

0.102(13) 

0.102(13) 

0.102(13) 

(61)880*0 

0.183(13) 

0.183(13) 

0.183(13) 

Atoms 

b 

S 

b 

a 

X 

a 

s 

H6C 

H6D 

H7C 

H7D 

00 

X 

a 

HI  1C 

HI  ID 

H12C 

H12D 

rH 

X 

H15C 

H15D 

H16C 

H16D 

HIT 

H18D 

H18E 

H18F 

H30 

H19A 

H19B 

H19C 

O 

cC‘ 

0.005(28) 

/-N 

00 

8 

d 

-0.007(31) 

0.066(32) 

-0.018(30) 

(61)0100- 

-0.010(19) 

0.034(20) 

0.034(20) 

-0.029(31) 

0.028(29) 

-0.015(21) 

-0.015(21) 

0.019(20) 

(03)6100 

-0.006(29) 

0.018(22) 

0.018(22) 

-0.005(22) 

-0.005(22) 

0.041(33) 

-0.035(17) 

-0.035(17) 

-0.035(17) 

o 

<C 

0.121(19) 

0.125(21) 

0.109(24) 

0.107(22) 

0.127(23) 

0.105(15) 

0.105(15) 

0.134(15) 

0.134(15) 

0.150(23) 

0.089(22) 

0.110(16) 

0.110(16) 

(91)6600 

(91)6600 

0.118(24) 

0.070(17) 

0.070(17) 

0.132(16) 

0.132(16) 

0.167(24) 

0.106(13) 

0.106(13) 

0.106(13) 

Atoms 

OIH 

H20 

HI 

H2 

H4 

V9H 

H6B 

H7A 

H7B 

00 

X 

a 

VllH 

P3 
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< 
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(S 

X 

H14 

H15A 

PQ 

in 

X 

1-H 

X 

PQ 

VO 

1-H 

tl! 

H17 

V81H 

H18B 

H18C 

Table  18.  Multipole  populations  (e  )  of  Hydrogen  atoms  of  17>^estradiol*V^methanol. 


Bond _ jO(rc)  V  jO(rc) _ ^ _ ^ _ ^ _ ^ _ ^ _ ^ _ ^ 

01 -C3 _ 1.992  -18.666  1.3731  0.8342  0.5389  -17.30  -13.85 _ 1249 _ 0.25 

31’ -C3’ _ 2.101 _ -18.114  1.3692 _ 0.8106 _ 0.5586 _ -17.16 _ -15.42 _ 14.46 _ 031 

31-HlO  2.379  -31.187  0.9702  0.7400  0.2302  -37.16  -36.51  42.48  0.02 
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Figure  7.  The  Laplacian  of  the  total  electron  density  of  atoms  at  rest  in  the 

H12A  -  C12  -  H12B  plane  of  17;^estradiol«‘/2methanol.  Contour  intervals  are  5 
eA'^  starting  at  5  eA'^  (solid  blue  lines),  -2  eA'^  starting  at  -2  eA'^  (dotted  red 
lines),  and  the  dashed  line  plots  0  eA'^. 
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Figure  27. 


Geometry  of  hydrogen  bonding  interactions  of  molecule  1  of  llfi- 
estradiol*‘/imethanol. 
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Appendix  I. 


Crystal  structure  and  charge  density  analysis  of  estrone 


1 


C(12) 


Figure  la.  The  estrone  molecide  with  50%  probability  thermal  ellipsoids 
for  non-hydrogen  atoms. 


Experimental 


Compound  name: 

estrone 

Chemical  formula: 

C18H22O2 

Chemical  formula  weight: 

270.36 

Crystal  density: 

1.263 

Crystal  size,  mm 

0.2x0.25x0.30 

p,  mm'* 

0.05 

Unit  cell  parameters: 

a  (A) 

7.7618(1) 

b(A) 

9.9536(1) 

c(A) 

18.4055(1) 

V(A) 

1421.97(2) 

z 

4 

F(OOO) 

584 

Space  group 

P2i2,2i 

Temperature: 

92.0(2) 

Radiation(A): 

AgK«(A.=0.56086A) 

hmin/hmax 

-20/20 

kttimdCrnax 

-26/25 

Imin/lmax 

-48/48 

0inin/0max 

1.75/48.26 

Intensity  decay 

0% 

Total  No.  of  reflections 

90695 

No.  of  independent  reflections 

16278 

Rint 

0.0231 

No.  of  reflections  used 

8875 

Criterion 

I>3cy 

Refinement  on 

F^ 

Total  No.  of  parameters 

810 

Weighting  scheme 

w=l/cy^ 

Multipole  refinement: 

R(F2) 

0.0279 

Rw(F') 

0.0366 

S 

1.3337 

K-  refinement: 

R 

0.0445 

Rw 

0.0611 

S 

1.54 

4 


in  estrone 


Atomic  coordinates  and  Uiso 


Atom 

X 

y 

X 

Uiso 

0(1) 

-1.34702  (4) 

-0.80621(3) 

-0.31144  (2) 

0.022 

0(2) 

-0.81887(4) 

-0.28634  (3) 

0.15628  (2) 

0.021 

C(l) 

-1.27691  (5) 

-0.55423  (4) 

-0.16813  (2) 

0.017 

C(2) 

-1.35941(5) 

-0.62550  (4) 

-0.22323  (2) 

0.018 

C(3) 

-1.27233  (5) 

-0.73037(4) 

-0.25798  (2) 

0.017 

C(4) 

-1.10329  (5) 

-0.76064  (4) 

-0.23767  (2) 

0.017 

C(5) 

-1.02040  (5) 

-0.68873  (4) 

-0.18217  (2) 

0.015 

C(6) 

-0.83709(5) 

-0.72801(4) 

-0.16266  (2) 

0.017 

C(7) 

-0.74637  (5) 

-0.62131(4) 

-0.11661  (2) 

0.017 

C(8) 

-0.86413  (4) 

-0.57930  (3) 

-0.05415(2) 

0.014 

C(9) 

-1.02082  (5) 

-0.50316(4) 

-0.08568  (2) 

0.015 

C(10) 

-1.10798  (5) 

-0.58421(4) 

-0.14580  (2) 

0.015 

C(ll) 

-1.14594  (5) 

-0.45811  (4) 

-0.02538(2) 

0.019 

C(12) 

-1.05678  (5) 

-0.37607(4) 

0.03510  (2) 

0.019 

C(13) 

-0.90195  (5) 

-0.45307(4) 

0.06482  (2) 

0.015 

C(14) 

-0.77799(5) 

-0.48930  (4) 

0.00219(2) 

0.015 

C(15) 

-0.61274  (5) 

-0.53505(4) 

0.04134  (2) 

0.020 

C(16) 

-0.59941(5) 

-0.43575(4) 

0.10593  (2) 

0.021 

C(17) 

-0.78054  (5) 

-0.37816(4) 

0.11509(2) 

0.017 

C(18) 

-0.96009(6) 

-0.57650  (4) 

0.10964 (2) 

0.022 

H(1A) 

-1.4610 (9) 

-0.7709  (7) 

-0.3207  (3) 

0.038  (2) 

H(l) 

-1.3493  (8) 

-0.4749  (6) 

-0.1419(3) 

0.049(2) 

H(2) 

-1.4904 (9) 

-0.6010  (7) 

-0.2387  (3) 

0.045(2) 

H(4) 

-1.0358 (9) 

-0.8413 (7) 

-0.2651  (4) 

0.052 (2) 

H(6A) 

-0.8420 (9) 

-0.8227 (6) 

-0.1328(3) 

0.054  (2) 

H(6B) 

-0.7634  (8) 

-0.7498 (7) 

-0.2118  (3) 

0.052  (2) 

H(7A) 

-0.6238  (8) 

-0.6610 (6) 

-0.0970  (3) 

0.057(2) 

H(7B) 

-0.7152 (8) 

-0.5326  (6) 

-0.1490  (3) 

0.049(2) 

H(8) 

-0.9113 (8) 

-0.6713 (7) 

-0.0278(3) 

0.052  (2) 

H(9) 

-0.9648 (8) 

-0.4117  (6) 

-0.1094  (3) 

0.052  (2) 

H(llA) 

-1.2105 (8) 

-0.5456  (7) 

-0.0021(3) 

0.043 (2) 

H(llB) 

-1.2491(9) 

-0.3978  (7) 

-0.0491(4) 

0.054  (2) 

H(12A) 

-1.1501(8) 

-0.3501 (6) 

0.0772(3) 

0.049(2) 

H(12B) 

-1.0124  (8) 

-0.2810  (6) 

0.0121(3) 

0.047  (2) 

H(14) 

-0.7469(7) 

-0.3938  (7) 

-0.0248(3) 

0.048  (2) 

H(15A) 

-0.6236 (9) 

-0.6387 (7) 

0.0602(3) 

0.051(2) 

H(15B) 

-0.5003  (10) 

-0.5324  (7) 

0.0057(4) 

0.056  (2) 

H(16A) 

-0.5572 (9) 

-0.4816  (7) 

0.1568(3) 

0.052 (2) 

H(16B) 

-0.5107 (9) 

-0.3521 (7) 

0.0974(4) 

0.063  (2) 

H(18A) 

-1.0328  (8) 

-0.6461  (6) 

0.0787(3) 

0.057  (2) 

H(18B) 

-1.0397  (8) 

-0.5423  (7) 

0.1525(3) 

0.058 (2) 

H(18C) 

-0.8571(9) 

-0.6322  (7) 

0.1324(3) 

0.065  (2) 

5 


Uy  values 

Atom 

0(1) 

0(2) 

C(l) 

C(2) 

C(3) 
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Figure  5a.  The  dynamic  Fourier  multipole  electron  density  map 
in  the  0(2)-C(17)-H(l  A)next_moiecuie  plane.  The  contour  interval  is  0.05  eA 
The  hydrogen  bonding  region  is  shown. 


Figure  5b.  The  static  Fourier  multipole  electron  density  map 
in  the  0(2)-C(17)-H(lA)next_moiecuie  plane.  The  contour  interval  is  0.05  eA 
The  hydrogen  bonding  region  is  shown. 


Figure  5  c.  The  Laplacian  of  the  electron  density  map  in  the 
0(2)-C(  1 7)-H(  1  A)next_molecule  plane.  The  contour  interval  is  2,4,8  x  eA’^ 

plus  100  eA'^  contour.  The  hydrogen  bonding  region  is  shown. 
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The  (3,-1)  bond  critical  points  in  the  estrone  crystal. 
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Figure  6.  The  electrostatic  potential  of  the  estrone  molecule. 
Surfaces  are  +0.5  eA"'  (blue)  and -0.1 5  eA"'  (red). 
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The  multipole  population  parameters  for  hydrogen  atoms  in  estrone 
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Department  of  Chemistry,  University  of  Toledo, 
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The  scientific  £q>proach  to  prevendng  or  curing  disease  requires  understanding  the 
mechanism  by  which  the  (hsease  state  is  initiated  or  propogated.  A  significant  body  of 
knowledge  exists  to  indicate  that  a  large  number  of  breast  cancers  are  initiated  by  estrogens 
-  the  &mily  of  female  hormones.  It  has  been  hypodiesized  that  the  mechanism  involves 
first  binding  of  the  steroid  molecule  by  the  estrogen  recq>tor  in  the  nucleus  of  target  breast 
cells,  followed  by  secondary  interaction  which  partially  reorganizes  die  receptor  pocket. 
This,  in  turn,  modifies  the  surface  of  the  proteifi  such  diat  the  genes  necessary  for  tumor 
development  are  activated.  The  purpose  of  the  present  study  is  to  obtain  information  on  the 
fundamental  properties  of  estrogen  molecules  tl^  are  responsible  for  the  binding  event  and 
for  the  subsequent  reorganization  of  the  receptor. 

When  two  molecules  approach  one  another,  be  they  small  organic  molecules  in  solution  or  a 
steroid  molecule  approaching  a  protein,  they  sense  each  other’s  electrostatic  potential. 

Much  as  two  magnets  will  feel  each  other’s  magnetic  field  and  orient  diemselves  to 
maximize  the  attraction,  so  the  two  molecules  will  seek  to  align  themselves  in  such  a  way  as 
to  iq^proach  negative  regions  to  positive  regions  and  vice  versa.  Perfect  complementarity  of 
the  potentials  will  maximize  the  primary  binding  interaction  and  the  secondary 
reorganization  step  necessary  for  gene  expression.  Molecules  with  different  electrostatic 
potentials  will  have  different  binding  affinities  and  different  activities  with  respect  to  gene 
expression.  These  parameters  are  known  for  many  compoimds,  whereas  the  electrostatic 
potentials  are  not. 

The  electrostatic  potential  is  a  complicated  three  dimensional  function  that  is  determined  by 
the  geometrical  structure  of  a  molecule  and  the  distribution  of  the  electrons  and  nuclei  that 
make  up  that  molecule.  By  using  very  accurate  X-ray  diffiaction  techniques,  the  complete 
description  of  the  geometrical  structure  and  the  distribution  of  all  of  the  electrons  in  a 
molecule  may  be  determined.  This  information  may  then  be  used  to  calculate  the 
electrostatic  potential  at  any  point  in  or  around  a  molecule.  We  report  the  first  studies  of 
this  type  on  a  number  of  natural  and  unnatizral  estrogen  derivatives. 
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Electrostatic  potential  isosurfaces  for  the  C3-hydroxy  Electrostatic  potential  isosurfaces  for  the  C3-hydroxy  Electrostatic  potential  isosurfaces  for  the  C3-hydroxy  Electrostatic  potential  isosurfaces  for  the  C3~hydroxy 

group  (left)  and  C17-hydroxy  group  (right);  blue  1.0  eA  ^  group  (left)  and  C17-hydroxy  group  (right);  blue  0.75  eA  *.  group  (left)  and  Cl7“Carbonyl  group  (right);  group  (left)  and  C17-carbonyl  group  (right);  blue  0.75  eA  *. 

red -0.2  eA'*  red-0.4eA'‘  red -0.15  eA  *  red -0.2  eA’*  red -0.15  eA  ,  blue  0.5  eA  .  red -0.075  eA*^  red -0.35  eA'^ 


Appendix  K. 

Electrostatic  potential  isosurfaces  for  several  hydroxyl  groups  of  selected  estrogens 
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Figure  1.  Electrostatic  potential  isosurfaces  of  the  01  hydroxy  group  of  17;^estradiol  • 
urea.  The  blue  surface  in  each  picture  represents +1.0eA'^  The  red  surfaces 
represent  the  following  potentials:  A  -0.05  eA'\  B  -0.10  eA'\  C  -0.15  eA'\  D  - 
0.20  eA-’. 
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Figure  2.  Electrostatic  potential  isosurfaces  of  the  02  hydroxy  group  of  17;ff-estradiol  • 
urea.  The  blue  surface  in  each  picture  represents  +1 .0  eA'\  The  red  surfaces 
represent  the  following  potentials:  A  -0.05  eA‘\  B  -0.10  eA‘\  C  -0.15  eA'\  D  - 
0.20  eA'^  E  -0.25  eA’',  F  -0.30  eA'\  G-0.35  eA‘*. 
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Figure  3.  Electrostatic  potential  isosurfaces  of  the  01  hydroxy  group  of  ITy^estradiol  • 

methanol.  The  blue  surface  in  each  pictare  represents  +1.0  eA"'.  The  red  surfaces 
represent  the  following  potentials:  A  -0.05  eA'\  B  -0.10  eA'\  C  -0.15  eA'\  D  - 
0.20  eA-’. 
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Figure  4.  Electrostatic  potential  isosurfaces  of  the  01  ’  hydroxy  group  of  17y0-estradiol  • 

methanol.  The  blue  surface  in  each  picture  represents  +1.0  eA‘^  The  red  surfaces 
represent  the  following  potentials:  A  -0.05  eA'\  B  -0.10  eA'\  C  -0.15  eA'',  D  - 
0.20  eA'^ 
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Figure  5.  Electrostatic  potential  isosurfaces  of  the  02  hydroxy  group  of  17>9-estradiol  • 

methanol.  The  blue  surface  in  each  picture  represents  +1.0  eA'*.  The  red  surfaces 
represent  the  following  potentials:  A  -0.05  eA‘\  B  -0.10  eA"',  C  -0.15  eA‘\  D  - 
0.20  eA‘‘. 
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Figure  6.  Electrostatic  potential  isosurfaces  of  the  02’  hydroxy  group  of  17>0-estradiol  • 

methanol.  The  blue  surface  in  each  picture  represents  +1.0  eA‘^  The  red  surfaces 
represent  the  following  potentials:  A  -0.05  eA'*,  B  -0.10  eA‘',  C  -0.15  eA"',  D  - 
0.20  eA-’. 


Figure  7.  Electrostatic  potential  isosurfaces  of  the  01  hydroxy  group  of  17a-estradiol  •  ¥2 
H2O.  The  blue  surface  in  each  picture  represents  +1.0  eA  '.  The  red  surfaces 
represent  the  following  potentials:  A  -0.05  eA‘\  B  -0.10  eA‘\  C  -0.15  eA'\  D  - 
0.20  eA'^ 
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)tential  isosurfaces  of  the  02  hydroxy  group  of  17a-estradiol  •  Vi 
surface  in  each  picture  represents  +1.0  eA"’.  The  red  surfaces 
)llowing  potentials:  A  -0.05  eA'*,  B  -0.10  eA’',  C  -0.15  eA'*,  D  - 
.25  eA'*,  F  -0.30  eA'’,  G  -0.35  eA'*,  H  -0.40  eA'\  I  -0.45  eA’^ 


